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

CTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACC  >  NC_000913/2332676‑2332969
                                                                                                                                                    |                                                                                                                                                   
cTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAt                                                                                                                                                  >  2:1763116/1‑151 (MQ=255)
  aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGttt                                                                                                                                                >  1:584296/1‑151 (MQ=255)
  aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATTGGCGGGGGGGGGGttt                                                                                                                                                >  1:428001/1‑151 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCTGCGGCGGGGCGAtttt                                                                                                                                               >  1:1039601/1‑150 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGtttt                                                                                                                                               >  2:1516695/1‑150 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGtttt                                                                                                                                               >  1:1558325/1‑150 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTa                                                                                                                                               >  2:1340567/1‑151 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTa                                                                                                                                               >  1:868766/1‑151 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGGGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGAGTTa                                                                                                                                               >  2:1451057/1‑151 (MQ=255)
     aCTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGAtttta                                                                                                                                              >  2:887130/1‑148 (MQ=255)
        gTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGttttaagc                                                                                                                                           >  1:331578/1‑145 (MQ=255)
        gTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGGTTTAACGc                                                                                                                                           >  2:1250206/1‑150 (MQ=255)
        gTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGc                                                                                                                                           <  2:331686/150‑1 (MQ=255)
         tATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGcc                                                                                                                                          <  1:1112637/150‑1 (MQ=255)
          aTTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGCGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTAAGCca                                                                                                                                         >  1:882521/1‑150 (MQ=255)
          aTTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAGTTTAAGCca                                                                                                                                         >  1:1062006/1‑150 (MQ=255)
             gtgCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTAAGCCACAtt                                                                                                                                     >  1:1306779/1‑150 (MQ=255)
             gtgCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGGTTTTACGCcacaag                                                                                                                                     >  1:1298904/1‑149 (MQ=255)
              tgCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATg                                                                                                                                     <  2:197657/150‑1 (MQ=255)
                cGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGAc                                                                                                                                   <  2:491334/150‑1 (MQ=255)
                cGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGAc                                                                                                                                   <  2:1559425/150‑1 (MQ=255)
                  gCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTNNNNNNGGCGGGGGNNATTTACCGCCCNAtacct                                                                                                                                  >  1:594743/1‑146 (MQ=255)
                  gCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACAAGAccc                                                                                                                                 >  2:160962/1‑149 (MQ=255)
                     gTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTACGCCACATGACCCCTgg                                                                                                                             >  2:885530/1‑151 (MQ=255)
                      tACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCCCCTGaaccccgg                                                                                                                             >  2:273283/1‑143 (MQ=255)
                       aCTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTTTCGGCGGCGGGGGGATTTAACGCCACATGACCGCTgggg                                                                                                                           >  2:592729/1‑149 (MQ=255)
                         ttAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTAAAGCCACATGACCGCTggggg                                                                                                                          >  2:1372256/1‑147 (MQ=255)
                              gAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCCCATGACCCCTGGCGTCGGgca                                                                                                                    >  2:945136/1‑151 (MQ=255)
                                gTAATGCGCCAGCCCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTGTATCGCCACCAGAACGCTGGCGGCGGgcac                                                                                                                   >  1:1611480/1‑150 (MQ=255)
                                gTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCCCTGGCGGCGGacac                                                                                                                   >  2:1632823/1‑150 (MQ=255)
                                gTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATAACCGCTGGCGCCGGgcaa                                                                                                                   >  2:1117586/1‑149 (MQ=255)
                                gTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGATTTAACGCCACATGACCCCTGGcggcgggaac                                                                                                                   >  2:653093/1‑146 (MQ=255)
                                 tAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTACACCACATGACCGCCTGCGTCGGGCACCg                                                                                                                 >  1:36928/1‑151 (MQ=255)
                                 tAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGacac                                                                                                                   <  1:1294602/149‑1 (MQ=255)
                                   aTGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGGTTTTAAGCCACAAGACCGCCGGGGTCGGGCAccc                                                                                                                 >  1:961306/1‑148 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTAAAGCCAAAAGACCGCTTGGGGCGGGCACCCtttt                                                                                                             >  2:584371/1‑150 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGGTTAACGCCCCATGACCGCTGGCGGCGGGCACCCtttt                                                                                                             >  2:382964/1‑150 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCCCCGGGGGCGGGCAcccgtt                                                                                                              >  2:736908/1‑146 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAAAGCCACATGACCGCCGGCGTCGGACAcccgttt                                                                                                             >  2:918600/1‑146 (MQ=255)
                                      cgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTACGCCACCTGACCGCTGGGGTCGGGCAcccgttc                                                                                                             >  1:812992/1‑145 (MQ=255)
                                               gcgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTACGCCACAATACCGCCGGCGGCGGGCAACGGTTTGCTGccctt                                                                                                    >  1:1034570/1‑147 (MQ=255)
                                                cggtCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTcc                                                                                                   <  1:664636/147‑1 (MQ=255)
                                                cgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCCGGCGTCGGGCACCGTTTCGCTTCCtttc                                                                                                   >  1:1175795/1‑148 (MQ=255)
                                                cgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTcc                                                                                                   <  2:1693207/150‑1 (MQ=255)
                                                cgTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTcc                                                                                                   <  1:243436/150‑1 (MQ=255)
                                                       tGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGacac                                                                                           <  1:886643/151‑1 (MQ=255)
                                                       tGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGacac                                                                                           <  1:1597847/151‑1 (MQ=255)
                                                        gAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGGGGCGGGGGGGTTTTAGGCCACATGACCCCTGGCGGCGGGCACCGCTTCGCTGCCCTCCCACGACACg                                                                                          >  1:20420/1‑151 (MQ=255)
                                                        gAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGGTTTTACGCCCCATGACCGCTTGCGTCGGGCCACCTTTCGCCTCCCTCCAGCggaacg                                                                                          >  2:1382936/1‑146 (MQ=255)
                                                          ggATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTACGCCACATGACCGCCGGCGGCGGGCACCCTTTCGCTGCCTTCCCGCGACACggg                                                                                        >  2:790873/1‑151 (MQ=255)
                                                             tctcTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTACGCCCCAAGACCCCCGGCGGCGGGCACCCTTTTGCTGCCCTCCCGCGACACggggga                                                                                     >  2:842580/1‑149 (MQ=255)
                                                              ctctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTAAGCCACATGACCCCCGGGGTCGGGAACCCTTTCGCTGCCTTCCAGCGACACGGGGgaa                                                                                    >  2:431167/1‑151 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTACGCCACATGACCCCTGGCGTCGGGCACCCGTTTGCTGCCTTCCCGCGAACCGGGGGAACg                                                                                  >  2:574561/1‑151 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTACGCCCCAAGACCGCCGGCGTCGGGCACCCTTTTGCCGCCCTCCCGCGACACGGGGgaaag                                                                                  >  2:1441194/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATTAACCCCGGCGGCGGGAACCCGTTCGCTGCCTTCCCGCGACACGGGGgaaag                                                                                  >  2:873167/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACAAGACCGCTGGCGGCGGACACCCGTTCGCCGCCTTCCAGCGACACGGGGgaaag                                                                                  >  2:1550884/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGAGATTTAAAGCCACATGACCGCTGGCGTCGGACACCGGTTCGCTGCCTTCCAGCGACACGGGGgaaag                                                                                  >  2:824817/1‑149 (MQ=255)
                                                                ctTTGGGGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCCGCCGCGGGGGGATTTTAGCCCACGTAATCCCCGGCGTTTGTGACCGTTTGGCTTGCTTTCCGCCGCACGGGGGAAca                                                                                  >  2:801983/1‑150 (MQ=255)
                                                                  ttGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGATCGGTTTATCGGCGGCGGGGGGATTTTACGCCCCATGACCGCCGGGCTCGGGCACCCTTTTGCTGCCTTCCCGCGAACCGGGGGAACGt                                                                                 >  2:1339285/1‑150 (MQ=255)
                                                                   tggtggtggTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTcc                                                                               <  2:821642/151‑1 (MQ=255)
                                                                        gtggtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGGTTAACGCCACATGACCGCTGGCGTCGGGCACCGGTTCGCTGCCTTTCAGCGACACCGGGGAACGGCCGGCaa                                                                          >  1:66987/1‑150 (MQ=255)
                                                                              aCGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGGGTCGGGCACCGTTTTGCTGCCTTCCCGCGACACGGGGGAACGGCCGGCCTAACCCa                                                                    >  2:982354/1‑151 (MQ=255)
                                                                              aCGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGc                                                                     <  2:1391430/150‑1 (MQ=255)
                                                                              aCGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGc                                                                     <  1:104063/150‑1 (MQ=255)
                                                                                   aTACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTACGCCCCATGACCGCTGGGGTCGGGCACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGGCCGGCATAAccccaaacg                                                               >  2:226168/1‑143 (MQ=255)
                                                                                        ccgcgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTAc                                                           <  1:157552/150‑1 (MQ=255)
                                                                                        ccgcgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTAc                                                           <  1:1643656/150‑1 (MQ=255)
                                                                                         cgcgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCg                                                         <  1:1737063/151‑1 (MQ=255)
                                                                                           cgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGAc                                                       <  1:247921/151‑1 (MQ=255)
                                                                                                  aCTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGAt                                               <  2:945666/151‑1 (MQ=255)
                                                                                                  aCTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGAt                                               <  2:458283/151‑1 (MQ=255)
                                                                                                  aCTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGAt                                               <  2:1527691/151‑1 (MQ=255)
                                                                                                   cTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATa                                              <  2:1043375/151‑1 (MQ=255)
                                                                                                   cTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATa                                              <  2:1393426/151‑1 (MQ=255)
                                                                                                     ggCGGTGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGCGGCGGACCCCGTTTCGCTGCCTTTCAGCGAAACGGGGAAACGGCCGGCATAAAGCAAAACGCTACCGAACCTTGTTAAt                                            >  2:843398/1‑151 (MQ=255)
                                                                                                         gagaCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGc                                         <  2:753707/150‑1 (MQ=255)
                                                                                                         gagaCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGc                                         <  1:1440543/150‑1 (MQ=255)
                                                                                                           gaCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGAACGCTGGCGTCGGACACCGGTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCCGCATAACGCAAATCGTTACCGACCCTTTATAATTGCTGa                                      >  1:1093701/1‑151 (MQ=255)
                                                                                                             cGGTTAACAAATAGCGGCTCGGTTTATCGGCGGGGGGGGGAATTAAAGCCAAATCAACGCCGCGGGCGGAAACCCGTTTGCTGGATTCCCGCGGCAACGGGTAACCTCCGGCATAACGCAAATCGGTACTGACCATATGAAAATGCTGAg                                     >  1:1423584/1‑150 (MQ=255)
                                                                                                             cGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCCACCATTGATAATTGCTGAgg                                    >  2:1454250/1‑150 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGGGGGGAGATTTTACGCCACAAGACCGCGGGCGTCGGACACCGTTTCGCTTGCTTTCAGCGGCAAGGGGGAAAATCCGGCATAAAGGCAATAGGTACCGACCCTTGATTATTGCTGAgcc                                   >  2:1062972/1‑149 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGGCCGCTGGCGTCCGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAAAGCAAATCGTTACCGACCCTTGATAATTGCTGAgcg                                   >  2:203618/1‑150 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGGCAAGGGGAAACGTCCCGCATAAAGCAAATCCTTACCGACCCTTGGTAATTTCCGAggg                                   >  2:251025/1‑148 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGGTTTGCTGCCTTTCAGCGACACGGGGAAAAGGCCGGCATAACGCAAATCGTTACCGACCATTTATAATTTCTGAgcg                                   >  1:1599420/1‑150 (MQ=255)
                                                                                                                    cAAATAGCGGCTCGGTTTATCGGCGGTGGGGCGATTTAACCCCACATTCCCGGTCGGGGCCGACACCGTTTTGCCGCCTTCCAGCGGCACGGGGACACGGTCCGCATAACGCAAAATCTTACCGCCCAATTATAATTGCTCAGCGCCCaaa                             >  1:1764624/1‑151 (MQ=255)
                                                                                                                    cAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTACGCCACATGACCGCTGGCCGCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAAAGTCCGGCATAACGCAAATCGTTACCGACCCTTGATAATTGGTGAGGGCGCcaa                             >  2:351071/1‑151 (MQ=255)
                                                                                                                             gCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTg                     <  2:1096099/150‑1 (MQ=255)
                                                                                                                                  gTTTATCGGCGGCGGGGCGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGGAACGTCCGGCATAACGCAAATCGTTACCCACCATTTATAATTGCTGAGCGCGCAAACTTAACTGCAcgcg               >  1:1348247/1‑151 (MQ=255)
                                                                                                                                  gTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACCCGGGGAAAAGTCCGACATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAAATTAACTGCAcgcg               >  2:1746054/1‑151 (MQ=255)
                                                                                                                                  gTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCAcgcg               >  1:1704896/1‑151 (MQ=255)
                                                                                                                                  gTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCAcgcg               >  1:1445733/1‑151 (MQ=255)
                                                                                                                                          gcggcggGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCAtttt       >  1:700375/1‑151 (MQ=255)
                                                                                                                                            ggcggGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCAtttttt     >  1:896926/1‑151 (MQ=255)
                                                                                                                                                 gggAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTAcc  <  2:67013/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACC  >  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: