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

GACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGC  >  NC_000913/1129821‑1130290
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            
gACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGATTTTGCCGCCCGTCCTTTCGTTACCg                                                                                                                                                                                                                                           >  1:104809/1‑238 (MQ=255)
     gtCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGCTTTTTCCGCCCCGCTGTTCGTTTACGGcgcg                                                                                                                                                                                                                                      >  1:293384/1‑238 (MQ=255)
     gtCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGAAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACTTCACACTGGTGCTGGTGGAGCCGGCTGCTTTTGCCTCCCAGCTTCTCGGTTCCcgcgcgc                                                                                                                                                                                                                                     >  1:928064/1‑238 (MQ=255)
             aTCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGGTGTTTTTGCCGCCCGGCTGTTCGTTTCCGGCGCGTCCGTgtgt                                                                                                                                                                                                                            >  1:138854/1‑236 (MQ=255)
                ggtttttttGCAAGTTTTGTTGCTGTTTGTTGATCAGCGCTTCGGCAATTTTGCCGGTGTCAATTTTTAGTTCACCGTTTCGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTcg                                                                                                                                                                                                                         <  2:107607/231‑1 (MQ=255)
                      tctGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGGGGAGGCGGTTGTTTTTGCCGCCCCGCGGTTCTTTTACGCCGCGGCACGGGGTTcgccgggt                                                                                                                                                                                                                   >  1:573654/1‑234 (MQ=255)
                          gcaGTTTTTCTGCTTTTCTTTGATTCGCGGATTGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGaa                                                                                                                                                                                                               <  2:447027/237‑1 (MQ=255)
                          cAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTTTTTTTTCCGCCCCGCTGTTCGTTTCCCTCGGGTCACGGGGTTTGCGCGGtggga                                                                                                                                                                                                               >  1:721263/1‑235 (MQ=255)
                          cAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTTCCTCCCGGCTGGTCGTTTCCCGCCCGTCAGTGGGTTCTCGCGGGTGac                                                                                                                                                                                                               >  1:39980/1‑238 (MQ=255)
                          cAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATAACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTGTTTTTTTCCGCCGCGCGTTTTGTTTCCGGCCGGCCTGTGGTTTcggcgggtggta                                                                                                                                                                                                               >  1:249289/1‑229 (MQ=255)
                               cTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACg                                                                                                                                                                                                             <  1:918227/236‑1 (MQ=255)
                                  gctgcGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGCTGTTTTTGTCCACCCGCCGTTCGTTTCCGGCGCGTCCGTGGGTTTTCGCGGtgttacggggcc                                                                                                                                                                                                        >  1:494920/1‑227 (MQ=255)
                                   ctgcGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTTTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTa                                                                                                                                                                                                      <  1:176317/239‑1 (MQ=255)
                                   ccgcgcTTCGTTGATCAGGGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCCCCGTTACGAATCGCCAGTTTTAACGCCTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTa                                                                                                                                                                                                      <  1:916592/237‑1 (MQ=255)
                                     gcgtTTCGTTTGTCCGCGGATCGGCAATTTTTGCGGGGTCCATTTTTAGTTCACCGTTGCGAATCGCCGGTTTTAACGGTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACa                                                                                                                                                                                                    <  2:158995/239‑1 (MQ=255)
                                        cTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCCGTTCGTTACCGGCGCGGCCGTGGGTTTGCGCGGGTGGACGGGGCtttaaggg                                                                                                                                                                                                 >  1:211468/1‑233 (MQ=255)
                                           cGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTa                                                                                                                                                                                                      <  1:713294/231‑1 (MQ=255)
                                                   aGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGa                                                                                                                                                                                      <  1:258315/239‑1 (MQ=255)
                                                     cgcATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGGTGTTTTTTCCCTCCCGGTTGTCTTTTCCCGCGCCCCCGTGGTTTTGCGGCGTTGTACGGCGCTTACAGGGTTCaacccgccct                                                                                                                                                                                    >  2:713294/1‑230 (MQ=255)
                                                      ggatggccAATTGTGCGTGTATCCTTTTGTAGTGCACCGTTAACAGTCTCACTTTTTAACGCTTCTACCCGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCGGGTTTGCTTGCGCGTCGCTTAACGTCATCCTGGTCCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGt                                                                                                                                                                                   <  2:701443/232‑1 (MQ=255)
                                                        aTCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTg                                                                                                                                                                                  <  1:677405/238‑1 (MQ=255)
                                                          cGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTTTTTTTTCCGCCCGGCCGTTCGTTTCCCGCGCGTCCGCGGGTTTTCCGGGGTTGACCGGGCTTTCCAGCGTCCGAGGCGGAGTggg                                                                                                                                                                                >  2:277942/1‑236 (MQ=255)
                                                           ggCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGAt                                                                                                                                                                              <  1:216785/239‑1 (MQ=255)
                                                             taaTTTTGCGGGTTTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGGTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTGTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATc                                                                                                                                                                             <  2:158165/237‑1 (MQ=255)
                                                              aaTTTTGCCGGTTTCCTTTTTTAGTTCACCGTTCCGGATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGTAGTGCGATCaa                                                                                                                                                                           <  2:890482/239‑1 (MQ=255)
                                                              aaTGTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTGCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCaa                                                                                                                                                                           <  1:341232/239‑1 (MQ=255)
                                                                    gCCGGTTTCCTTTTTTATTTCACCGTTACGAAGCGTCAGTTTTAACGCTTCGACACGTTCAAGAGTGATATCCCTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCa                                                                                                                                                                     <  2:223596/239‑1 (MQ=255)
                                                                    gCCGGTGTCCATTTTTAGTTCCCCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCa                                                                                                                                                                     <  1:713660/239‑1 (MQ=255)
                                                                     tgcaggtCCATTTTTATTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCa                                                                                                                                                                     <  1:871227/233‑1 (MQ=255)
                                                                              aTTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTc                                                                                                                                                           <  1:589782/239‑1 (MQ=255)
                                                                                         ccctgtACGAATCGCCAGTTTTAACCCTTCGAAACGTTCAAGATTGATCTCACTGCTGCCGGTTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAggg                                                                                                                                                <  1:416048/234‑1 (MQ=255)
                                                                                         cccGTTTCGAATCGGCAGTTTTAACCGTTCGACACGGTCAAGGTTGATTACTCTGCGTCCGGGTTGCATCAGGTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAggg                                                                                                                                                <  2:211468/238‑1 (MQ=255)
                                                                                         aCCGTTACGCATCGCCCTTTTTAAGCCTTCGAGCGGTTCAAGGTTGGTATCACTGCTGCCGGGTTTCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAggg                                                                                                                                                <  2:138854/239‑1 (MQ=255)
                                                                                         aCCGTTACGAATCGCTAGTTTTAACGCTTCGACACGTTCAAGATTGATACCACTGGTGCCGGGTTGCATCCGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAggg                                                                                                                                                <  1:49875/239‑1 (MQ=255)
                                                                                           cGTTCCGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGgcg                                                                                                                                              <  1:575357/239‑1 (MQ=255)
                                                                                                         attttgAACGCTTCGACACGTTCAAGATTGATATCTCTGCTGGCTGGTTGCATCAGTTTTGCTTGCGAGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGGTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATgtgt                                                                                                                                <  1:388843/237‑1 (MQ=255)
                                                                                                           ttttAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCtt                                                                                                                                           <  1:63863/226‑1 (MQ=255)
                                                                                                                                  ggttGTTATCCCTGCGGCGGGGTTGCATCAGTTTTGCTGGGGGGTCGCTTAACGTCACACTGGTGCTGGGGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACggg                                                                                                       <  2:215808/237‑1 (MQ=255)
                                                                                                                                   agtGATATCCCTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACggg                                                                                                       <  1:390216/236‑1 (MQ=255)
                                                                                                                                    ttGATATCACTGCTGCCGGGTTGCATCAGTTTGGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGGTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTa                                                                                                     <  1:1003855/239‑1 (MQ=255)
                                                                                                                                            cttgcttcTGGTTTGTATCTGTTTTGTTTGTGCGTCTCTTAACGTCACGCTGGGGCTGGTGGTGTCGTTTGTTTGTGCCGCCCGGCTGTTCGTTACCGGCGCGGTAGTGGTTCCCCGCGGTTGAACGGTGCTTGCAGGCTTCAGAGGCGAAGTGCGATCAATTCTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTaa                                                                                             <  2:128767/237‑1 (MQ=255)
                                                                                                                                            actgtgtgagggTTACAACATTTTGGGTTGCGAGTCGCTTAAGGTCCCACTGGTGCTGGTGGAGGCGGTTGTTTTTGCAGTCGGGCTGTACGTTTCCGGCGCGTCAGTGGTTTGGCGCGGTTGAACGGTGCTTACAGGATTAAGAGGCGAAGTGCGATAAATACTCATGGTTTCTTCCTCATTCAGGGCGCTTGTATCATGTGTTGCGTATTAATCGGCAAGGGACGGGTGATCTTTaa                                                                                             <  2:665263/230‑1 (MQ=255)
                                                                                                                                            aCTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTaa                                                                                             <  1:560792/239‑1 (MQ=255)
                                                                                                                                            aCTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTaa                                                                                             <  1:940272/239‑1 (MQ=255)
                                                                                                                                                    ggggTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTAc                                                                                     <  1:393112/238‑1 (MQ=255)
                                                                                                                                                    cGGGTTGCATCAGTTTTGCTTGCGGGTCGCTTAACGTCACCCTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTAc                                                                                     <  1:334351/239‑1 (MQ=255)
                                                                                                                                                    cGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTAc                                                                                     <  1:198525/239‑1 (MQ=255)
                                                                                                                                                                   ttgcttgcGCGTCGCTTAACGGCACGCTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCTGCGCGGCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAata                                                                      <  1:800695/239‑1 (MQ=255)
                                                                                                                                                                     ggctgcGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATAtt                                                                    <  1:277942/236‑1 (MQ=255)
                                                                                                                                                                         gcggGTCGGTTAACGTCACACGGGTGCTGGTGGAGGCGGTTGTTTGTGCCGCCCGGGTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTTCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCa                                                                <  1:1039838/239‑1 (MQ=255)
                                                                                                                                                                         gcggGTCGCTTAACGTTCCACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCa                                                                <  1:818839/239‑1 (MQ=255)
                                                                                                                                                                           gcgTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGGCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAt                                                               <  1:875125/238‑1 (MQ=255)
                                                                                                                                                                                         tcctttggtgcgggggGGGGGTGTTTTTTTTGCCGCACGGGTGTTTGTTTCCGGCGCGTGAGTGTTTTCGCGCGGTTGAACGGTGTTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCAGTGAGGGCGCTTTTAACATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAAcgc                                                <  2:494920/234‑1 (MQ=255)
                                                                                                                                                                                         cacaCTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGt                                                                                                                                                                                   <  1:968272/108‑1 (MQ=255)
                                                                                                                                                                                                    tggtggAGGCGGTTTTTGTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                                     <  1:458248/239‑1 (MQ=255)
                                                                                                                                                                                                    tgggggAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCCGCGCGTCAGTGGTTTCGCGCGTTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                                     <  1:269411/239‑1 (MQ=255)
                                                                                                                                                                                                    tcggggtgggtgttgttTTTGCCGCCCGGGTGTTCGTTACCTGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTATATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                                     <  2:625821/228‑1 (MQ=255)
                                                                                                                                                                                                     ggtggGGGGGTTTTTTTTGGCCGCCCGGGTGTTCGTTACCGGGGCGTCAGTGGTTTCGCGGGGTTTAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                                     <  2:400830/238‑1 (MQ=255)
                                                                                                                                                                                                     ggtggAGGCGGTTGTTTTTGCCGCCCGGCTTTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                                     <  1:948835/238‑1 (MQ=255)
                                                                                                                                                                                                     gggggAGGGGGTTTTTTTTGCCGCCCGGGTGTTTGTTACCGGCGCGGCAGTGGTTTCGCGCGGTTGAACGGTGCTTTCAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTAcc                                     <  2:204917/235‑1 (MQ=255)
                                                                                                                                                                                                                         ataagggaaaggTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATCCTCAGGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCa                  <  1:558839/227‑1 (MQ=255)
                                                                                                                                                                                                                                ggCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCa          <  1:835092/239‑1 (MQ=255)
                                                                                                                                                                                                                                        gTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTTCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTTATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGc  <  1:25496/239‑1 (MQ=255)
                                                                                                                                                                                                                                        gTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGc  <  1:952358/239‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                            
GACGTGTCATACGATCAGTTACTCTGCAAGTCTTGCTGCGCTTCGTTGATCAGCGCATCGGCAATTTTGCCGGTGTCCATTTTTAGTTCACCGTTACGAATCGCCAGTTTTAACGCTTCGACACGTTCAAGATTGATATCACTGCTGCCGGGTTGCATCAGTTTTGCTTGCGCGTCGCTTAACGTCACACTGGTGCTGGTGGAGGCGGTTGTTTTTGCCGCCCGGCTGTTCGTTACCGGCGCGTCAGTGGTTTCGCGCGGTTGAACGGTGCTTACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCATCTGCATCAACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGC  >  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: