Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,965,960 A→T 7.6% intergenic (‑148/+116) fliE ← / ← atoC flagellar hook‑basal body complex protein FliE/two‑component system DNA‑binding transcriptional activator AtoC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,965,9600AT 7.6% 862.9 / 46.3 305intergenic (‑148/+116)fliE/atoCflagellar hook‑basal body complex protein FliE/two‑component system DNA‑binding transcriptional activator AtoC
Reads supporting (aligned to +/- strand):  ref base A (99/166);  new base T (11/12);  total (127/178)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.82e-01

Only 100 of 311 total aligned reads displayed.

ATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAACCCG  >  NC_002947/4965811‑4966109
                                                                                                                                                     |                                                                                                                                                     
  gACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGTCGAATGACCTACCTAAAGCATAAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCaga                                                                                                                                                     <  1:3563742/150‑1 (MQ=255)
   aCCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCagag                                                                                                                                                    <  1:8611670/150‑1 (MQ=255)
      cTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCATGCGTCGCCATATGGCCGCAGAGCgg                                                                                                                                                 >  1:6675914/1‑150 (MQ=255)
           ctctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTCTGCGGcccta                                                                                                                                            <  2:3257283/150‑3 (MQ=255)
            tctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGATAGGGGCGCCGAGCGGcccccg                                                                                                                                           >  1:1884234/1‑148 (MQ=255)
             ctcCTGCGGCCGCATTTTTTTGACCGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGc                                                                                                                                          <  1:7207139/150‑1 (MQ=255)
             ctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGc                                                                                                                                          <  2:4300555/150‑1 (MQ=255)
              tcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCg                                                                                                                                         >  2:4569820/1‑150 (MQ=255)
              tcCTGCGGCCGCATTTTTTTGACAGTCAGGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTCTGCGGccctatcg                                                                                                                                         <  1:300269/150‑6 (MQ=255)
                 tGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCATGCGGCGCGATAGGGCCGCACAGCGGCCCCCGCGAta                                                                                                                                      >  2:4013636/1‑149 (MQ=255)
                  gCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCt                                                                                                                                     <  1:1105731/150‑1 (MQ=255)
                     gCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGGCGTGCGGCGCGATTGGGCCCCAGAGCGGCCCCCCCGATCGTTg                                                                                                                                  >  1:1072639/1‑150 (MQ=255)
                     gCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTg                                                                                                                                  >  2:3284930/1‑150 (MQ=255)
                     gCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCATGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTg                                                                                                                                  >  2:3886053/1‑150 (MQ=255)
                      ccGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAACGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgc                                                                                                                                 <  2:1141323/150‑1 (MQ=255)
                      ccGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgc                                                                                                                                 <  2:1055107/150‑1 (MQ=255)
                        gCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcgg                                                                                                                               <  1:507520/150‑1 (MQ=255)
                          aTTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                                             >  1:8027610/1‑150 (MQ=255)
                          aTTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                                             <  2:8696209/150‑1 (MQ=255)
                            ttttttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAAGTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGa                                                                                                                           <  2:5954971/150‑1 (MQ=255)
                              ttttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGAACTTTGCGGCGGAGc                                                                                                                         >  2:572372/1‑150 (MQ=255)
                               tttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCt                                                                                                                        >  2:1939156/1‑150 (MQ=255)
                                ttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGGCGGGGGGGGCCCGTGGGGCGCAGAGCGGCCCCGGCGCTCGTTGGGGGGGCGCTg                                                                                                                       >  1:341786/1‑150 (MQ=255)
                                ttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTg                                                                                                                       <  1:911478/150‑1 (MQ=255)
                                   aCAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAg                                                                                                                    >  2:520908/1‑150 (MQ=255)
                                   aCAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGCTAGGGCCGCCGCGCGGCCCCGGCGCGCGTTGCGGCGGCGCGGAAg                                                                                                                    >  1:4590000/1‑150 (MQ=255)
                                      gTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGAGGCGGCGCAGAAGAcc                                                                                                                 >  2:2894821/1‑150 (MQ=255)
                                         aTGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTgg                                                                                                              <  2:8027610/150‑1 (MQ=255)
                                           gcggcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCCGCGATCATTGCGGCGGAGCTGAAGACCTgggg                                                                                                            >  2:6762882/1‑150 (MQ=255)
                                            cggcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGAt                                                                                                                                       >  2:1699834/1‑122 (MQ=255)
                                            cggcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGc                                                                                                           <  1:4177273/150‑1 (MQ=255)
                                            cggcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGc                                                                                                           <  1:2853191/150‑1 (MQ=255)
                                              gcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCCTGGGGCGCCGGAGGGCCCCGAAGCGGCCCCCGCGATGGGTTGGGCGGGGGAGAAGACCTGGGGGCg                                                                                                         >  1:266297/1‑150 (MQ=255)
                                               cgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGc                                                                                                        <  1:2591286/150‑1 (MQ=255)
                                                gAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGGCGTTCGTCCGGAATGGGCCCCCCCGCGGCCCCCGCCCGCATTGGGGCGGGGCGGGAGACCACGGGgcgcc                                                                                                       >  1:1670227/1‑149 (MQ=255)
                                                  aTGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGCTAGGGCCGCCGAGCGGCCCCCGCGATCTTTGCGGCGGTGCTGAAGACCTGGGGCCGCGTc                                                                                                     >  2:4799482/1‑150 (MQ=255)
                                                    gACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTcgc                                                                                                   <  2:1758082/150‑1 (MQ=255)
                                                        tacctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGccc                                                                                               <  1:4102374/150‑1 (MQ=255)
                                                           ctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGGCGTGGGGCGCGGTAGGGCCGCCCAGCCGCCCCCCCCCGCTTTTGGGCGGGGCGGGAGACCAGGGGGGGCGGCGCGccccca                                                                                            >  1:3175090/1‑148 (MQ=255)
                                                           ctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGGCGTGCGGCGCGCTAGGGGCGCCCGGCGGCCCCCGCGGTCGTTGGGGCGGAGGTGGAGAACAGGGGGCGCCTCGCGcccccg                                                                                            >  1:8364538/1‑148 (MQ=255)
                                                            taAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCAtt                                                                                           >  2:3574243/1‑150 (MQ=255)
                                                             aaaGCAACACAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTc                                                                                          <  1:5280230/150‑1 (MQ=255)
                                                             aaaGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTc                                                                                          <  2:797757/150‑1 (MQ=255)
                                                              aaGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCATGCGTCGCGATAGGGGCGCAGAGCGGCCCCCGCGAACGTTGCGGCGGGGCTGAAGACCAGGGGCCGCGTCGcccccccaacg                                                                                         >  1:1341467/1‑145 (MQ=255)
                                                                gCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGGCGCGATAGGGGCGCAGCGCGGCCCCCGCGGGCTTTGCGGCGGGGCGGTAGACCAGGGGCCGCGTCGCGCCCCaatagcc                                                                                       >  2:3942423/1‑144 (MQ=255)
                                                                 caacaaAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGa                                                                                      <  1:6842016/150‑1 (MQ=255)
                                                                     aaaGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATATGGCCGCAGAGCGGCCCCCGCGATCTTTGCGGCGGAGCGGAAGACCTGGGGCCGCGTCGCGCCCCCTTCGCGAcaca                                                                                  >  2:5420294/1‑150 (MQ=39)
                                                                      caGGTGCCAGCCTTCGCAACTCATGTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTGGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACaa                                                                                 <  2:7000261/149‑1 (MQ=255)
                                                                      aaGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACaa                                                                                 <  2:3364263/150‑1 (MQ=255)
                                                                      aaGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACaa                                                                                 <  1:4875183/150‑1 (MQ=255)
                                                                      aaGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACaa                                                                                 <  2:2311583/150‑1 (MQ=255)
                                                                       aGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAg                                                                                <  2:1225185/150‑1 (MQ=255)
                                                                       aGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGGCGCGCTAGGGCCCCAGAGCGGCCCCCGCGCGCTGTGCGGCGGGGCTGGAGACCCGGGGCCGCGGCGCGCCCCCATCCCGACACAAg                                                                                >  1:6029176/1‑150 (MQ=255)
                                                                        ggTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAgg                                                                               <  2:2104886/150‑1 (MQ=255)
                                                                          tGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCAGGGGCCGCTTCGCGCCCCATTCGCGACACAACGcc                                                                             >  2:1440773/1‑150 (MQ=255)
                                                                            ccAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCACGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGt                                                                           >  1:2267915/1‑150 (MQ=255)
                                                                             cAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGtt                                                                          <  2:5525484/150‑1 (MQ=255)
                                                                                 cTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCca                                                                      <  2:2268445/150‑1 (MQ=255)
                                                                                 cTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCca                                                                      <  2:5813802/150‑1 (MQ=255)
                                                                                 cTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCCGCCCCCGCGATCTTTTCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCca                                                                      >  1:5862490/1‑150 (MQ=255)
                                                                                    cGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACag                                                                   <  1:3292351/150‑1 (MQ=255)
                                                                                     gCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACaga                                                                  <  2:3130459/150‑1 (MQ=255)
                                                                                     gCAACTCATGTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACaga                                                                  <  2:8339297/150‑1 (MQ=255)
                                                                                       aaCTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTTCGGCGGAGCTGAAGACCAGGGGCCCCTTCGCGCCCCATTAGCGACACAAGGCCGTTACCACAGAgg                                                                >  2:4201861/1‑150 (MQ=39)
                                                                                           cATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCcg                                                            <  2:634139/150‑1 (MQ=255)
                                                                                            aTTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCcgc                                                           <  2:2541375/150‑1 (MQ=255)
                                                                                            aTTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCAATCGCGACACAAGGCCGTTCCCACAGAGGTCcgc                                                           >  2:8755532/1‑150 (MQ=255)
                                                                                               tCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGcc                                                        <  1:5643309/150‑1 (MQ=255)
                                                                                                 ccAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAg                                                      >  2:4990102/1‑150 (MQ=255)
                                                                                                   agagATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGtt                                                    <  2:1171667/150‑1 (MQ=255)
                                                                                                            acaATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAgatgat                                           <  1:4406565/150‑1 (MQ=255)
                                                                                                            acaATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCGGCGCGCCCCAGTCGCCACACAAGGCCGTTCCCACACAGGTCCGCCCCAGTTTCAcatgat                                           >  2:7086146/1‑150 (MQ=37)
                                                                                                              aaTTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGAttt                                         <  1:4569820/150‑1 (MQ=255)
                                                                                                              aaTTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTCGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGAttt                                         <  2:2817291/150‑1 (MQ=255)
                                                                                                                  ccTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCa                                     <  2:2061418/150‑1 (MQ=255)
                                                                                                                    tGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGt                                   <  2:1578998/150‑1 (MQ=255)
                                                                                                                    tGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCCGCCATCTTTTCGGCGGAGCTGAAGAACAGGGGCCGCTTCGCGCCCCAATAGCGACACAACGCCGTTCCCACACAGGTCCGCGCCACTTGCAGATGAATTGCCAGt                                   >  2:3506275/1‑150 (MQ=37)
                                                                                                                      cGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCa                                 <  2:4514838/150‑1 (MQ=255)
                                                                                                                      cGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCa                                 <  1:100368/150‑1 (MQ=255)
                                                                                                                        ggAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGc                               <  2:6672177/150‑1 (MQ=255)
                                                                                                                         gAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAACGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGcc                              >  1:6982434/1‑150 (MQ=255)
                                                                                                                             ggCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAAc                          <  1:4674346/150‑1 (MQ=255)
                                                                                                                              gCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACa                         <  1:6483505/150‑1 (MQ=255)
                                                                                                                              gCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACa                         >  1:5051355/1‑150 (MQ=255)
                                                                                                                              gCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACa                         <  1:788878/150‑1 (MQ=255)
                                                                                                                               ccTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAg                        >  2:7809919/1‑150 (MQ=255)
                                                                                                                                   gCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTg                    <  1:754783/150‑1 (MQ=255)
                                                                                                                                     gTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCGTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGc                  >  2:3652825/1‑150 (MQ=255)
                                                                                                                                      tCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGcc                 >  1:6174146/1‑150 (MQ=255)
                                                                                                                                        gcgATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTc               <  1:3886053/150‑1 (MQ=255)
                                                                                                                                        gcgATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTc               <  1:1939156/150‑1 (MQ=255)
                                                                                                                                          gATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGa             >  2:7354016/1‑150 (MQ=255)
                                                                                                                                           ccggggCCGCTCTGCGGCCCTGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGAc            >  2:8473720/4‑150 (MQ=39)
                                                                                                                                            cggggCCGCTCTGCGGCCCTGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACg           >  1:6403552/3‑150 (MQ=255)
                                                                                                                                             aGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGt          >  1:7361077/1‑150 (MQ=255)
                                                                                                                                              gggCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTc         <  1:1174109/150‑1 (MQ=255)
                                                                                                                                               ggCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCa        >  2:4085748/1‑150 (MQ=255)
                                                                                                                                                  cGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAAc     <  2:953641/150‑1 (MQ=255)
                                                                                                                                                  cGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAAc     <  2:6916765/150‑1 (MQ=255)
                                                                                                                                                    cAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAAccc   <  2:4626032/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAACCCG  >  NC_002947/4965811‑4966109

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: