Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,947,358 T→A 6.9% intergenic (+44/‑161) betX → / → yjcH choline/betaine/carnitine ABC transporter substrate‑binding protein BetX/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,947,3580TA 6.9% 634.8 / 22.5 219intergenic (+44/‑161)betX/yjcHcholine/betaine/carnitine ABC transporter substrate‑binding protein BetX/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (110/88);  new base A (6/8);  total (116/103)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 219 total aligned reads displayed.

GGCGAGGTGATGCTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGTGCTGCT  >  NC_002947/1947210‑1947506
                                                                                                                                                     |                                                                                                                                                    
ggCGAGGTGATGCTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTc                                                                                                                                                     >  2:488403/1‑150 (MQ=255)
     ggTGATGCTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCAcccgta                                                                                                                                                <  2:6404028/150‑5 (MQ=255)
      gtgtTGGTGGCGATGCATGTAGGGGCCAATCCTGAAGCGGCCGGCAAGGTCTTGGTCGGCAAGGAACCTTCTCGTGTTTAGGATTTGTATGGCAATTAAAAGCCTTCCAGGGCCAATGTTTTTTCTTCGGGGGTTAACCCGCTCCTACgg                                                                                                                                               <  2:3879175/150‑1 (MQ=255)
      gTGATGCTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACgg                                                                                                                                               <  2:3493393/150‑1 (MQ=255)
      gTGATGCTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACgg                                                                                                                                               <  2:6429670/150‑1 (MQ=255)
       tGATGCTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACggg                                                                                                                                              >  1:1088213/1‑150 (MQ=255)
        gATGCTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGt                                                                                                                                             <  2:1501355/150‑1 (MQ=255)
        gATGCTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGt                                                                                                                                             <  1:5342915/150‑1 (MQ=255)
        gATGCTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGt                                                                                                                                             <  1:5459689/150‑1 (MQ=255)
        gATGCTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCAcccgtagga                                                                                                                                             <  1:5665551/150‑8 (MQ=255)
           gCTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCg                                                                                                                                          <  1:4814205/150‑1 (MQ=255)
            cTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCgg                                                                                                                                         <  1:6014109/150‑1 (MQ=255)
            cTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCgg                                                                                                                                         <  1:5311305/150‑1 (MQ=255)
            cTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCgg                                                                                                                                         <  1:2551627/150‑1 (MQ=255)
              ggCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGt                                                                                                                                       <  2:5834042/150‑1 (MQ=255)
                 gATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTa                                                                                                                                    <  1:5575829/150‑1 (MQ=255)
                 gATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTa                                                                                                                                    >  2:806534/1‑150 (MQ=255)
                  aTCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTAc                                                                                                                                   <  2:6465366/150‑1 (MQ=255)
                  aTCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCACCCGTAGGAGCGGGTTTAc                                                                                                                                   <  2:5588408/150‑1 (MQ=39)
                    ccAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTAccc                                                                                                                                 >  2:1191184/1‑150 (MQ=255)
                     cAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCcg                                                                                                                                >  1:5769490/1‑150 (MQ=255)
                     cAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCcg                                                                                                                                >  1:1217874/1‑150 (MQ=255)
                      aGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCcgc                                                                                                                               >  2:1908701/1‑150 (MQ=255)
                      aGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCcgc                                                                                                                               >  2:3912868/1‑150 (MQ=255)
                      aGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCcgc                                                                                                                               >  1:4097507/1‑150 (MQ=255)
                       ggAAGGGGACTATCCTGATGCGGCCGCCTCGGACTGGGTCGCCAAGCACCGTGTTCTTGTGAAGGAGTGGACTGGAAATTAAATGCCTTCCTGGGCCAATGTTGTTTCTTGGGGGGTAAACCCGCTCATCCTGGTGCGGGTTTACCcgcg                                                                                                                              <  2:4630797/150‑1 (MQ=255)
                        gAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGa                                                                                                                             >  1:2767238/1‑150 (MQ=255)
                        gAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGa                                                                                                                             >  1:998151/1‑150 (MQ=255)
                         aaGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGaa                                                                                                                            <  2:4952638/150‑1 (MQ=255)
                              gCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGt                                                                                                                       >  1:5986699/1‑150 (MQ=255)
                              gCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGGTTACCCGCGAACGGGt                                                                                                                       >  2:17190/1‑150 (MQ=255)
                                     cTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACt                                                                                                                >  1:4749204/1‑150 (MQ=255)
                                      tGACGCGGCCGCCACGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAACGCCTTCCCGGGCCACTGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACtt                                                                                                               <  2:3345384/150‑1 (MQ=39)
                                      tGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACtt                                                                                                               >  1:866136/1‑150 (MQ=255)
                                          gCGGCCGCCACGGCCTGGGTCGCCAAGCACCCTGATGGTGTGAAGGCGTGGACTGCAAAGTAAACGCCTTCCCGGGCCAATGTTTTTTCTTCGGGGTTAAACCCGCACCTAGGGGTGCGTTTTTACCCGCGAACGGGTCATCACTTTGtt                                                                                                           <  2:1963864/150‑1 (MQ=255)
                                          gCGGCCGCCACGGACTGGGGCGCCAAGCACCCTGATGGTGTGACGGAGTGGAGTGGCAAGTAAACGCCGTCCTGGGCCAATGTTTTTTCTTCGCGGGTAAACCCGCTCCTAGGGGTGGGGGTTTACCCGCGAACGGGTCATCACTTTGtt                                                                                                           <  2:6629670/150‑1 (MQ=255)
                                               cgccAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGATGTTTCTTCGCGGGTAAACCCGCTCCTAAGGGTGCGGGTTTACCCGCGAACGGGTCATCACATTGTTACGCa                                                                                                      >  1:4642060/1‑150 (MQ=255)
                                                  cAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGCCATCACTTTGTTACGCAATc                                                                                                   >  1:692289/1‑150 (MQ=255)
                                                   aaGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAAACg                                                                                                  >  1:4574641/1‑150 (MQ=255)
                                                         tGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTaaaa                                                                                            >  1:4313019/1‑150 (MQ=255)
                                                               gCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGt                                                                                      >  2:2809201/1‑150 (MQ=255)
                                                                ccATGCACCTTGTTGGTTTTATGGAGTGGCCTGCCAATAAAAAGCCTGCCCGGGCAATTGTTGTTTCTTGGCGGGAAAACCCGCCCCTCCGGGGGGGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGtt                                                                                     <  1:6728151/150‑1 (MQ=255)
                                                                ccAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCATGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTCCGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGtt                                                                                     <  2:5369112/150‑1 (MQ=255)
                                                                ccAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTCCCTGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGtt                                                                                     <  2:512713/150‑1 (MQ=255)
                                                                        ccTGATGGTGTGAAGGTGTGGACTGGCATGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTCCGGGTGGGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTaa                                                                             <  2:1025014/150‑1 (MQ=255)
                                                                        ccTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTaa                                                                             <  2:5535696/150‑1 (MQ=255)
                                                                              tctgtgATGGTGTGGACGGGCATGTAAAAGCCTTCCAGGGCCACTGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGGGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTaa                                                                       <  2:1564006/148‑1 (MQ=39)
                                                                               gtgtGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAg                                                                      <  2:2767238/150‑1 (MQ=255)
                                                                                  tgAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAgaga                                                                   >  1:1564006/1‑147 (MQ=255)
                                                                                       gAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTg                                                              >  2:699107/1‑150 (MQ=255)
                                                                                       gAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTg                                                              >  2:2999207/1‑150 (MQ=255)
                                                                                                ggCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTc                                                     <  2:2760207/150‑1 (MQ=255)
                                                                                                 gCAAGTAAACGCCTCCCAGGGCCACTGTTGTTTCTTCGCGGGTAACCCCGCTCCTCCGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTcc                                                    <  2:2392532/150‑1 (MQ=255)
                                                                                                 gCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTcc                                                    >  2:617979/1‑150 (MQ=255)
                                                                                                 gCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGGCTGGTTGGCGTcc                                                    >  1:632079/1‑150 (MQ=255)
                                                                                                   aaGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCAt                                                  >  1:5059980/1‑150 (MQ=255)
                                                                                                       aaaaGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGgc                                              >  2:966013/1‑150 (MQ=255)
                                                                                                        aaaGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGgca                                             >  1:749171/1‑150 (MQ=255)
                                                                                                        aaaGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGgca                                             >  2:5868202/1‑150 (MQ=255)
                                                                                                         aaGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGgcag                                            >  1:3744677/1‑150 (MQ=255)
                                                                                                            ccTCCCCGGGCCATTGTTGTTTCTCCGCGGCTAAACCCGCCCCCCCGGGGGGGGTTTAACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCAt                                         <  1:3205411/150‑1 (MQ=255)
                                                                                                            ccTCCCCGGGCAATTGTTTTTTCTTCGGGGGTAACCCCGCTCCTCGGGGTGGGGGTTACCCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCAt                                         <  1:388637/150‑1 (MQ=255)
                                                                                                               tCCTGGGCCATTGTTTTTTCTTCGCGGTTAAACCCGCTCCTCCGGGTGGGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATaaa                                      <  2:3156902/150‑1 (MQ=255)
                                                                                                                cctggcccATTGTTGTTTCTGCGCGGTAACACCCCCTCCTCGGGTTGCGGGTTCACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAg                                     <  1:6325803/144‑1 (MQ=255)
                                                                                                                  aGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTa                                   >  1:2008525/1‑150 (MQ=255)
                                                                                                                   gggCAAATTTTGTTTCTTCGCGGGTAAACCCGCTCCTCCGGGTGCGGTTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTaa                                  <  2:4306324/150‑1 (MQ=255)
                                                                                                                     gcctttgttgtTTCTTCGCGGTTAAACCCGCTCCTAGGGGTGCGGTTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAgg                                <  2:2272786/145‑1 (MQ=255)
                                                                                                                     gactttgttttTTCTGCGGGTTAAACCCCGCTCCTCGGGGTGGGTTTTTCCCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAgg                                <  1:3536466/145‑1 (MQ=255)
                                                                                                                     gCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAgg                                <  2:5622529/150‑1 (MQ=255)
                                                                                                                     cacatttttgtttgtTGGGGGGAAACCCCGCGCCTGGGGGTGGGGTTTAACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAgg                                <  1:3683353/143‑1 (MQ=255)
                                                                                                                       atttgttttTTCTTGGGGGTTAAACCCGCTCCTCGGGGTGGGGTTTAACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGgt                              <  1:1359462/147‑1 (MQ=255)
                                                                                                                         aTGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGgtgt                            >  2:4747488/1‑150 (MQ=255)
                                                                                                                             tgttTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGt                        >  2:6801537/1‑150 (MQ=255)
                                                                                                                              gttTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGtt                       >  2:914756/1‑150 (MQ=255)
                                                                                                                               tttCTTCGCGGGTAAACCCGCTCCTACGGGTGGGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTc                      <  2:1505786/150‑1 (MQ=255)
                                                                                                                               tttCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTc                      <  2:4415488/150‑1 (MQ=255)
                                                                                                                                ccgttcGCGGGTAAACCCGCACCCGTAGGAGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTcc                     >  2:6116940/4‑150 (MQ=38)
                                                                                                                                       cgGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGt              >  2:492959/1‑150 (MQ=255)
                                                                                                                                        gggTAACCCCGCTCCTCGGGGTGGGGTTTAACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTa             <  1:4134504/150‑1 (MQ=255)
                                                                                                                                        gggTAAACCCGCACCCGTAGGAGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTa             >  1:2252319/1‑150 (MQ=39)
                                                                                                                                           tAAACCCGCACCCGTAGGAGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCt          >  2:108243/1‑150 (MQ=39)
                                                                                                                                             aaCCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCtgt        >  1:4426417/1‑150 (MQ=255)
                                                                                                                                             aaCCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCtgt        >  2:6283796/1‑150 (MQ=255)
                                                                                                                                             aaCCCGCCCCTACGGGTGCGGGTTCACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCtgt        <  2:6011636/150‑1 (MQ=255)
                                                                                                                                              aCCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCtgtg       >  1:4507713/1‑150 (MQ=255)
                                                                                                                                                cccGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTaa                                                                             <  1:1238405/79‑1 (MQ=255)
                                                                                                                                                cccGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGtgc      >  1:4175794/1‑150 (MQ=255)
                                                                                                                                                cccGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGtgc      >  2:4128437/1‑150 (MQ=255)
                                                                                                                                                cccGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGtgc      >  2:6734582/1‑150 (MQ=255)
                                                                                                                                                cccGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGtgc      >  1:4474961/1‑150 (MQ=255)
                                                                                                                                                cccGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGtgc      >  2:5738852/1‑150 (MQ=255)
                                                                                                                                                 ccGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGtgct     >  2:3172129/1‑150 (MQ=255)
                                                                                                                                                 ccGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGtgct     >  1:1745314/1‑150 (MQ=255)
                                                                                                                                                 ccGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGtgct     >  2:5134009/1‑150 (MQ=255)
                                                                                                                                                 ccGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGtgct     >  2:3732155/1‑150 (MQ=255)
                                                                                                                                                 ccGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGtgct     >  2:2450958/1‑150 (MQ=255)
                                                                                                                                                   gCTCCTCGGGGTGGGGTTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGtgctgc   <  1:5407205/150‑1 (MQ=255)
                                                                                                                                                   gCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCg                                                >  2:1600896/1‑105 (MQ=255)
                                                                                                                                                   gCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGtgcggc   >  1:2431141/1‑150 (MQ=255)
                                                                                                                                                    cTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGTgctgct  <  1:6141772/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGCGAGGTGATGCTGGCGATCCAGGAAGGGGCCAAGCCTGAAGCGGCCGCCAAGGACTGGGTCGCCAAGCACCCTGATCGTGTGAAGGAGTGGACTGGCAAGTAAAAGCCTTCCAGGGCCAATGTTGTTTCTTCGCGGGTAAACCCGCTCCTACGGGTGCGGGTTTACCCGCGAACGGGTCATCACTTTGTTACGCAATCGGTAAAACACCGTTGCATCTAAGACTAAGGTCGTCTGGTTGGCGTCCATCGGCAGCATAAAGTAAGGGTGTGGGTTCCATCCCGTATCTGTGCTGCT  >  NC_002947/1947210‑1947506

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: