Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,576,764 A→G 5.0% F168S (TTC→TCC)  motB ← flagellar motor rotation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,576,7640AG 5.0% 1040.7 / 25.0 318F168S (TTC→TCC) motBflagellar motor rotation protein
Reads supporting (aligned to +/- strand):  ref base A (119/183);  new base G (11/5);  total (130/188)
Fisher's exact test for biased strand distribution p-value = 3.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 326 total aligned reads displayed.

CCAGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAGACGCTCGCGC  >  NC_002947/5576616‑5576913
                                                                                                                                                    |                                                                                                                                                     
 cAGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAAc                                                                                                                                                     <  2:8747612/150‑1 (MQ=255)
 cAGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAAc                                                                                                                                                     >  1:1409511/1‑150 (MQ=255)
  aGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGaa                                                                                                                                                      <  2:91172/148‑1 (MQ=255)
  aGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACa                                                                                                                                                    >  1:7657670/1‑150 (MQ=255)
       ccGTCTTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACAtcggt                                                                                                                                               <  2:4174820/150‑1 (MQ=255)
         gTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCCGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACAtcggtcg                                                                                                                                             >  1:3033107/1‑150 (MQ=255)
         gTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACAtcggtcg                                                                                                                                             >  1:296134/1‑150 (MQ=255)
           aTTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATcggtcggt                                                                                                                                           >  1:5257334/1‑150 (MQ=255)
               gCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTtctc                                                                                                                                       >  1:1177616/1‑150 (MQ=255)
               gCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTtctc                                                                                                                                       >  2:7343331/1‑150 (MQ=255)
               gCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTtctc                                                                                                                                       <  1:2330979/150‑1 (MQ=255)
               gCCAGTGCCGGCATACGACTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTtctc                                                                                                                                       <  1:3509054/150‑1 (MQ=255)
                ccAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCg                                                                                                                                      <  2:4984751/150‑1 (MQ=255)
                 cAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCgg                                                                                                                                     <  2:3399812/150‑1 (MQ=255)
                        ggCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCAt                                                                                                                              <  1:2451629/150‑1 (MQ=255)
                             aCGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGTCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCt                                                                                                                         <  1:4939111/150‑1 (MQ=255)
                             aCGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCt                                                                                                                         <  2:7640896/150‑1 (MQ=255)
                             aCGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCt                                                                                                                         <  1:2359049/150‑1 (MQ=255)
                              cGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTg                                                                                                                        <  1:6824135/150‑1 (MQ=255)
                                gCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGa                                                                                                                      <  2:199583/150‑1 (MQ=255)
                                 cTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGAt                                                                                                                     <  2:6348156/150‑1 (MQ=255)
                                      cATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCAGCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCa                                                                                                                <  1:5547025/150‑1 (MQ=255)
                                      cATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCa                                                                                                                <  1:6192321/150‑1 (MQ=255)
                                       aTCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAg                                                                                                               >  1:5468111/1‑150 (MQ=255)
                                       aTCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAg                                                                                                               >  2:2531498/1‑150 (MQ=255)
                                               ggCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCt                                                                                                       <  1:3255530/150‑1 (MQ=255)
                                                 ccGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGa                                                                                                     <  1:8856511/150‑1 (MQ=255)
                                                 ccGCTGATGCTGACCTTCTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGa                                                                                                     <  1:7514953/150‑1 (MQ=255)
                                                       aTGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGGAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGAtc                                                                                               <  2:6839286/150‑1 (MQ=255)
                                                         gCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGAtctc                                                                                             <  2:56018/150‑1 (MQ=255)
                                                          cTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCg                                                                                            >  2:7585128/1‑150 (MQ=255)
                                                            gACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGaa                                                                                          >  2:6952176/1‑150 (MQ=255)
                                                             aCCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAAc                                                                                         <  2:2274975/150‑1 (MQ=255)
                                                             aCCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAAc                                                                                         <  1:7580339/150‑1 (MQ=255)
                                                              ccTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACa                                                                                        <  1:6536255/150‑1 (MQ=255)
                                                                  gTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGAt                                                                                    <  1:5916503/150‑1 (MQ=255)
                                                                   ttCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATc                                                                                   >  1:6066117/1‑150 (MQ=255)
                                                                   ttCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATc                                                                                   >  2:1419401/1‑150 (MQ=255)
                                                                      ggCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTgg                                                                                <  2:1907554/150‑1 (MQ=255)
                                                                      ggCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTgg                                                                                <  2:8169374/150‑1 (MQ=255)
                                                                      ggCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTgg                                                                                <  1:1042357/150‑1 (MQ=255)
                                                                      ggCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTgg                                                                                <  1:360144/150‑1 (MQ=255)
                                                                       gCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGt                                                                               >  2:1720261/1‑150 (MQ=255)
                                                                        cACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTc                                                                              <  2:2927253/150‑1 (MQ=255)
                                                                        cACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTc                                                                              >  1:5199991/1‑150 (MQ=255)
                                                                         aCCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTcc                                                                             >  1:3981468/1‑150 (MQ=255)
                                                                          ccGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCt                                                                            >  2:8990433/1‑150 (MQ=255)
                                                                            gCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTg                                                                          <  1:8139164/150‑1 (MQ=255)
                                                                             ccTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGa                                                                         <  1:8518611/150‑1 (MQ=255)
                                                                                   tGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCt                                                                   <  2:8466608/150‑1 (MQ=255)
                                                                                   tGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCt                                                                   <  1:3734421/150‑1 (MQ=255)
                                                                                   tGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCt                                                                   <  1:5569051/150‑1 (MQ=255)
                                                                                      tgtCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCa                                                                <  1:2574636/150‑1 (MQ=255)
                                                                                             cATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCgg                                                         <  1:4627888/150‑1 (MQ=255)
                                                                                             cATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCgg                                                         <  1:964353/150‑1 (MQ=255)
                                                                                             cATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCgg                                                         <  1:7892765/150‑1 (MQ=255)
                                                                                              aTGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGAGACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCggg                                                        <  2:5159349/150‑1 (MQ=255)
                                                                                                 gCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTtc                                                     <  2:6377679/150‑1 (MQ=255)
                                                                                                    agcagGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTctct                                                  >  2:3015017/1‑150 (MQ=255)
                                                                                                     gcagGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCtt                                                 <  2:3227810/150‑1 (MQ=255)
                                                                                                     gcagGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCtt                                                 <  2:3891219/150‑1 (MQ=255)
                                                                                                      cagGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTTAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTc                                                >  2:5956403/1‑150 (MQ=255)
                                                                                                        ggATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCa                                              <  2:276679/150‑1 (MQ=255)
                                                                                                            atCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCtt                                          <  1:8363642/150‑1 (MQ=255)
                                                                                                             tCTTCGAAGCAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTg                                         <  1:8939201/150‑1 (MQ=255)
                                                                                                               ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTGGGTTCTCGGCATCCATGATCTGGATAGGCAGGCCGTCGTGAGTGATCTCGAACAGGATCTGGTCGTTGAAGTTGTGCAGCTGCGGGTTGTCTTCCTCCTTGtt                                       >  2:4813577/1‑150 (MQ=255)
                                                                                                               ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCgg                                                         >  1:641045/1‑132 (MQ=255)
                                                                                                               ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGtt                                       <  2:1394105/150‑1 (MQ=255)
                                                                                                               ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGtt                                       >  2:1229601/1‑150 (MQ=255)
                                                                                                               ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGtt                                       >  2:5821240/1‑150 (MQ=255)
                                                                                                               ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGtt                                       >  1:4482959/1‑150 (MQ=255)
                                                                                                               ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGtt                                       >  1:5665193/1‑150 (MQ=255)
                                                                                                               ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGtt                                       <  2:3852481/150‑1 (MQ=255)
                                                                                                               ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGtt                                       >  1:4444229/1‑150 (MQ=255)
                                                                                                               ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGtt                                       >  2:346558/1‑150 (MQ=255)
                                                                                                               ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGtt                                       >  1:1342199/1‑150 (MQ=255)
                                                                                                               ttCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGtt                                       >  1:3852536/1‑150 (MQ=255)
                                                                                                                  gAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTg                                    >  1:2869947/1‑150 (MQ=255)
                                                                                                                    aGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAGCTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCa                                  <  2:6652691/150‑1 (MQ=255)
                                                                                                                    aGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCa                                  <  2:4309812/150‑1 (MQ=255)
                                                                                                                    aGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCa                                  <  1:6851943/150‑1 (MQ=255)
                                                                                                                        gggCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTc                              <  1:7344685/150‑1 (MQ=255)
                                                                                                                        gggCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTc                              <  1:5059777/150‑1 (MQ=255)
                                                                                                                         ggCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTcc                             <  2:3537171/150‑1 (MQ=255)
                                                                                                                          gCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTcat                            >  2:6535808/1‑148 (MQ=255)
                                                                                                                            tGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTgc                          <  1:4690308/150‑1 (MQ=255)
                                                                                                                                 gcgcgcgcTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCaa                     <  2:296134/150‑1 (MQ=255)
                                                                                                                                       gcTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTc               <  1:4177322/150‑1 (MQ=255)
                                                                                                                                        cgaCCGATGTTCGACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTcc              >  1:7120930/3‑150 (MQ=255)
                                                                                                                                         tACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGACCTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCACCAGCTCCa             >  2:5494906/1‑150 (MQ=255)
                                                                                                                                         tACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCa             <  1:484241/150‑1 (MQ=255)
                                                                                                                                         gaCCGATGTTCGACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCa             >  2:7780104/2‑150 (MQ=255)
                                                                                                                                         gaCCGATGTTCGACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCa             >  2:4575114/2‑150 (MQ=255)
                                                                                                                                          aCCGATGTTCGACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAg            >  2:7561861/1‑150 (MQ=255)
                                                                                                                                          aCCGATGTTCGACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAg            >  2:8441492/1‑150 (MQ=255)
                                                                                                                                          aCCGATGTTCGACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAg            >  2:5533420/1‑150 (MQ=255)
                                                                                                                                            cGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAGAc          >  2:212061/1‑150 (MQ=255)
                                                                                                                                             gATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAGAcg         <  2:1815293/150‑1 (MQ=255)
                                                                                                                                                  cGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAGAcgctcgc    >  2:8019231/1‑150 (MQ=255)
                                                                                                                                                   gAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAGAcgctcgcg   >  1:859969/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCAGTTACCGTATTCGCCAGTGCCGGCATACGGCTTGGCATCGGTATGGCCGCTGATGCTGACCTTGTTCGGCACCGCCTTGATGGTGTCAGCCATGGCCAGCAGGATATCTTCGAAGTAGGGCTGCAGGCGCGCGCTACCGATGTCGAACATCGGTCGGTTCTCGGCATCCATGATCTGGATACGCAGGCCGTCCTGAGTGATCTCGAACAGGATCTGGTCCTTGAACTTCTGCAGCTGCGGGTTCTCTTCCACCTTGTTCTGCAGTTCCTGCAGCAACAGCTCCAGACGCTCGCGC  >  NC_002947/5576616‑5576913

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