Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,130,055 A→T 16.8% V25E (GTA→GAA)  flgM ← anti‑sigma factor for FliA (sigma 28)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,130,0550AT16.8% 149.2 / 3.0 64V25E (GTA→GAA) flgManti‑sigma factor for FliA (sigma 28)
Reads supporting (aligned to +/- strand):  ref base A (2/50);  new base T (10/2);  total (12/52)
Fisher's exact test for biased strand distribution p-value = 2.68e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGC  >  NC_000913/1129821‑1130290
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            
gACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑ATTTTGCCGCCCGTCCTTTCGTTACCg                                                                                                                                                                                                                                           >  1:104809/1‑238 (MQ=255)
     gtCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑CTTTTTCCGCCCCGCTGTTCGTTTACGGcgcg                                                                                                                                                                                                                                      >  1:293384/1‑238 (MQ=255)
     gtCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGAAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACTTCACACTGGTGCTGGTGGAGCCGGCTG‑CTTTTGCCTCCCAGCTTCTCGGTTCCcgcgcgc                                                                                                                                                                                                                                     >  1:928064/1‑238 (MQ=255)
             aTCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGGTG‑TTTTTGCCGCCCGGCTGTTCGTTTCCGGCGCGTCCG‑Tgtgt                                                                                                                                                                                                                            >  1:138854/1‑236 (MQ=255)
                ggtttttttGCAAGTTTTGTTGCTGTTTGTTGATCAGCGCTTCGGCAATTTTGCCGGTGTCAATTTTTAGTTCACCGTTTCGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTcg                                                                                                                                                                                                                         <  2:107607/231‑1 (MQ=255)
                      tctGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGGGGAGGCGGTTG‑TTTTTGCCGCCCCGCGGTTCTTTTACGCCGCGGCAC‑GGGGTTcgccgggt                                                                                                                                                                                                                   >  1:573654/1‑234 (MQ=255)
                          gcaGTTTTTCTGCTTTTCTTTGATTCGCGGATTGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGaa                                                                                                                                                                                                               <  2:447027/237‑1 (MQ=255)
                          cAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTT‑TTTTTTCCGCCCCGCTGTTCGTTTCCCTCGGGTCAC‑GGGGTTTGCGCGGtggga                                                                                                                                                                                                               >  1:721263/1‑235 (MQ=255)
                          cAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTTCCTCCCGGCTGGTCGTTTCCCGCCCGTCAG‑TGGGTTCTCGCGGGTGac                                                                                                                                                                                                               >  1:39980/1‑238 (MQ=255)
                          cAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATAACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTGT‑TTTTTTCCGCCGCGCGTTTTGTTTCCGGCCGGCCTG‑TGGTTTcggcgggtggta                                                                                                                                                                                                               >  1:249289/1‑229 (MQ=255)
                               cTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACg                                                                                                                                                                                                             <  1:918227/236‑1 (MQ=255)
                                  gctgcGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGCTG‑TTTTTGTCCACCCGCCGTTCGTTTCCGGCGCGTCCG‑TGGGTTTTCGCGGtgttacggggcc                                                                                                                                                                                                        >  1:494920/1‑227 (MQ=255)
                                   ctgcGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTTTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTa                                                                                                                                                                                                      <  1:176317/239‑1 (MQ=255)
                                   ccgcgcTTCGTTGATCAGGGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCCCCGTTACGAATCGCCAGTTTTAACGCCTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTa                                                                                                                                                                                                      <  1:916592/237‑1 (MQ=255)
                                     gcgtTTCGTTTGTCCGCGGATCGGCAATTTTTGCGGGGTCCATTTTTAGTTCACCGTTGCGAATCGCCGGTTTTAACGGTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACa                                                                                                                                                                                                    <  2:158995/239‑1 (MQ=255)
                                        cTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCCGTTCGTTACCGGCGCGGCCG‑TGGGTTTGCGCGGGTGGACGGGGCtttaaggg                                                                                                                                                                                                 >  1:211468/1‑233 (MQ=255)
                                           cGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTa                                                                                                                                                                                                      <  1:713294/231‑1 (MQ=255)
                                                   aGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGa                                                                                                                                                                                      <  1:258315/239‑1 (MQ=255)
                                                     cgcATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGGTGTTTTTTCCCTCCCGGTTG‑TCTTTTCCCGCGCCCCCG‑TGGTTTTGCGGCGTTGTACGGCGCTTACAGGGTTCaacccgccct                                                                                                                                                                                    >  2:713294/1‑230 (MQ=255)
                                                      ggatggccAATTGTGCGTGTATCCTTTTGTAGTGCACCGTTAACAGTCTCACTTTTTAACGCTTCTACCCGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCGGGTTTGCTTGCGCGTCGCTTAACGTCATCCTGGTCCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGt                                                                                                                                                                                   <  2:701443/232‑1 (MQ=255)
                                                        aTCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTg                                                                                                                                                                                  <  1:677405/238‑1 (MQ=255)
                                                          cGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTT‑TTTTTTCCGCCCGGCCGTTCGTTTCCCGCGCGTCCGCGGGTTTTCCGGGGTTGACCGG‑GCTTTCCAGCGTCCGAGGCGGAGTggg                                                                                                                                                                                >  2:277942/1‑236 (MQ=255)
                                                           ggCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGAt                                                                                                                                                                              <  1:216785/239‑1 (MQ=255)
                                                             taaTTTTGCGGGTTTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGGTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTGTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATc                                                                                                                                                                             <  2:158165/237‑1 (MQ=255)
                                                              aaTTTTGCCGGTTTCCTTTTTTAGTTCACCGTTCCGGATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGTAGTGCGATCaa                                                                                                                                                                           <  2:890482/239‑1 (MQ=255)
                                                              aaTGTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTGCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCaa                                                                                                                                                                           <  1:341232/239‑1 (MQ=255)
                                                                    gCCGGTTTCCTTTTTTATTTCACCGTTACGAAGCGTCAGTTTTAACGCTTCGACACGTTCAAGAGTGATATCCCTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCa                                                                                                                                                                     <  2:223596/239‑1 (MQ=255)
                                                                    gCCGGTGTCCATTTTTAGTTCCCCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCa                                                                                                                                                                     <  1:713660/239‑1 (MQ=255)
                                                                     tgcaggtCCATTTTTATTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCa                                                                                                                                                                     <  1:871227/233‑1 (MQ=255)
                                                                              aTTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTc                                                                                                                                                           <  1:589782/239‑1 (MQ=255)
                                                                                         ccctgtACGAATCGCCAGTTTTAACCCTTCGAAACGTTCAAGATTGATCTCACTGCTGCCGGTTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAggg                                                                                                                                                <  1:416048/234‑1 (MQ=255)
                                                                                         cccGTTTCGAATCGGCAGTTTTAACCGTTCGACACGGTCAAGGTTGATTACTCTGCGTCCGGGTTGCATCAGGTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAggg                                                                                                                                                <  2:211468/238‑1 (MQ=255)
                                                                                         aCCGTTACGCATCGCCCTTTTTAAGCCTTCGAGCGGTTCAAGGTTGGTATCACTGCTGCCGGGTTTCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAggg                                                                                                                                                <  2:138854/239‑1 (MQ=255)
                                                                                         aCCGTTACGAATCGCTAGTTTTAACGCTTCGACACGTTCAAGATTGATACCACTGGTGCCGGGTTGCATCCGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAggg                                                                                                                                                <  1:49875/239‑1 (MQ=255)
                                                                                           cGTTCCGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGgcg                                                                                                                                              <  1:575357/239‑1 (MQ=255)
                                                                                                         attttgAACGCTTCGACACGTTCAAGATTGATATCTCTGCTGGCTGGTTGCATCAGTTTTGCTTGCGAGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGGTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATgtgt                                                                                                                                <  1:388843/237‑1 (MQ=255)
                                                                                                           ttttAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCtt                                                                                                                                           <  1:63863/226‑1 (MQ=255)
                                                                                                                                  ggttGTTATCCCTGCGGCGGGGTTGCATCAGTTTTGCTGGGGGGTCGCTTAACGTCACACTGGTGCTGGGGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACggg                                                                                                       <  2:215808/237‑1 (MQ=255)
                                                                                                                                   agtGATATCCCTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACggg                                                                                                       <  1:390216/236‑1 (MQ=255)
                                                                                                                                    ttGATATCACTGCTGCCGGGTTGCATCAGTTTGGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑GTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTa                                                                                                     <  1:1003855/239‑1 (MQ=255)
                                                                                                                                            cttgcttcTGGTTTGTATCTGTTTTGTTTGTGCGTCTCTTAACGTCACGCTGGGGCTGGTGGTGTCGTTTG‑TTTGTGCCGCCCGGCTGTTCGTTACCGGCGCGGTAG‑TGGTTCCCCGCGGTTGAACGGTGCTTGCAGGCTTCAGAGGCGAAGTGCGATCAATTCTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTaa                                                                                             <  2:128767/237‑1 (MQ=255)
                                                                                                                                            actgtgtgagggTTACAACATTTTGGGTTGCGAGTCGCTTAAGGTCCCACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCAGTCGGGCTGTACGTTTCCGGCGCGTCAG‑TGGTTTGGCGCGGTTGAACGGTGCTTACAGGATTAAGAGGCGAAGTGCGATAAATACTCATGGTTTCTTCCTCATTCAGGGCGCTTGTATCATGTGTTGCGTATTAATCGGCAAGGGACGGGTGATCTTTaa                                                                                             <  2:665263/230‑1 (MQ=255)
                                                                                                                                            aCTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTaa                                                                                             <  1:560792/239‑1 (MQ=255)
                                                                                                                                            aCTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTaa                                                                                             <  1:940272/239‑1 (MQ=255)
                                                                                                                                                    ggggTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTAc                                                                                     <  1:393112/238‑1 (MQ=255)
                                                                                                                                                    cGGGTTGCATCAGTTTTGCTTGCGGGTCGCTTAACGTCACCCTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTAc                                                                                     <  1:334351/239‑1 (MQ=255)
                                                                                                                                                    cGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTAc                                                                                     <  1:198525/239‑1 (MQ=255)
                                                                                                                                                                   ttgcttgcGCGTCGCTTAACGGCACGCTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCTGCGCGGCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAata                                                                      <  1:800695/239‑1 (MQ=255)
                                                                                                                                                                     ggctgcGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATAtt                                                                    <  1:277942/236‑1 (MQ=255)
                                                                                                                                                                         gcggGTCGGTTAACGTCACACGGGTGCTGGTGGAGGCGGTTG‑TTTGTGCCGCCCGGGTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTTCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCa                                                                <  1:1039838/239‑1 (MQ=255)
                                                                                                                                                                         gcggGTCGCTTAACGTTCCACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCa                                                                <  1:818839/239‑1 (MQ=255)
                                                                                                                                                                           gcgTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGGCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAt                                                               <  1:875125/238‑1 (MQ=255)
                                                                                                                                                                                         tcctttggtgcgggggGGGGGTGTTT‑TTTTTGCCGCACGGGTGTTTGTTTCCGGCGCGTGAG‑TGTTTTCGCGCGGTTGAACGGTGTTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCAGTGAGGGCGCTTTTAACATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAAcgc                                                <  2:494920/234‑1 (MQ=255)
                                                                                                                                                                                         cacaCTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGt                                                                                                                                                                                   <  1:968272/108‑1 (MQ=255)
                                                                                                                                                                                                    tggtggAGGCGGTTT‑TTGTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                                     <  1:458248/239‑1 (MQ=255)
                                                                                                                                                                                                    tgggggAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCCGCGCGTCAG‑TGGTTTCGCGCGTTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                                     <  1:269411/239‑1 (MQ=255)
                                                                                                                                                                                                    tcggggtgggtgttg‑ttTTTGCCGCCCGGGTGTTCGTTACCTGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTATATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                                     <  2:625821/228‑1 (MQ=255)
                                                                                                                                                                                                     ggtggGGGGGTTTT‑TTTTGGCCGCCCGGGTGTTCGTTACCGGGGCGTCAG‑TGGTTTCGCGGGGTTTAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                                     <  2:400830/238‑1 (MQ=255)
                                                                                                                                                                                                     ggtggAGGCGGTTG‑TTTTTGCCGCCCGGCTTTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                                     <  1:948835/238‑1 (MQ=255)
                                                                                                                                                                                                     gggggAGGGGGTTT‑TTTTTGCCGCCCGGGTGTTTGTTACCGGCGCGGCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTTCAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                                     <  2:204917/235‑1 (MQ=255)
                                                                                                                                                                                                                         ataagggaaaggTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATCCTCAGGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCa                  <  1:558839/227‑1 (MQ=255)
                                                                                                                                                                                                                                ggCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCa          <  1:835092/239‑1 (MQ=255)
                                                                                                                                                                                                                                        gTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTTCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTTATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGc  <  1:25496/239‑1 (MQ=255)
                                                                                                                                                                                                                                        gTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGc  <  1:952358/239‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            
GACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTG‑TTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAG‑TGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGC  >  NC_000913/1129821‑1130290

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