Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,481,522 T→C 5.3% S75G (AGC→GGC)  rhlB ← ATP‑dependent RNA helicase RhlB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,481,5220TC 5.3% 547.6 / 11.5 171S75G (AGC→GGC) rhlBATP‑dependent RNA helicase RhlB
Reads supporting (aligned to +/- strand):  ref base T (92/69);  new base C (3/6);  total (95/76)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 173 total aligned reads displayed.

AACCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGCGTCA  >  NC_002947/1481375‑1481671
                                                                                                                                                   |                                                                                                                                                     
 aCCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGc                                                                                                                                                    >  2:3545122/1‑150 (MQ=255)
  cccAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGcc                                                                                                                                                   <  1:2428600/150‑1 (MQ=255)
  cccAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGcc                                                                                                                                                   <  2:4792751/150‑1 (MQ=255)
   ccAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCt                                                                                                                                                  >  2:1075032/1‑150 (MQ=255)
   ccAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCt                                                                                                                                                  >  1:2638148/1‑150 (MQ=255)
   ccAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCt                                                                                                                                                  >  1:321495/1‑150 (MQ=255)
   ccAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCt                                                                                                                                                  >  1:3460338/1‑150 (MQ=255)
   ccAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCt                                                                                                                                                  <  2:582916/150‑1 (MQ=255)
      ggACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCCAGCCtgtg                                                                                                                                               <  1:2070867/150‑1 (MQ=255)
        aCCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCCAGCCTGTgga                                                                                                                                             <  2:2208915/150‑2 (MQ=255)
         ccTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCCTGCCTGTggta                                                                                                                                            <  1:427244/150‑3 (MQ=255)
         ccTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCCAGCCTGTggaa                                                                                                                                            <  1:3720420/150‑3 (MQ=255)
           tGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttg                                                                                                                                          <  2:4179010/150‑1 (MQ=255)
           tGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttg                                                                                                                                          <  2:4178715/150‑1 (MQ=255)
            gcgcCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttgg                                                                                                                                         >  2:4080909/1‑150 (MQ=255)
             cgcCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttggg                                                                                                                                        <  2:2603316/150‑1 (MQ=255)
              gcCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCGTCCACAGGCTGGCCTGTggcttgggc                                                                                                                                       >  2:2488826/1‑150 (MQ=255)
              gcCTGGATCGGCGAGCAGTACGGGAAGCCCAGGTCGTGGATGGAGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttgggc                                                                                                                                       <  2:4642419/150‑1 (MQ=255)
               ccTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTggcttgggct                                                                                                                                      >  2:1191950/1‑150 (MQ=255)
                     tCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTGTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGtt                                                                                                                                >  1:305754/1‑150 (MQ=255)
                     tCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGtt                                                                                                                                >  2:4748663/1‑150 (MQ=255)
                      cGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTc                                                                                                                               <  1:5069674/150‑1 (MQ=255)
                        gCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTcgc                                                                                                                             >  1:1622446/1‑150 (MQ=255)
                             cAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTcgcggcgc                                                                                                                        <  2:4204334/150‑1 (MQ=255)
                                  cGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGttt                                                                                                                   <  2:1153846/150‑1 (MQ=255)
                                     gAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTc                                                                                                                <  1:2474807/150‑1 (MQ=255)
                                     gAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTc                                                                                                                >  1:3764903/1‑150 (MQ=255)
                                       aGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCt                                                                                                              <  2:156569/150‑1 (MQ=255)
                                       aGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCt                                                                                                              <  2:3233440/150‑1 (MQ=255)
                                        gCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCtt                                                                                                             <  1:1535257/150‑1 (MQ=255)
                                         cccAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTg                                                                                                            <  2:2073679/150‑1 (MQ=255)
                                         cccAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTg                                                                                                            <  2:1393478/150‑1 (MQ=255)
                                          ccAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTgg                                                                                                           >  1:1838821/1‑150 (MQ=255)
                                               tCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGt                                                                                                      <  1:1268985/150‑1 (MQ=255)
                                               tCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGt                                                                                                      >  1:3693652/1‑150 (MQ=255)
                                                 gTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCGTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTg                                                                                                    >  2:2855045/1‑150 (MQ=255)
                                                      tGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCgg                                                                                               >  2:3837389/1‑150 (MQ=255)
                                                        gCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCa                                                                                             <  1:2356935/150‑1 (MQ=255)
                                                                 aGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTc                                                                                    <  2:4263525/150‑1 (MQ=255)
                                                                  gCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTcg                                                                                   >  2:1226232/1‑150 (MQ=255)
                                                                   cTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTcgc                                                                                  >  2:3936624/1‑150 (MQ=255)
                                                                   cTCGTTGGACAACTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTcgc                                                                                  <  2:543902/150‑1 (MQ=255)
                                                                      gTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCgcggc                                                                               <  1:850475/150‑1 (MQ=255)
                                                                      gTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCgcggc                                                                               <  2:1037274/150‑1 (MQ=255)
                                                                       ttGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGAt                                                                                                  >  2:3755023/1‑130 (MQ=255)
                                                                          gACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCa                                                                           >  1:4837762/1‑150 (MQ=255)
                                                                          gACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGCCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCa                                                                           >  2:182608/1‑150 (MQ=255)
                                                                           aCAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAg                                                                          <  2:625308/150‑1 (MQ=255)
                                                                               cTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTg                                                                      >  2:1846297/1‑150 (MQ=255)
                                                                                ttGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGc                                                                                                                                   >  2:1994910/1‑88 (MQ=255)
                                                                                ttGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGc                                                                                                                                   <  1:1994910/88‑1 (MQ=255)
                                                                                ttGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTg                                                                                         <  2:3632368/130‑1 (MQ=255)
                                                                                   aaGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCa                                                                  <  1:3479111/150‑1 (MQ=255)
                                                                                    aGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAg                                                                 <  1:2904714/150‑1 (MQ=255)
                                                                                        gTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCtt                                                             >  1:4286371/1‑150 (MQ=255)
                                                                                         tGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTg                                                            >  1:4606053/1‑150 (MQ=255)
                                                                                         tGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTg                                                            >  1:3567518/1‑150 (MQ=255)
                                                                                         gggAAGCGGGTTTTGCCGTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTg                                                            <  2:2915073/149‑1 (MQ=255)
                                                                                              gCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGc                                                       >  1:1006692/1‑150 (MQ=255)
                                                                                              gCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGc                                                       >  1:1006689/1‑150 (MQ=255)
                                                                                               cGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCg                                                      <  2:4614635/150‑1 (MQ=255)
                                                                                                gggTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCgg                                                     <  2:3304697/150‑1 (MQ=255)
                                                                                                gggTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCgg                                                     >  2:326595/1‑150 (MQ=255)
                                                                                                gggTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCgg                                                     >  1:3473615/1‑150 (MQ=255)
                                                                                                  gTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCgggg                                                   <  1:3639367/150‑1 (MQ=255)
                                                                                                   ttttGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGt                                                  >  1:1357750/1‑150 (MQ=255)
                                                                                                   ttttGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGt                                                  >  1:3222714/1‑150 (MQ=255)
                                                                                                    tttGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGtt                                                 <  2:3565526/150‑1 (MQ=255)
                                                                                                      tGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTgg                                               <  1:5094845/150‑1 (MQ=255)
                                                                                                          ttCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATg                                           >  1:209469/1‑150 (MQ=255)
                                                                                                          ttCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATg                                           >  2:2506666/1‑150 (MQ=255)
                                                                                                           tCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATgg                                          <  2:4353021/150‑1 (MQ=255)
                                                                                                            ccTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATggg                                         <  1:2193360/150‑1 (MQ=255)
                                                                                                            ccTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATggg                                         >  1:1927317/1‑150 (MQ=255)
                                                                                                            ccTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATggg                                         <  1:1191950/150‑1 (MQ=255)
                                                                                                              tGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGc                                                                                                          >  1:4016771/1‑83 (MQ=255)
                                                                                                                  gCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCgg                                   <  2:3071873/150‑1 (MQ=255)
                                                                                                                    tCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCGTGGGGGCGGGGTTGGCATGGGGTGCGGcc                                 >  1:1316031/1‑150 (MQ=255)
                                                                                                                         cACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGttt                            <  2:805185/150‑1 (MQ=255)
                                                                                                                         aaCGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGttt                            >  2:2959765/2‑150 (MQ=255)
                                                                                                                          aCGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGtttt                           >  1:4592551/1‑150 (MQ=255)
                                                                                                                           cGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTg                          >  1:4110078/1‑150 (MQ=255)
                                                                                                                              aGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCcg                       >  1:2221391/1‑150 (MQ=255)
                                                                                                                                  ttCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCtt                                                             >  1:2251341/1‑108 (MQ=255)
                                                                                                                                  ttCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCtt                                                             <  2:2251341/108‑1 (MQ=255)
                                                                                                                                  ttCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTGCCGcgggc                  >  1:4704556/1‑150 (MQ=255)
                                                                                                                                    cGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGcgggcc                 <  1:3829099/150‑1 (MQ=255)
                                                                                                                                    cGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGcgggcc                 >  2:4587055/1‑150 (MQ=255)
                                                                                                                                     gggCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGcgggccg                >  1:2411615/1‑150 (MQ=255)
                                                                                                                                       caagccACAGGCCAGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGcgggccggg              >  1:3857465/5‑150 (MQ=255)
                                                                                                                                         ttCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCt            >  1:3257570/1‑150 (MQ=255)
                                                                                                                                            cacaGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTCGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTcc         >  1:1708577/1‑150 (MQ=255)
                                                                                                                                             acaGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCg        <  2:4964478/150‑1 (MQ=255)
                                                                                                                                             acaGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCg        >  1:902887/1‑150 (MQ=255)
                                                                                                                                             acaGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCg        >  2:3260966/1‑150 (MQ=255)
                                                                                                                                               aGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGc      <  2:1480675/150‑1 (MQ=255)
                                                                                                                                                 gCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGCGt    <  1:4842958/150‑1 (MQ=255)
                                                                                                                                                 gCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGCGt    <  1:3034202/150‑1 (MQ=255)
                                                                                                                                                   tGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGCGTCa  >  1:4716580/1‑150 (MQ=255)
                                                                                                                                                   tGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGCGTCa  >  2:3729606/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AACCCAGGACCTGCGCCTGGATCGGCGTGCAGTACGGGAAGCCCAGGTCGTGGATGGCGTGCATCAGCTCGTTGGACAGCTTGAAGTCGTGGAAGCGGGTTTTGCCTTCCTGCGGCTCGACCACGAAGTCTTCGGGCTTCCACAGGCTGGCCTGTGGCTTGGGCTTGCGTTCGCGGCGCGGTTTGTCCTTGGCCGGTTGATCGGCAGGTGGCGTCGCGGCGGCAGGTTGCTCAGCCTTGGGGGCGGGGTTGGCATGGGGTGCGGCCGGTTTTGCCGCGGGCCGGGCTTCCGGCGTCA  >  NC_002947/1481375‑1481671

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: