Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,332,8230AG56.7% 54.3 / 107.6 93pseudogene (1580/4605 nt)yfaSpseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (12/27);  new base G (51/0);  total (66/27)
Fisher's exact test for biased strand distribution p-value = 5.68e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-10

CTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACC  >  NC_000913/2332676‑2332969
                                                                                                                                                    |                                                                                                                                                   
cTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGAt                                                                                                                                                  >  2:1763116/1‑151 (MQ=255)
  aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGttt                                                                                                                                                >  1:584296/1‑151 (MQ=255)
  aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATTGGCGGG‑GGGGGGGttt                                                                                                                                                >  1:428001/1‑151 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCTGCGGC‑GGGGCGAtttt                                                                                                                                               >  1:1039601/1‑150 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGtttt                                                                                                                                               >  2:1516695/1‑150 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGtttt                                                                                                                                               >  1:1558325/1‑150 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGTTTa                                                                                                                                               >  2:1340567/1‑151 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTa                                                                                                                                               >  1:868766/1‑151 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGGGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGC‑GGGGGGAGTTa                                                                                                                                               >  2:1451057/1‑151 (MQ=255)
     aCTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGC‑GGGGGGAtttta                                                                                                                                              >  2:887130/1‑148 (MQ=255)
        gTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGttttaagc                                                                                                                                           >  1:331578/1‑145 (MQ=255)
        gTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGGTTTAACGc                                                                                                                                           >  2:1250206/1‑150 (MQ=255)
        gTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGc                                                                                                                                           <  2:331686/150‑1 (MQ=255)
         tATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGcc                                                                                                                                          <  1:1112637/150‑1 (MQ=255)
          aTTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGCGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGTTTTAAGCca                                                                                                                                         >  1:882521/1‑150 (MQ=255)
          aTTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGAGTTTAAGCca                                                                                                                                         >  1:1062006/1‑150 (MQ=255)
             gtgCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGTTTTAAGCCACAtt                                                                                                                                     >  1:1306779/1‑150 (MQ=255)
             gtgCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGC‑GGGGGGGTTTTACGCcacaag                                                                                                                                     >  1:1298904/1‑149 (MQ=255)
              tgCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATg                                                                                                                                     <  2:197657/150‑1 (MQ=255)
                cGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGAc                                                                                                                                   <  2:491334/150‑1 (MQ=255)
                cGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGAc                                                                                                                                   <  2:1559425/150‑1 (MQ=255)
                  gCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTNNNNNNGGCGGGGGNNATTTACCGCCCNAtacct                                                                                                                                  >  1:594743/1‑146 (MQ=255)
                  gCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAACGCCACAAGAccc                                                                                                                                 >  2:160962/1‑149 (MQ=255)
                     gTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTTACGCCACATGACCCCTgg                                                                                                                             >  2:885530/1‑151 (MQ=255)
                      tACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAACGCCCCCTGaaccccgg                                                                                                                             >  2:273283/1‑143 (MQ=255)
                       aCTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTTTCGGCGGC‑GGGGGGATTTAACGCCACATGACCGCTgggg                                                                                                                           >  2:592729/1‑149 (MQ=255)
                         ttAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGTTTAAAGCCACATGACCGCTggggg                                                                                                                          >  2:1372256/1‑147 (MQ=255)
                              gAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAACGCCCCATGACCCCTGGCGTCGGgca                                                                                                                    >  2:945136/1‑151 (MQ=255)
                                gTAATGCGCCAGCCCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGTGTATCGCCACCAGAACGCTGGCGGCGGgcac                                                                                                                   >  1:1611480/1‑150 (MQ=255)
                                gTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAACGCCACATGACCCCTGGCGGCGGacac                                                                                                                   >  2:1632823/1‑150 (MQ=255)
                                gTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATAACCGCTGGCGCCGGgcaa                                                                                                                   >  2:1117586/1‑149 (MQ=255)
                                gTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGC‑GGGGGGATTTAACGCCACATGACCCCTGGcggcgggaac                                                                                                                   >  2:653093/1‑146 (MQ=255)
                                 tAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGTTTTACACCACATGACCGCCTGCGTCGGGCACCg                                                                                                                 >  1:36928/1‑151 (MQ=255)
                                 tAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGacac                                                                                                                   <  1:1294602/149‑1 (MQ=255)
                                   aTGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGC‑GGGGGGGTTTTAAGCCACAAGACCGCCGGGGTCGGGCAccc                                                                                                                 >  1:961306/1‑148 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGTTTAAAGCCAAAAGACCGCTTGGGGCGGGCACCCtttt                                                                                                             >  2:584371/1‑150 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGGTTAACGCCCCATGACCGCTGGCGGCGGGCACCCtttt                                                                                                             >  2:382964/1‑150 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAACGCCACATGACCCCCGGGGGCGGGCAcccgtt                                                                                                              >  2:736908/1‑146 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAAAGCCACATGACCGCCGGCGTCGGACAcccgttt                                                                                                             >  2:918600/1‑146 (MQ=255)
                                      cgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGTTTTACGCCACCTGACCGCTGGGGTCGGGCAcccgttc                                                                                                             >  1:812992/1‑145 (MQ=255)
                                               gcgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGTTTTACGCCACAATACCGCCGGCGGCGGGCAACGGTTTGCTGccctt                                                                                                    >  1:1034570/1‑147 (MQ=255)
                                                cggtCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTcc                                                                                                   <  1:664636/147‑1 (MQ=255)
                                                cgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAACGCCACATGACCGCCGGCGTCGGGCACCGTTTCGCTTCCtttc                                                                                                   >  1:1175795/1‑148 (MQ=255)
                                                cgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTcc                                                                                                   <  2:1693207/150‑1 (MQ=255)
                                                cgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTcc                                                                                                   <  1:243436/150‑1 (MQ=255)
                                                       tGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGacac                                                                                           <  1:886643/151‑1 (MQ=255)
                                                       tGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGacac                                                                                           <  1:1597847/151‑1 (MQ=255)
                                                        gAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGGGGC‑GGGGGGGTTTTAGGCCACATGACCCCTGGCGGCGGGCACCGCTTCGCTGCCCTCCCACGACACg                                                                                          >  1:20420/1‑151 (MQ=255)
                                                        gAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGC‑GGGGGGGTTTTACGCCCCATGACCGCTTGCGTCGGGCCACCTTTCGCCTCCCTCCAGCggaacg                                                                                          >  2:1382936/1‑146 (MQ=255)
                                                          ggATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTTACGCCACATGACCGCCGGCGGCGGGCACCCTTTCGCTGCCTTCCCGCGACACggg                                                                                        >  2:790873/1‑151 (MQ=255)
                                                             tctcTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTTACGCCCCAAGACCCCCGGCGGCGGGCACCCTTTTGCTGCCCTCCCGCGACACggggga                                                                                     >  2:842580/1‑149 (MQ=255)
                                                              ctctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTTAAGCCACATGACCCCCGGGGTCGGGAACCCTTTCGCTGCCTTCCAGCGACACGGGGgaa                                                                                    >  2:431167/1‑151 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGTTTTACGCCACATGACCCCTGGCGTCGGGCACCCGTTTGCTGCCTTCCCGCGAACCGGGGGAACg                                                                                  >  2:574561/1‑151 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTTACGCCCCAAGACCGCCGGCGTCGGGCACCCTTTTGCCGCCCTCCCGCGACACGGGGgaaag                                                                                  >  2:1441194/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAACGCCACATTAACCCCGGCGGCGGGAACCCGTTCGCTGCCTTCCCGCGACACGGGGgaaag                                                                                  >  2:873167/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAACGCCACAAGACCGCTGGCGGCGGACACCCGTTCGCCGCCTTCCAGCGACACGGGGgaaag                                                                                  >  2:1550884/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGC‑GGGGAGATTTAAAGCCACATGACCGCTGGCGTCGGACACCGGTTCGCTGCCTTCCAGCGACACGGGGgaaag                                                                                  >  2:824817/1‑149 (MQ=255)
                                                                ctTTGGGGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCCGCCGC‑GGGGGGATTTTAGCCCACGTAATCCCCGGCGTTTGTGACCGTTTGGCTTGCTTTCCGCCGCACGGGGGAAca                                                                                  >  2:801983/1‑150 (MQ=255)
                                                                  ttGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGATCGGTTTATCGGCGGC‑GGGGGGATTTTACGCCCCATGACCGCCGGGCTCGGGCACCCTTTTGCTGCCTTCCCGCGAACCGGGGGAACGt                                                                                 >  2:1339285/1‑150 (MQ=255)
                                                                   tggtggtggTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTcc                                                                               <  2:821642/151‑1 (MQ=255)
                                                                        gtggtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGGGTTAACGCCACATGACCGCTGGCGTCGGGCACCGGTTCGCTGCCTTTCAGCGACACCGGGGAACGGCCGGCaa                                                                          >  1:66987/1‑150 (MQ=255)
                                                                              aCGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAACGCCACATGACCGCTGGGGTCGGGCACCGTTTTGCTGCCTTCCCGCGACACGGGGGAACGGCCGGCCTAACCCa                                                                    >  2:982354/1‑151 (MQ=255)
                                                                              aCGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGc                                                                     <  2:1391430/150‑1 (MQ=255)
                                                                              aCGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGc                                                                     <  1:104063/150‑1 (MQ=255)
                                                                                   aTACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTTACGCCCCATGACCGCTGGGGTCGGGCACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGGCCGGCATAAccccaaacg                                                               >  2:226168/1‑143 (MQ=255)
                                                                                        ccgcgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTAc                                                           <  1:157552/150‑1 (MQ=255)
                                                                                        ccgcgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTAc                                                           <  1:1643656/150‑1 (MQ=255)
                                                                                         cgcgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCg                                                         <  1:1737063/151‑1 (MQ=255)
                                                                                           cgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGAc                                                       <  1:247921/151‑1 (MQ=255)
                                                                                                  aCTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGAt                                               <  2:945666/151‑1 (MQ=255)
                                                                                                  aCTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGAt                                               <  2:458283/151‑1 (MQ=255)
                                                                                                  aCTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGAt                                               <  2:1527691/151‑1 (MQ=255)
                                                                                                   cTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGATa                                              <  2:1043375/151‑1 (MQ=255)
                                                                                                   cTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGATa                                              <  2:1393426/151‑1 (MQ=255)
                                                                                                     ggCGGTGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAACGCCACATGACCGCTGGCGGCGGACCCCGTTTCGCTGCCTTTCAGCGAAACGGGGAAACGGCCGGCATAAAGCAAAACGCTACCGAACCT‑TGTTAAt                                            >  2:843398/1‑151 (MQ=255)
                                                                                                         gagaCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGATAATTGc                                         <  2:753707/150‑1 (MQ=255)
                                                                                                         gagaCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGATAATTGc                                         <  1:1440543/150‑1 (MQ=255)
                                                                                                           gaCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGAACGCTGGCGTCGGACACCGGTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCCGCATAACGCAAATCGTTACCGACCCT‑TTATAATTGCTGa                                      >  1:1093701/1‑151 (MQ=255)
                                                                                                             cGGTTAACAAATAGCGGCTCGGTTTATCGGCGGG‑GGGGGGAATTAAAGCCAAATCAACGCCGCGGGCGGAAACCCGTTTGCTGGATTCCCGCGGCAACGGGTAACCTCCGGCATAACGCAAATCGGTACTGACCATATGA‑AAATGCTGAg                                     >  1:1423584/1‑150 (MQ=255)
                                                                                                             cGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCCACCAT‑TGATAATTGCTGAgg                                    >  2:1454250/1‑150 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGG‑GGGGAGATTTTACGCCACAAGACCGCGGGCGTCGGACACCGTTTCGCTTGCTTTCAGCGGCAAGGGGGAAAATCCGGCATAAAGGCAATAGGTACCGACCCT‑TGATTATTGCTGAgcc                                   >  2:1062972/1‑149 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAACGCCACATGGCCGCTGGCGTCCGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAAAGCAAATCGTTACCGACCCT‑TGATAATTGCTGAgcg                                   >  2:203618/1‑150 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGGCAAGGGGAAACGTCCCGCATAAAGCAAATCCTTACCGACCCT‑TGGTAATTTCCGAggg                                   >  2:251025/1‑148 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGGTTTGCTGCCTTTCAGCGACACGGGGAAAAGGCCGGCATAACGCAAATCGTTACCGACCAT‑TTATAATTTCTGAgcg                                   >  1:1599420/1‑150 (MQ=255)
                                                                                                                    cAAATAGCGGCTCGGTTTATCGGCGGT‑GGGGCGATTTAACCCCACATTCCCGGTCGGGGCCGACACCGTTTTGCCGCCTTCCAGCGGCACGGGGACACGGTCCGCATAACGCAAAATCTTACCGCCCAA‑TTATAATTGCTCAGCGCCCaaa                             >  1:1764624/1‑151 (MQ=255)
                                                                                                                    cAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGGGATTTTACGCCACATGACCGCTGGCCGCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAAAGTCCGGCATAACGCAAATCGTTACCGACCCT‑TGATAATTGGTGAGGGCGCcaa                             >  2:351071/1‑151 (MQ=255)
                                                                                                                             gCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGATAATTGCTGAGCGCGCAAACTTAACTg                     <  2:1096099/150‑1 (MQ=255)
                                                                                                                                  gTTTATCGGCGGC‑GGGGCGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGGAACGTCCGGCATAACGCAAATCGTTACCCACCAT‑TTATAATTGCTGAGCGCGCAAACTTAACTGCAcgcg               >  1:1348247/1‑151 (MQ=255)
                                                                                                                                  gTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACCCGGGGAAAAGTCCGACATAACGCAAATCGTTACCGACCAT‑TGATAATTGCTGAGCGCGCAAAATTAACTGCAcgcg               >  2:1746054/1‑151 (MQ=255)
                                                                                                                                  gTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGATAATTGCTGAGCGCGCAAACTTAACTGCAcgcg               >  1:1704896/1‑151 (MQ=255)
                                                                                                                                  gTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGATAATTGCTGAGCGCGCAAACTTAACTGCAcgcg               >  1:1445733/1‑151 (MQ=255)
                                                                                                                                          gcggc‑ggGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCAtttt       >  1:700375/1‑151 (MQ=255)
                                                                                                                                            ggc‑ggGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCAtttttt     >  1:896926/1‑151 (MQ=255)
                                                                                                                                                 gggAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTAcc  <  2:67013/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGC‑GGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCAT‑TGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACC  >  NC_000913/2332676‑2332969

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