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

CTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACC  >  NC_000913/2332676‑2332969
                                                                                                                                                   |                                                                                                                                                    
cTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑t                                                                                                                                                  >  2:1052138/1‑151 (MQ=255)
cTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑t                                                                                                                                                  >  1:699030/1‑151 (MQ=255)
  aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑ttt                                                                                                                                                >  2:1478728/1‑151 (MQ=255)
  aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑ttt                                                                                                                                                >  2:1061487/1‑151 (MQ=255)
  aTTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCCGTTTATCGGCGGCCGGGGGG‑ttt                                                                                                                                                >  2:275676/1‑151 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGG‑TTTa                                                                                                                                               >  2:1407700/1‑151 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGG‑tttt                                                                                                                                               >  1:278559/1‑150 (MQ=255)
   ttACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTAAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCCGGGGGG‑tttt                                                                                                                                               >  2:541462/1‑150 (MQ=255)
   atACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGCGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGGTAACAAAAAGCGGCTCGGTTTATCGGCGGGGGGGGGG‑TTTa                                                                                                                                               >  1:1293085/2‑151 (MQ=255)
     aCTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGG‑TTTaa                                                                                                                                              >  2:1353663/1‑150 (MQ=255)
     aCTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑ttttaa                                                                                                                                             >  2:1082144/1‑148 (MQ=255)
     aCTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑tttta                                                                                                                                              >  1:869430/1‑148 (MQ=255)
     aCTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGCCGGCGGGGTGA‑TTTaaa                                                                                                                                             >  1:1172564/1‑150 (MQ=255)
       tGTATTGTGCGGCAGTACTTAATGAGTAATGCGGCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGc                                                                                                                                           <  1:159995/151‑1 (MQ=255)
           ttGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑GTTTAAGCccca                                                                                                                                       >  1:1431566/1‑151 (MQ=255)
                cGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGAc                                                                                                                                   <  1:1073418/150‑1 (MQ=255)
                      tACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGG‑TTTTACGCCACATGACCGCCgg                                                                                                                             >  2:1332561/1‑150 (MQ=255)
                      tACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCCCATGAcccccgg                                                                                                                             >  2:1047809/1‑145 (MQ=255)
                      tACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTg                                                                                                                              <  1:630239/149‑1 (MQ=255)
                       aCTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGG‑TTTTAAGCCCCATTAcccccgggg                                                                                                                           >  2:1363310/1‑144 (MQ=255)
                       aCTTAAGGAGTAATGCGCCAGACCGCGTGCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCg                                                                                                                           <  1:1230483/151‑1 (MQ=255)
                           aaTGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTAACGCCACATGACCGCTGGCGTCgg                                                                                                                       >  2:1288673/1‑151 (MQ=255)
                             tGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGac                                                                                                                     <  2:374163/151‑1 (MQ=255)
                             tGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGac                                                                                                                     <  1:236991/151‑1 (MQ=255)
                              gAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGG‑TTTTAAGCCACATGACCCCCGGCGTCGGgca                                                                                                                    >  2:1059307/1‑151 (MQ=255)
                                 tAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGaca‑c                                                                                                                  <  1:934579/149‑1 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGG‑TTTTACGCCCCATAAACCCTGGCGCCCGGCA‑cc                                                                                                                 >  1:240165/1‑146 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGG‑TTTTACGCCACATGACCGCTGGCGGCGGGCA‑ACGtttt                                                                                                            >  2:1246081/1‑150 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTAACGCCCCATGACCGCTGGCGGCGGGCA‑CCGgttt                                                                                                            >  1:1054124/1‑150 (MQ=255)
                                     gcgcCAGACCGCGTTCAATGAGGATCTCTTTGGCGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCCCGGTTTATCGGCGGGGGGGACG‑TTTAAACCCACATGAACCCCGGCGTccgggacccgttt                                                                                                             >  2:1421812/1‑139 (MQ=255)
                                       gcCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTAACGCCCCATGACCGCCGGCGTCGGGCA‑CCCttttg                                                                                                           >  2:321035/1‑148 (MQ=255)
                                            aCCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGGTAACAAATAGCGGCTCGGTTTATCGGCGGGGGGGGGAGTTTAAC‑ACACAAGACAGCCGGCGTCGGACA‑CCCGTTCGCTGcc                                                                                                      >  1:1343252/1‑150 (MQ=255)
                                                 gTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAAGAACGGCTCGGTTTATGGGCGGCGGGGGGG‑ATTAAAGCCACATGACCGCTGGCGTCGGAAA‑CCCGTTCGCTGCCTTCCAg                                                                                                >  1:841367/1‑151 (MQ=255)
                                                      aTGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTAACGCCACATGACCCCCGGCGTCGGACA‑CCGTTTTGCTGCCCTCCCGCGac                                                                                            >  2:1479343/1‑150 (MQ=255)
                                                           gATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACgg                                                                                        <  1:13193/149‑1 (MQ=255)
                                                           gATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGGTTATCGGCGGCGGGGGGA‑TTTAAAGCCACATTACCGCTGGGGGCGGGCC‑CCGGTTCGCTGCCTTTCAGCGACACgggg                                                                                      >  1:1144001/1‑151 (MQ=255)
                                                             tctcTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTAACGCCCCCTTACCCCCGGCGTCGGGCA‑CCCTTTCGCTGCCTTCCCGCGACACggggga                                                                                    >  2:906766/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGG‑TTTAACCCCCCATAGCCCCCGGCGTCGGACC‑CCGCTTTGCGTCCTTTCAGCGACACGGGGgaaag                                                                                 >  1:1008885/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGG‑TTTAAAGGCACATGACCGCCGGGGTCGGACA‑CCGTTTCGCTGCCTTTCACCCACACGGGGaaaag                                                                                 >  2:1236811/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTTACGCCACATGACCCCTCGCGTCGGACA‑CCGTTTTGCTGCCCTCCCGCGACACGGGGaaaag                                                                                 >  2:846823/1‑149 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTAACGCCCCATGACCGCCGGCGTCGGAAA‑CCGTTTTGCTGCCTTCCAGCGACACGGGGAAACg                                                                                 >  2:86085/1‑151 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCCGCCGGGCGA‑TTTTATGCCCCATAACCGCCGGCGCCGCGAA‑CCGTGTTGCTGCCTTCCACCCGAACCGGGAAACg                                                                                 >  1:286336/1‑151 (MQ=255)
                                                                ctTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAAAAGCGGCCCGGTTTTTCGGCGGCGGGGGGA‑TTTAACGCCCACTGACCCCTGCCGCCGGAAACCCGTTTCG‑TGCCTTCCCACAACCCGGGGaaaag                                                                                 >  2:647922/1‑149 (MQ=255)
                                                                         tggtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCAt                                                                         <  2:20058/150‑1 (MQ=255)
                                                                          ggtgACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑ATTTACGCCCCAAGACCGCTGGGGTGCGACA‑CCGCTTCGCCGCCCTTCAGCGACACCGGGAAACGTCCGGCCTaa                                                                       >  2:1180647/1‑151 (MQ=255)
                                                                               cgcgATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTAAAGCCCCATGACCGCTGGCGTCGGACA‑CCCTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCaa                                                                  >  2:1242838/1‑151 (MQ=255)
                                                                               cgcgATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑GTTTACACCACATTAACCCTGGCGGCGGGCC‑CCGGTTTGGTGGCTTTCCGGGAAACGGGGGAAAGTTCgggccaacgcca                                                                  >  1:129233/1‑141 (MQ=255)
                                                                                gcgATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGGTTCGCTGCCTTCCAGCGACACGGGGGAACGGCCGGCATAACGCaaa                                                                 >  2:367263/1‑151 (MQ=255)
                                                                                gcgATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCCCATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGGAACGGCCCGCATAACGCaaa                                                                 >  2:913727/1‑151 (MQ=255)
                                                                                 cgATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGG‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTTGCTGCCTTCCAGCGGCACCGGGGAACGTCCGGCATAACGCAAAt                                                                >  1:63246/1‑151 (MQ=255)
                                                                                   aTACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCg                                                              <  2:914326/151‑1 (MQ=255)
                                                                                           cgccgTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGAc                                                      <  2:869901/151‑1 (MQ=255)
                                                                                             aggTCACTGGCGAAGACGATTAACAAAGAGCGGCTCAGGTTATCGGCGGTGGGAAGA‑TTTAACGCCAAATAACCGCTGGCGTCTGACA‑CCGTTTCGCTGCCTTCCAACGACACGGGGAAACGTCCGGCAAAACGCAAATCGTTACCGGcc                                                     <  1:750447/148‑1 (MQ=255)
                                                                                                 cACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAACGACACGGGGAAAAGGCCGGCATAACGCAAATCGTTACCGACCATTGa                                                >  2:157062/1‑151 (MQ=255)
                                                                                                   cTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATa                                              <  2:843488/151‑1 (MQ=255)
                                                                                                         gagaCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGc                                         <  2:63279/150‑1 (MQ=255)
                                                                                                            aCGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACCCGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCCTTGATAATTGCTGAg                                     >  2:1465933/1‑151 (MQ=255)
                                                                                                             aggTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTAACGCCACATAACCGCTGGCGTCGGGAA‑CCGTTTAGCTGCCTTCCAGGGAAAAGGGGAAACGTCCGGCATAAAGCCAATCGCTTACGAACATTGGTAATTGGTGAg                                     >  2:480853/2‑150 (MQ=255)
                                                                                                              ggTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAAAGCCACATTACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGGCCGGCATAACGCAAATCGTTTCCGACCCTTGATAATTGCTGAgcg                                   >  1:624545/1‑151 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGTG‑TTTTAAGCCACAAAACCGCTGGGGCGGGGCT‑CCGGTTTGCTGGCTTCCAGCCACACGGGGGAACGTTCGGGCTAACGCAAATCGTTTCCCCCCACTGAAAATTgtctagcgg                                  >  1:301870/1‑143 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTTACGCCACATGACCGCTGGGGTCGGACA‑ACGTTTCGCTTCCTTCCAGCGAAAAGGGGAAACCGCCGGCATAACGGAAAATGTTAACGACCATTGGTAAATGCTGGgcg                                   >  2:1514659/1‑150 (MQ=255)
                                                                                                               gTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAgcg                                   >  2:511299/1‑150 (MQ=255)
                                                                                                                 tAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGGAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAgcgcg                                 >  2:244778/1‑150 (MQ=255)
                                                                                                                      aaTAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACt                           <  2:466047/151‑1 (MQ=255)
                                                                                                                        tAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTa                         <  1:1499884/151‑1 (MQ=255)
                                                                                                                          gCGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTAACGCCACATGACCGCTGGCGTCGGGCA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGGCCGGCATAACGCAAATCGTTACCGACCATTTATAATTGCTGAGCGCGCCAACTTaa                        >  1:1214170/1‑150 (MQ=255)
                                                                                                                           cGGCTCGGTTTATCGGCGGCGGGGGGA‑TTTAACGCCCCATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCTCGCAAACTTaa                        >  1:755516/1‑149 (MQ=255)
                                                                                                                           cGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAAc                       <  1:805535/150‑1 (MQ=255)
                                                                                                                             gCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTg                     <  1:5314/150‑1 (MQ=255)
                                                                                                                             gCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTg                     <  2:1144542/150‑1 (MQ=255)
                                                                                                                             gCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTg                     <  2:1403069/150‑1 (MQ=255)
                                                                                                                              cTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGc                    <  2:1293718/150‑1 (MQ=255)
                                                                                                                                  gTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCCTTGATAATTGCTGAGCGCGCAAACTTAACTGCAcgcg               >  1:931643/1‑151 (MQ=255)
                                                                                                                                  gTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCAcgcg               >  2:382206/1‑151 (MQ=255)
                                                                                                                                      aTCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCa           >  2:233845/1‑151 (MQ=255)
                                                                                                                                               ggggAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTAc   >  1:269350/1‑150 (MQ=255)
                                                                                                                                                gggAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTAcc  <  1:1128664/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTATTACTGTATTGTGCGGCAGTACTTAATGAGTAATGCGCCAGACCGCGTTCAATGAGGATCTCTTTGGTGGTGGTGACGCGATACGCCGCGCCGTCACTGGCGGAGACGGTTAACAAATAGCGGCTCGGTTTATCGGCGGCGGGGAGA‑TTTAACGCCACATGACCGCTGGCGTCGGACA‑CCGTTTCGCTGCCTTCCAGCGACACGGGGAAACGTCCGGCATAACGCAAATCGTTACCGACCATTGATAATTGCTGAGCGCGCAAACTTAACTGCACGCGGGCATTTTTTACC  >  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: ‑