Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,332,8230AG54.7% 50.1 / 91.8 77pseudogene (1580/4605 nt)yfaSpseudogene, bacterial alpha2‑macroglobulin YfaS variant family; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (12/22);  new base G (41/0);  total (55/22)
Fisher's exact test for biased strand distribution p-value = 1.06e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.37e-07

CTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACC  >  NC_000913/2332676‑2332969
                                                                                                                                                   |                                                                                                                                                    
cTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAt                                                                                                                                                  >  2:1052138/1‑151 (MQ=255)
cTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAt                                                                                                                                                  >  1:699030/1‑151 (MQ=255)
  aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAttt                                                                                                                                                >  2:1478728/1‑151 (MQ=255)
  aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAttt                                                                                                                                                >  2:1061487/1‑151 (MQ=255)
  aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCCGTTTATCGGCGGCCGGGGGGttt                                                                                                                                                >  2:275676/1‑151 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTa                                                                                                                                               >  2:1407700/1‑151 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGGtttt                                                                                                                                               >  1:278559/1‑150 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTAAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCCGGGGGGtttt                                                                                                                                               >  2:541462/1‑150 (MQ=255)
   atACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGCGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGGTAACAAAAAGCGGCTCGGTTTATCGGCGGGGGGGGGGTTTa                                                                                                                                               >  1:1293085/2‑151 (MQ=255)
     aCTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTaa                                                                                                                                              >  2:1353663/1‑150 (MQ=255)
     aCTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAttttaa                                                                                                                                             >  2:1082144/1‑148 (MQ=255)
     aCTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAtttta                                                                                                                                              >  1:869430/1‑148 (MQ=255)
     aCTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGCCGGCGGGGTGATTTaaa                                                                                                                                             >  1:1172564/1‑150 (MQ=255)
       tGTATTGTGCGGCAGTACTTAATGAGTAATGCGGCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGc                                                                                                                                           <  1:159995/151‑1 (MQ=255)
           ttGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAGTTTAAGCccca                                                                                                                                       >  1:1431566/1‑151 (MQ=255)
                cGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGAc                                                                                                                                   <  1:1073418/150‑1 (MQ=255)
                      tACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTACGCCACATGACCGCCgg                                                                                                                             >  2:1332561/1‑150 (MQ=255)
                      tACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCCCATGAcccccgg                                                                                                                             >  2:1047809/1‑145 (MQ=255)
                      tACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTg                                                                                                                              <  1:630239/149‑1 (MQ=255)
                       aCTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGGTTTTAAGCCCCATTAcccccgggg                                                                                                                           >  2:1363310/1‑144 (MQ=255)
                       aCTTAAGGAGTAATGCGCCAGACCGCGTGCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCg                                                                                                                           <  1:1230483/151‑1 (MQ=255)
                           aaTGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGCGTCgg                                                                                                                       >  2:1288673/1‑151 (MQ=255)
                             tGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGac                                                                                                                     <  2:374163/151‑1 (MQ=255)
                             tGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGac                                                                                                                     <  1:236991/151‑1 (MQ=255)
                              gAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTAAGCCACATGACCCCCGGCGTCGGgca                                                                                                                    >  2:1059307/1‑151 (MQ=255)
                                 tAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGacac                                                                                                                  <  1:934579/149‑1 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTACGCCCCATAAACCCTGGCGCCCGGCAcc                                                                                                                 >  1:240165/1‑146 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTTACGCCACATGACCGCTGGCGGCGGGCAACGtttt                                                                                                            >  2:1246081/1‑150 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCCCATGACCGCTGGCGGCGGGCACCGgttt                                                                                                            >  1:1054124/1‑150 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGCGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCCCGGTTTATCGGCGGGGGGGACGTTTAAACCCACATGAACCCCGGCGTccgggacccgttt                                                                                                             >  2:1421812/1‑139 (MQ=255)
                                       gcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCCCATGACCGCCGGCGTCGGGCACCCttttg                                                                                                           >  2:321035/1‑148 (MQ=255)
                                            aCCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGGTAACAAATAGCGGCTCGGTTTATCGGCGGGGGGGGGAGTTTAACACACAAGACAGCCGGCGTCGGACACCCGTTCGCTGcc                                                                                                      >  1:1343252/1‑150 (MQ=255)
                                                 gTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAAGAACGGCTCGGTTTATGGGCGGCGGGGGGGATTAAAGCCACATGACCGCTGGCGTCGGAAACCCGTTCGCTGCCTTCCAg                                                                                                >  1:841367/1‑151 (MQ=255)
                                                      aTGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCCCCGGCGTCGGACACCGTTTTGCTGCCCTCCCGCGac                                                                                            >  2:1479343/1‑150 (MQ=255)
                                                           gATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACgg                                                                                        <  1:13193/149‑1 (MQ=255)
                                                           gATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGATTTAAAGCCACATTACCGCTGGGGGCGGGCCCCGGTTCGCTGCCTTTCAGCGACACgggg                                                                                      >  1:1144001/1‑151 (MQ=255)
                                                             tctcTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCCCCTTACCCCCGGCGTCGGGCACCCTTTCGCTGCCTTCCCGCGACACggggga                                                                                    >  2:906766/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTAACCCCCCATAGCCCCCGGCGTCGGACCCCGCTTTGCGTCCTTTCAGCGACACGGGGgaaag                                                                                 >  1:1008885/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTAAAGGCACATGACCGCCGGGGTCGGACACCGTTTCGCTGCCTTTCACCCACACGGGGaaaag                                                                                 >  2:1236811/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTTACGCCACATGACCCCTCGCGTCGGACACCGTTTTGCTGCCCTCCCGCGACACGGGGaaaag                                                                                 >  2:846823/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCCCATGACCGCCGGCGTCGGAAACCGTTTTGCTGCCTTCCAGCGACACGGGGAAACg                                                                                 >  2:86085/1‑151 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCCGCCGGGCGATTTTATGCCCCATAACCGCCGGCGCCGCGAACCGTGTTGCTGCCTTCCACCCGAACCGGGAAACg                                                                                 >  1:286336/1‑151 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAAAAGCGGCCCGGTTTTTCGGCGGCGGGGGGATTTAACGCCCACTGACCCCTGCCGCCGGAAACCCGTTTCGTGCCTTCCCACAACCCGGGGaaaag                                                                                 >  2:647922/1‑149 (MQ=255)
                                                                         tggtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCAt                                                                         <  2:20058/150‑1 (MQ=255)
                                                                          ggtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGAATTTACGCCCCAAGACCGCTGGGGTGCGACACCGCTTCGCCGCCCTTCAGCGACACCGGGAAACGTCCGGCCTaa                                                                       >  2:1180647/1‑151 (MQ=255)
                                                                               cgcgATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAAAGCCCCATGACCGCTGGCGTCGGACACCCTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCaa                                                                  >  2:1242838/1‑151 (MQ=255)
                                                                               cgcgATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGAGTTTACACCACATTAACCCTGGCGGCGGGCCCCGGTTTGGTGGCTTTCCGGGAAACGGGGGAAAGTTCgggccaacgcca                                                                  >  1:129233/1‑141 (MQ=255)
                                                                                gcgATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGCGTCGGACACCGGTTCGCTGCCTTCCAGCGACACGGGGGAACGGCCGGCATAACGCaaa                                                                 >  2:367263/1‑151 (MQ=255)
                                                                                gcgATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCCCATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGGAACGGCCCGCATAACGCaaa                                                                 >  2:913727/1‑151 (MQ=255)
                                                                                 cgATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGGTTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTTGCTGCCTTCCAGCGGCACCGGGGAACGTCCGGCATAACGCAAAt                                                                >  1:63246/1‑151 (MQ=255)
                                                                                   aTACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCg                                                              <  2:914326/151‑1 (MQ=255)
                                                                                           cgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGAc                                                      <  2:869901/151‑1 (MQ=255)
                                                                                             aggTCACTGGCGAAGACGATTAACAAAGAGCGGCTCAGGTTATCGGCGGTGGGAAGATTTAACGCCAAATAACCGCTGGCGTCTGACACCGTTTCGCTGCCTTCCAACGACACGGGGAAACGTCCGGCAAAACGCAAATCGTTACCGGcc                                                     <  1:750447/148‑1 (MQ=255)
                                                                                                 cACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAACGACACGGGGAAAAGGCCGGCATAACGCAAATCGTTACCGACCATTGa                                                >  2:157062/1‑151 (MQ=255)
                                                                                                   cTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATa                                              <  2:843488/151‑1 (MQ=255)
                                                                                                         gagaCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGc                                         <  2:63279/150‑1 (MQ=255)
                                                                                                            aCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACCCGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCCTTGATAATTGCTGAg                                     >  2:1465933/1‑151 (MQ=255)
                                                                                                             aggTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATAACCGCTGGCGTCGGGAACCGTTTAGCTGCCTTCCAGGGAAAAGGGGAAACGTCCGGCATAAAGCCAATCGCTTACGAACATTGGTAATTGGTGAg                                     >  2:480853/2‑150 (MQ=255)
                                                                                                              ggTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAAAGCCACATTACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGGCCGGCATAACGCAAATCGTTTCCGACCCTTGATAATTGCTGAgcg                                   >  1:624545/1‑151 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGTGTTTTAAGCCACAAAACCGCTGGGGCGGGGCTCCGGTTTGCTGGCTTCCAGCCACACGGGGGAACGTTCGGGCTAACGCAAATCGTTTCCCCCCACTGAAAATTgtctagcgg                                  >  1:301870/1‑143 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTTACGCCACATGACCGCTGGGGTCGGACAACGTTTCGCTTCCTTCCAGCGAAAAGGGGAAACCGCCGGCATAACGGAAAATGTTAACGACCATTGGTAAATGCTGGgcg                                   >  2:1514659/1‑150 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAgcg                                   >  2:511299/1‑150 (MQ=255)
                                                                                                                 tAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGGAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAgcgcg                                 >  2:244778/1‑150 (MQ=255)
                                                                                                                      aaTAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACt                           <  2:466047/151‑1 (MQ=255)
                                                                                                                        tAGCGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTa                         <  1:1499884/151‑1 (MQ=255)
                                                                                                                          gCGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCACATGACCGCTGGCGTCGGGCACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGGCCGGCATAACGCAAATCGTTACCGACCATTTATAATTGCTGAGCGCGCCAACTTaa                        >  1:1214170/1‑150 (MQ=255)
                                                                                                                           cGGCTCGGTTTATCGGCGGCGGGGGGATTTAACGCCCCATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCTCGCAAACTTaa                        >  1:755516/1‑149 (MQ=255)
                                                                                                                           cGGCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAAc                       <  1:805535/150‑1 (MQ=255)
                                                                                                                             gCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTg                     <  1:5314/150‑1 (MQ=255)
                                                                                                                             gCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTg                     <  2:1144542/150‑1 (MQ=255)
                                                                                                                             gCTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTg                     <  2:1403069/150‑1 (MQ=255)
                                                                                                                              cTCGGTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGc                    <  2:1293718/150‑1 (MQ=255)
                                                                                                                                  gTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCCTTGATAATTGCTGAGCGCGCAAACTTAACTGCAcgcg               >  1:931643/1‑151 (MQ=255)
                                                                                                                                  gTTTATCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCAcgcg               >  2:382206/1‑151 (MQ=255)
                                                                                                                                      aTCGGCGGCGGGGAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCa           >  2:233845/1‑151 (MQ=255)
                                                                                                                                               ggggAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTAc   >  1:269350/1‑150 (MQ=255)
                                                                                                                                                gggAGATTTAACGCCACATGACCGCTGGCGTCGGACACCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTAcc  <  1:1128664/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 < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: