Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,965,962 A→T 7.7% intergenic (‑150/+114) 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,9620AT 7.7% 890.8 / 43.9 305intergenic (‑150/+114)fliE/atoCflagellar hook‑basal body complex protein FliE/two‑component system DNA‑binding transcriptional activator AtoC
Reads supporting (aligned to +/- strand):  ref base A (107/167);  new base T (12/11);  total (127/178)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

Only 100 of 315 total aligned reads displayed.

GACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAACCCGGC  >  NC_002947/4965813‑4966111
                                                                                                                                                     |                                                                                                                                                     
gACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGTCGAATGACCTACCTAAAGCATAAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCaga                                                                                                                                                       <  1:3563742/150‑1 (MQ=255)
  cccTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGGAGCGATAGGGCCGCAGAGc                                                                                                                                                     >  2:3961369/1‑150 (MQ=255)
     ttCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGc                                                                                                                                                  >  2:3148432/1‑150 (MQ=255)
          tctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGGCGTGCGGGGCGATATGGCCCCAGAGGGGcccccg                                                                                                                                             >  1:7364182/1‑148 (MQ=255)
          tctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGATAGGGGCGCCGAGCGGcccccg                                                                                                                                             >  1:1884234/1‑148 (MQ=255)
           ctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATGGGGCCGCTCTGCGGccctatc                                                                                                                                            <  2:4246101/150‑5 (MQ=255)
           ctcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGc                                                                                                                                            <  1:2049814/150‑1 (MQ=255)
            tcCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCTCTGCGGccctatcg                                                                                                                                           <  1:300204/150‑6 (MQ=255)
              cTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGAt                                                                                                                                         <  1:5063909/150‑1 (MQ=255)
                gCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCt                                                                                                                                       <  1:1105731/150‑1 (MQ=255)
                 cGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCtt                                                                                                                                      <  2:7869252/150‑1 (MQ=255)
                   gCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGGCGTGCGGCGCGATTGGGCCCCAGAGCGGCCCCCCCGATCGTTg                                                                                                                                    >  1:1072639/1‑150 (MQ=255)
                    ccGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCCGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgc                                                                                                                                   >  2:1066388/1‑150 (MQ=255)
                      gCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcgg                                                                                                                                 <  1:452454/150‑1 (MQ=255)
                        aTTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                                               <  2:8696209/150‑1 (MQ=255)
                        aTTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTgcggcg                                                                                                                               <  1:6438655/150‑1 (MQ=255)
                          ttttttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCGTTTCGGCGGa                                                                                                                             >  2:4559280/1‑150 (MQ=255)
                             tttGACAGTCAGGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCt                                                                                                                          <  1:797525/150‑1 (MQ=255)
                              ttGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGGCGGGGGGGGCCCGTGGGGCGCAGAGCGGCCCCGGCGCTCGTTGGGGGGGCGCTg                                                                                                                         >  1:341786/1‑150 (MQ=255)
                                gACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGaa                                                                                                                       <  2:2446535/150‑1 (MQ=255)
                                   aGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGCGCGGCCCCCGCGATCGTTGCGGCGGAGCGGAAGAc                                                                                                                    >  2:4881504/1‑150 (MQ=255)
                                       aTGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTgg                                                                                                                <  2:3601455/150‑1 (MQ=255)
                                       aTGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTgg                                                                                                                <  2:8027610/150‑1 (MQ=255)
                                        tGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTggg                                                                                                               <  2:2304727/150‑1 (MQ=255)
                                        tGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTggg                                                                                                               <  1:7763144/150‑1 (MQ=255)
                                        tGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGCCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTggg                                                                                                               <  2:5461128/150‑1 (MQ=255)
                                         gcggcgAATGACCTACCTACAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAAGTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTgggg                                                                                                              <  2:1786866/150‑1 (MQ=255)
                                         gcggcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCCGCGATCATTGCGGCGGAGCTGAAGACCTgggg                                                                                                              >  2:6762882/1‑150 (MQ=255)
                                           ggcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGcc                                                                                                            >  2:5723883/1‑150 (MQ=255)
                                            gcgAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCCTGGGGCGCCGGAGGGCCCCGAAGCGGCCCCCGCGATGGGTTGGGCGGGGGAGAAGACCTGGGGGCg                                                                                                           >  1:266297/1‑150 (MQ=255)
                                              gAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGGCGTTCGTCCGGAATGGGCCCCCCCGCGGCCCCCGCCCGCATTGGGGCGGGGCGGGAGACCACGGGgcgcc                                                                                                         >  1:1670227/1‑149 (MQ=255)
                                                  gACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTcgc                                                                                                     <  2:6181573/150‑1 (MQ=255)
                                                  gACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTcgc                                                                                                     <  1:6738598/150‑1 (MQ=255)
                                                    cctacctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCgcgc                                                                                                   <  1:4164543/150‑1 (MQ=255)
                                                     ctacctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGcc                                                                                                  <  1:2834652/150‑1 (MQ=255)
                                                     ctacctaAAGCAACAAAGGTGCCAGCCTTCGCAACTCATGTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGTGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGcc                                                                                                  <  1:6975679/150‑1 (MQ=255)
                                                            aaGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCATGCGTCGCGATAGGGGCGCAGAGCGGCCCCCGCGAACGTTGCGGCGGGGCTGAAGACCAGGGGCCGCGTCGcccccccaacg                                                                                           >  1:1341467/1‑145 (MQ=255)
                                                                    caGGTGCCAGCCTTCGCAACTCATGTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTGGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACaa                                                                                   <  2:7000261/149‑1 (MQ=255)
                                                                    aaGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACaa                                                                                   <  2:3364263/150‑1 (MQ=255)
                                                                     aGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAg                                                                                  <  2:1224969/150‑1 (MQ=255)
                                                                      ggTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAgg                                                                                 <  1:520908/150‑1 (MQ=255)
                                                                        ggCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGcc                                                                               <  2:5935790/149‑1 (MQ=255)
                                                                          ccAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGGCGCGATAGGGCCGCCGAGCGGCCCCCGCGATCTTTTCGGCGGCGCGGGAGAACGGGGGCCGCGGCGCGCCCCCATCGCGACACAACGCCGt                                                                             >  2:1593100/1‑150 (MQ=255)
                                                                           cAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGATAGGGCCGCACAGCGGCCCCGGCGAACGTTGCGGCGGAGCGGGAGACCAGGGGCCGCGTCGCGCCCCAATCGCGACACAAGGCCggt                                                                            >  2:6613477/1‑148 (MQ=255)
                                                                            aGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTc                                                                           <  2:199207/150‑1 (MQ=255)
                                                                             gCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTcc                                                                          <  1:6639261/150‑1 (MQ=255)
                                                                             gCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTcc                                                                          <  2:3315857/150‑1 (MQ=255)
                                                                                 tCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCcaca                                                                      <  2:1709879/150‑1 (MQ=255)
                                                                                  cGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACag                                                                     <  1:3292351/150‑1 (MQ=255)
                                                                                   gCAACTCATGTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACaga                                                                    <  2:8339297/150‑1 (MQ=255)
                                                                                     aaCTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGATATGGCCGCCGAGCGGCCCCCGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCGTCGCGCCCCCATCGCGACACAAGGCCCTTCCCACAGAgg                                                                  >  2:7699751/1‑150 (MQ=255)
                                                                                     aaCTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCATTCGTCGCGATAGGGCCGCAGAGCGGCCCCCGCGAACTTTTCGGCGGAGCTGAAGAACAGGGGCCGCTTCGCGCCCCATTAGCGACACAAGGCCGTTACCACAGAgg                                                                  >  2:3283074/1‑150 (MQ=255)
                                                                                       cTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTc                                                                <  2:3017971/150‑1 (MQ=255)
                                                                                          aTTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCcgc                                                             <  2:7162191/150‑1 (MQ=255)
                                                                                          aTTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCcgc                                                             <  2:2839047/150‑1 (MQ=255)
                                                                                          aTTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCAATCGCGACACAAGGCCGTTCCCACAGAGGTCcgc                                                             >  2:8755532/1‑150 (MQ=255)
                                                                                            ttCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCgcgc                                                           <  1:4300984/150‑1 (MQ=255)
                                                                                             tCCCAGAGATCACACAATTCCTGCGGGAGCGGCCATGCGTCGCGATATGGCCGCAGAGCGGCCCCCGCGCTCTTTTCGGCGGAGCTGAAAACCAGGGGCCGCTTAGCGCCCCAATCGCGACACAAGGCCGTTACCACACAGGTCCGCGcc                                                          >  1:5864896/1‑150 (MQ=255)
                                                                                              cccAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCa                                                         <  1:4982191/150‑1 (MQ=255)
                                                                                                 agagATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGtt                                                      <  1:460056/150‑1 (MQ=255)
                                                                                                 agagATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGtt                                                      <  2:1171657/150‑1 (MQ=255)
                                                                                                   agaTCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGTGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGc                                                    <  2:4031225/150‑1 (MQ=255)
                                                                                                    gaTCACACACTTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCa                                                   <  2:4680327/150‑1 (MQ=255)
                                                                                                    gaTCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCa                                                   <  1:8581669/150‑1 (MQ=255)
                                                                                                    gaTCACACAATTCCTGCGGGAGCGGCCGTGCGTCGCGCTAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCGGGGGCCGCGTCGCGCCCCAATCGCGCCCCAAGGCCGTTCCCACAGAGGTCCGCGCCCGTTGCa                                                   >  2:4180554/1‑150 (MQ=38)
                                                                                                      tCACACAATTCCTGCGGGAGCGGCCGTGCGGCGCGCTAGGGCCGCCGAGCGGCCCCCGCGCGCGTTGCGGCGGCGCGGAAGACCCGGGGCCGCGTCGCGCCCCCCGCGCCACCCCACGCCCGTCCCACACAGGTCCGCGCCCGGTGccgc                                                 >  2:5161297/1‑147 (MQ=255)
                                                                                                          acaATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAgatgat                                             <  2:2618912/150‑1 (MQ=255)
                                                                                                          acaATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCGGCGCGCCCCAGTCGCCACACAAGGCCGTTCCCACACAGGTCCGCCCCAGTTTCAcatgat                                             >  2:7086146/1‑150 (MQ=37)
                                                                                                            aaTTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGAttt                                           <  1:3148432/150‑1 (MQ=255)
                                                                                                              ttCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGc                                         <  1:7321213/150‑1 (MQ=255)
                                                                                                                  tGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGt                                     <  2:3678571/150‑1 (MQ=255)
                                                                                                                  tGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCCGCCATCTTTTCGGCGGAGCTGAAGAACAGGGGCCGCTTCGCGCCCCAATAGCGACACAACGCCGTTCCCACACAGGTCCGCGCCACTTGCAGATGAATTGCCAGt                                     >  2:3506275/1‑150 (MQ=37)
                                                                                                                    cGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCa                                   <  1:100368/150‑1 (MQ=255)
                                                                                                                    cGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCa                                   <  1:7020554/150‑1 (MQ=255)
                                                                                                                     gggAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAg                                  <  1:6134647/150‑1 (MQ=255)
                                                                                                                       gAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGcc                                >  2:5408811/1‑150 (MQ=255)
                                                                                                                       gAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGcc                                >  1:550927/1‑150 (MQ=255)
                                                                                                                        aGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCg                               <  1:904732/150‑1 (MQ=255)
                                                                                                                           ggCCTTGCGTCGCGATAGGGGCGCTGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAAc                            <  1:4674561/150‑1 (MQ=255)
                                                                                                                              cTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAgg                         >  2:5599710/1‑150 (MQ=255)
                                                                                                                              cTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAgg                         >  2:4378233/1‑150 (MQ=255)
                                                                                                                              cTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAgg                         >  2:2099494/1‑150 (MQ=255)
                                                                                                                                 gCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTg                      >  1:5695795/1‑150 (MQ=255)
                                                                                                                                 gCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTg                      <  1:2358315/150‑1 (MQ=255)
                                                                                                                                  cGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTgg                     <  2:2262017/150‑1 (MQ=255)
                                                                                                                                      gcgATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTc                 <  1:3886053/150‑1 (MQ=255)
                                                                                                                                        gATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGa               >  2:3891370/1‑150 (MQ=255)
                                                                                                                                         ccggggCCGCTCTGCGGCCCTGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACCGGCTGGCCTCGAc              >  2:1467693/4‑150 (MQ=38)
                                                                                                                                         ccggggCCGCTCTGCGGCCCTGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGAc              >  2:8473720/4‑150 (MQ=39)
                                                                                                                                          cggggCCGCTCTGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGAc              <  2:4582743/147‑1 (MQ=255)
                                                                                                                                          cggggCCGCTCTGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGAc              >  1:4582743/3‑149 (MQ=255)
                                                                                                                                            gggCCGCTCTGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTc           >  1:2644485/1‑150 (MQ=255)
                                                                                                                                            gggCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTc           <  1:1174109/150‑1 (MQ=255)
                                                                                                                                                cGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAAc       <  2:6916765/150‑1 (MQ=255)
                                                                                                                                                cGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAAc       <  2:953641/150‑1 (MQ=255)
                                                                                                                                                 gCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAAcc      >  1:6216047/1‑150 (MQ=255)
                                                                                                                                                   agagCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAACCCg    <  1:995225/150‑1 (MQ=255)
                                                                                                                                                   agagCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAACCCg    <  1:8272004/150‑1 (MQ=255)
                                                                                                                                                    gagCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAACCCgg   <  1:4013636/150‑1 (MQ=255)
                                                                                                                                                    gagCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGCTCCGCGCCAGATGCACATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAACCCgg   >  2:1316931/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GACCCTTCTCTCTCCTGCGGCCGCATTTTTTTGACAGTCATGCGGCGAATGACCTACCTAAAGCAACAAAGGTGCCAGCCTTCGCAACTCATTTCCCAGAGATCACACAATTCCTGCGGGAGCGGCCTTGCGTCGCGATAGGGCCGCAGAGCGGCCCCGGCGATCTTTGCGGCGGAGCTGAAGACCTGGGGCCGCTTCGCGCCCCATTCGCGACACAAGGCCGTTCCCACAGAGGTCCGCGCCAGTTGCAGATGATTTGCCAGTCAGCCGAACAGGCTGGCCTCGACGTCAAACCCGGC  >  NC_002947/4965813‑4966111

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: