Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,292,4580TA39.8% 114.8 / 42.8 83intergenic (+93/+51)galU/hnsglucose‑1‑phosphate uridylyltransferase/global DNA‑binding transcriptional dual regulator H‑NS
Reads supporting (aligned to +/- strand):  ref base T (37/13);  new base A (0/33);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 6.97e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTT  >  NC_000913/1292309‑1292602
                                                                                                                                                        |                                                                                                                                                
acCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAAt                                                                                                                                                  >  2:609931/1‑150 (MQ=255)
   cTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATaaa                                                                                                                                               >  1:1151387/1‑150 (MQ=255)
     tGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATaaaaa                                                                                                                                             <  2:251970/150‑1 (MQ=255)
      ggCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAAt                                                                                                                                           <  1:1019508/151‑1 (MQ=255)
        cACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATcc                                                                                                                                         <  2:1132309/151‑1 (MQ=255)
           ggAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCccgc                                                                                                                                      >  1:87729/1‑151 (MQ=255)
              aTTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCccgccg                                                                                                                                    >  1:443151/1‑150 (MQ=255)
               tttAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCt                                                                                                                                  >  1:279161/1‑151 (MQ=255)
                    aGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCgg                                                                                                                             >  2:43868/1‑151 (MQ=255)
                          gCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGAtttt                                                                                                                       >  2:203620/1‑151 (MQ=255)
                            ttGAAGAAGAGAGGGGCATAAAGAAGTAACATCCGTATCGGTGTTACCCACGAAACGGCGTTGAGCAACCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAACCCGGAAAATAAAAAATCCCGCCGCTGGCGGGATTTTa                                                                                                                      <  1:205440/150‑1 (MQ=255)
                              aaagaagAGAGGGCCATTAAATAGCACCCTCCGTATGTGTGTTACCCACGAACCGGCTTTGACCAACCGACGCTTTTTTTTTATATCTTATTCTTATTAAGTATTTTAAAACCGGACAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAg                                                                                                                    <  2:84131/149‑1 (MQ=255)
                               aagaagAGATGGGCATTAGAAAGTAACATCCTTATGGGTTTTACCCACAAAACGGTTTGGAGCAATCGACGCTTTTTTTTTATAGCTTTTTCTTATTAAATTTTTTAAACCCGGAAAAAAAAAAACCCCGCCGCTGGCGGGATTTTaa                                                                                                                     <  2:472135/148‑1 (MQ=255)
                                     agaTGGGAATAAAAAAGTAACACCCGTATCGGTTTTACCCAAAAAACGGCTTTGAGAAATCGACGCCTTTTTTTTATAGCTTATTTTTATTAAATTTTTTAAACCCGGAAAATAAAAACCCCCGCCGCTGGCGGGATTTTAAGCAAGTg                                                                                                              <  2:441119/149‑1 (MQ=255)
                                        tGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGGGCaa                                                                                                           >  2:768889/1‑149 (MQ=255)
                                         gggAATAAAGAATTAACATCCGTATCGGTTTTACCCACGAAACGGCTTTGAGCAATCGACCCTTTTTTTTTATAGTTTATTCTTTTTAATTTTTTTAAACCCGGACAATAAAAACCCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCt                                                                                                        <  2:1139026/151‑1 (MQ=255)
                                          ggCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCt                                                                                                        >  2:97854/1‑150 (MQ=255)
                                            cATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTAc                                                                                                      >  2:754171/1‑150 (MQ=255)
                                               tAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaaa                                                                                                   >  1:908023/1‑150 (MQ=255)
                                                aagaagTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACaaa                                                                                                   >  2:623647/1‑149 (MQ=255)
                                                 agaagTAACATCCGTATCGGTGTTATCCACGAAACGGCTTTGAGCAATCGGCGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTAAACCGGGAAAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGa                                                                                                <  1:971555/151‑1 (MQ=255)
                                                 aaaagTACCATCCTTATCGGTTTTTCCCACGAACCGGCGTGGAGCAATCGACGCCTTTTTTTTATAGCTTTTTCTTATAAAATTTTATAAACCCGGACAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGa                                                                                                <  2:656051/149‑1 (MQ=255)
                                                    agTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGatt                                                                                              >  1:876807/1‑150 (MQ=255)
                                                        acacCCGTCTGGGTTTTATCCACGAACCGGAGTGGAGCAATCGACGCCGTTTTTTTATATTTTATTCTTATAAAATTTTTTAAACCAGAACAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTg                                                                                          <  2:81642/150‑1 (MQ=255)
                                                           tCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCtt                                                                                       >  2:683506/1‑150 (MQ=255)
                                                           tCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAATGCAATCTACAAAAGATTATTGCtt                                                                                       >  2:864541/1‑150 (MQ=255)
                                                             cGTTTCGGTTTTTTCATCGAAATGGGTTTGAGCATCCGTCCCTTTTTTTATATACTTTTTCTTTATAAATTTCTCTTACCCAGAAAAAAAAAAAACACCGCAGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGAt                                                                                    <  2:685561/151‑1 (MQ=255)
                                                             cGTATCGGTTTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGAACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGa                                                                                     <  2:1110820/150‑1 (MQ=255)
                                                                aTCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGAt                                                                                    >  1:112338/1‑148 (MQ=255)
                                                                  gggTTTTACCCACAAAACGGCTTTGACCAACCGCCGCTGTTTTTTTATATTTTTTTTTTATAAAATTTTTTAAACCGGAAAAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGAGTATTGCTTGATCAGGa                                                                               <  2:1242183/150‑1 (MQ=255)
                                                                   ggTGTTATCCACGAAACGGCGTTGAGCAATCGCCGCTTTTTTTTTATAGCTTTTTCTTTTTAAATTGTTTAAACCCGGAAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGAt                                                                                    <  1:578020/145‑1 (MQ=255)
                                                                   ggTGTTATCCACAAAACGGCGTTGAGCATTCGACGCGGTTTTTTTATAGCTTATTCTTATTAAATTGTCTAAACCCGGAAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGaa                                                                              <  1:377741/151‑1 (MQ=255)
                                                                   ggTGTTACCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTTTCTAAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAgg                                                                                <  2:611912/149‑1 (MQ=255)
                                                                   ggTGTGATCCAAGAAGGGGCGTTGGCAAATGGGTGTGGTTTTTATATATTTTATTCTTATAAAATTTTTTAAAACGGGAAAAAAAAAATTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAgg                                                                                <  1:279060/149‑1 (MQ=255)
                                                                         aTCCACGAAACGGCGTTGAGCAATCGACGCCTTTTTTTTATAGTTTATTCTTATTAAATTGTCTTAAACCGAACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgt                                                                         <  1:196614/150‑1 (MQ=255)
                                                                            cACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgtcg                                                                       >  2:49709/1‑149 (MQ=255)
                                                                             acaaaaggggtttGAGAAATCGACCCTTTTTTTTTATACTTTTTTTTTATAAATTTTTCTAAACCCGAACAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcg                                                                          <  2:55420/134‑1 (MQ=255)
                                                                               gAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAggg                                                                   >  1:628798/1‑150 (MQ=255)
                                                                                aacgggggTTGAGCAATCGCCGCCTTTTTTTTATAGTTTATTCTTATAAAATTTTCTAAACCCGGAAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGa                                                                      <  2:962564/139‑1 (MQ=255)
                                                                                aaccGGCTTTGAGAAATCGACGCCTTTTTTTTATAGTTTTTTTTTATTAAATTTTTTTAAACGGGAAAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGAt                                                                 <  1:865524/148‑1 (MQ=255)
                                                                                aaaCGGCGTTGAGCAATCGCCGCCGTTTTTTTATAGCTTATTTTTATAAAATTTTCTAAACCCGAAAAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAAtcgtcg                                                                       <  2:244077/145‑1 (MQ=255)
                                                                                     gCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTAcc                                                            >  1:641643/1‑151 (MQ=255)
                                                                                        ttGAGCAATCGACGCTGTTTTTTTATAGTTTATTTTTATAAAATTGTCTAAACCCGAAAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCtt                                                          <  1:637573/150‑1 (MQ=255)
                                                                                        ttGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTg                                                         >  1:65100/1‑151 (MQ=255)
                                                                                        ttGAGCAATCGACGCCGTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGc                                                        >  1:328090/1‑151 (MQ=255)
                                                                                        ttGAGAAACCGACGCCTTTTTTTTATAGCTTATTTTTATTAATTTTTTTAAACCCGAAAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCtt                                                          <  2:865180/150‑1 (MQ=255)
                                                                                          gAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCt                                                       >  2:1151118/1‑151 (MQ=255)
                                                                                           aGCAACCCACGCGTTTTTTTTATACCTTATTTTTATAAAATTTTCTAAACCCGAAAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTc                                                      <  2:934871/151‑1 (MQ=255)
                                                                                              aaTCGCCGCTTTTTTTTTATAGTTTATTTTTATTAATTGGTCTTAAACCGGAAAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTc                                                      <  1:544407/148‑1 (MQ=255)
                                                                                                tcgacgaatttttttcATAGCTTTTTCTTATGAATTTTTTTAAAGCCGGAAAAAAAAAACTCCCGCCGCTGGCGGGATTTTAAGCATGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATcc                                                  <  2:12503/141‑1 (MQ=255)
                                                                                                 cGACGCCTTTTTTTTATAGCTTATTTTTATTAAATTGTCTAAAACCGAAAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCAt                                                <  1:1221068/151‑1 (MQ=255)
                                                                                                  gaggctttttttttATATTTTTTTTTTATTAATTTGTTTAAAACCGGAAAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCa                                                 <  2:400762/142‑1 (MQ=255)
                                                                                                     gCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCt                                            >  1:101537/1‑151 (MQ=255)
                                                                                                      ccGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTTAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCCTGCTCATCCATTGCtt                                           >  2:1188223/1‑151 (MQ=255)
                                                                                                      ccGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtt                                           >  2:1076084/1‑151 (MQ=255)
                                                                                                      ccGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAAATTATTGGTTGAACAGGGAATTGACCAGGGCTTTACCTTGGTCATCCATTGc                                             >  1:616475/1‑149 (MQ=255)
                                                                                                        gTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttt                                         >  2:113021/1‑151 (MQ=255)
                                                                                                         tttttttATTGATTATTTTTATAAATTTGTTTAAAACCGTAAAATAAAAATTCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTg                                              <  1:458512/145‑1 (MQ=255)
                                                                                                         tttttttATATCTTATTTTTATAAAATTGTTTAAACCCGGAAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtt                                           <  1:1132812/148‑1 (MQ=255)
                                                                                                         tttttttATAGTTTATTTTTATAAAATTGTCTAAACCCGGACAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTg                                              <  1:1176222/145‑1 (MQ=255)
                                                                                                           tttttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCtttttt                                       >  2:655513/1‑150 (MQ=255)
                                                                                                            ttttATATCTTATTTTTATTAAATTGTCTAAACCCGAAAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTg                                      <  2:366912/150‑1 (MQ=255)
                                                                                                             tttATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGa                                     >  1:350766/1‑150 (MQ=255)
                                                                                                             tttAAACCTTTTTTTTTTAAAATTTTCCTAAGCCGAACAAAAAAACATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTg                                      <  2:1024975/149‑1 (MQ=255)
                                                                                                              ttATAGTTTATTTTTTTTAATTTTTCTAAACCGGAAAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGAt                                    <  1:63243/150‑1 (MQ=255)
                                                                                                              ttATAGCTTTTTTTTATAAAATTGTCTAAAACCGGACAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGAt                                    <  2:423021/150‑1 (MQ=255)
                                                                                                                  agtttattcttattAAATTGTCTAAACCCGAAAAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACa                                <  1:1024364/147‑1 (MQ=255)
                                                                                                                  aGCTTATTTTTATTAAATTGTCTTAAACCGGACAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAg                               <  1:827712/151‑1 (MQ=255)
                                                                                                                     ttattcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGCCGAGGGATTTACCTTGCTTATCCATTGCTTTTTTGATTCCAGCt                             >  2:1017107/1‑150 (MQ=255)
                                                                                                                      ttttctttattaAATTGTTTAAACCCGGAAAATAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGc                              <  2:323984/143‑1 (MQ=255)
                                                                                                                        ttcttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCACTGGCGGGATTTTAAGCAAGGGCAGTCTACAAAAAATTATTTCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTTATCCATTTCTTTTTTTATTACAGCTgg                           >  2:42812/1‑149 (MQ=255)
                                                                                                                          cttattAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAg                         >  2:481285/1‑149 (MQ=255)
                                                                                                                            ttattaATTTGTCTTAAACCGGAAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAg                         <  2:885649/148‑1 (MQ=255)
                                                                                                                             tataaAATTGTCTAAACCCGAACAAAAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTAc                      <  2:1064144/150‑1 (MQ=255)
                                                                                                                              attaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACgg                    >  1:805878/1‑151 (MQ=255)
                                                                                                                              attaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACgg                    >  1:476932/1‑151 (MQ=255)
                                                                                                                               ttaAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGc                   >  2:88255/1‑151 (MQ=255)
                                                                                                                                  aaTTGTTTAAAACCGAAAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCtt                <  2:353047/151‑1 (MQ=255)
                                                                                                                                  aaTTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCtt                >  2:873295/1‑151 (MQ=255)
                                                                                                                                  aaTTGCATTAAACCGGAAAAAAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCt                 <  1:500931/150‑1 (MQ=255)
                                                                                                                                   ccttttctAAACCGATAAAAAAAAAATACCGACGCTGGTGGGATTTTAAGCAAGTGCAATCTACAAAAGAATATTGCTGGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGAGTACAGCGGGAGTTGGGCCTTg               <  2:1115198/144‑1 (MQ=255)
                                                                                                                                      gTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGcc            >  2:275387/1‑151 (MQ=255)
                                                                                                                                              ccGGACAATAAAAAACCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAg     <  1:99662/150‑1 (MQ=255)
                                                                                                                                                 tgACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCCTTGCTTTTTTGATTACAGCTGGAGTAAGGCCTTGGCCAGTCCAggct  >  2:90209/2‑149 (MQ=255)
                                                                                                                                                        |                                                                                                                                                
ACCCTTGGCACGGAATTTAAAGCCTGGCTTGAAGAAGAGATGGGCATTAAGAAGTAACATCCGTATCGGTGTTATCCACGAAACGGCGTTGAGCAATCGACGCCGTTTTTTTATAGCTTATTCTTATTAAATTGTCTTAAACCGGACAATAAAAAATCCCGCCGCTGGCGGGATTTTAAGCAAGTGCAATCTACAAAAGATTATTGCTTGATCAGGAAATCGTCGAGGGATTTACCTTGCTCATCCATTGCTTTTTTGATTACAGCTGGAGTACGGCCTTGGCCAGTCCAGGTT  >  NC_000913/1292309‑1292602

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: