Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,629,486 G→A 8.6% Q318Q (CAG→CAA mucB → sigma factor AlgU regulator MucB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,629,4860GA 8.6% 189.7 / 7.6 70Q318Q (CAG→CAAmucBsigma factor AlgU regulator MucB
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base A (2/4);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

ACCGACATGCGCACCCAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTTG  >  NC_002947/1629337‑1629630
                                                                                                                                                     |                                                                                                                                                
aCCGACATGCGCACCCAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAg                                                                                                                                                  >  1:95947/1‑150 (MQ=255)
    aCATGCGCACCCAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGgcgc                                                                                                                                              >  1:624618/1‑150 (MQ=255)
    aCATGCGCACCCAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGgcgc                                                                                                                                              <  1:1771198/150‑1 (MQ=255)
      aTGCGCACCCAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGgcgcgc                                                                                                                                            >  1:635463/1‑150 (MQ=255)
         cgcACCCAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAg                                                                                                                                         <  1:361465/150‑1 (MQ=255)
              ccAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATg                                                                                                                                    >  2:1234373/1‑150 (MQ=255)
              ccAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATg                                                                                                                                    <  2:476133/150‑1 (MQ=255)
              ccAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATg                                                                                                                                    >  1:476133/1‑150 (MQ=255)
                aGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGAGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGgc                                                                                                                                  >  1:603052/1‑150 (MQ=255)
                  cTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGgcgc                                                                                                                                <  2:1868590/150‑1 (MQ=255)
                        ccGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAAGCGCGCCATTACTGgcgcgcctgg                                                                                                                          <  2:795188/150‑5 (MQ=255)
                        ccGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAAGCGCGCCATTACTGgcgcgcctgg                                                                                                                          <  2:435909/150‑5 (MQ=255)
                            ccTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGt                                                                                                                      <  2:603052/150‑1 (MQ=255)
                             ctcatcAGGCTCTCGGCGCCTAACTACACCCAAGGCCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGCACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGAGCGCCAGTAATGGCGCGCTTTTGCGTg                                                                                                                     >  2:1013772/6‑150 (MQ=255)
                             cTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAAGCGCGCCATTACTGgcgcgcctg                                                                                                                           >  2:599244/1‑141 (MQ=255)
                             cTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAAGCGCGCCATTACTGgcgcgcctg                                                                                                                           <  1:599244/144‑4 (MQ=255)
                                     tctcGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCa                                                                                                             >  2:1181894/1‑150 (MQ=255)
                                      ctcGCGTCGCCTGAATACCCCCAAGGGAAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCaa                                                                                                            <  2:1257588/150‑1 (MQ=255)
                                      ctcGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCaa                                                                                                            <  1:831034/150‑1 (MQ=255)
                                       tcGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAg                                                                                                           >  2:1211692/1‑150 (MQ=255)
                                           gTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAg                                                                                                       <  2:95947/150‑1 (MQ=255)
                                              gCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCAAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTg                                                                                                    >  1:548653/1‑150 (MQ=255)
                                               ccTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGNTgg                                                                                                   <  2:1194114/150‑1 (MQ=255)
                                                cacaaTACACCCAAGGTCAAAGTGATGGGGACCGCGGTTGACGAAATTCCGCTGGGTACGGCAGAACGCGCCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTggg                                                                                                  >  2:1321047/4‑150 (MQ=255)
                                                   aaTACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAAGGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCAt                                                                                               >  1:98112/1‑150 (MQ=255)
                                                        acCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAAt                                                                                                                                     >  2:1285667/1‑107 (MQ=255)
                                                        acCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAAt                                                                                                                                     <  1:1285667/107‑1 (MQ=255)
                                                        acCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGcc                                                                                          >  2:463606/1‑150 (MQ=255)
                                                        acCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGcc                                                                                          >  1:1307514/1‑150 (MQ=255)
                                                          ccAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCCAGGCAGGTGGGCATGCGCCtg                                                                                        >  1:1271546/1‑150 (MQ=255)
                                                          ccAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGc                                                                                                 >  1:1764946/1‑141 (MQ=255)
                                                          ccAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGc                                                                                                 <  2:1764946/141‑1 (MQ=255)
                                                          ccAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtg                                                                                        <  1:847602/150‑1 (MQ=255)
                                                           cAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAgg                                                                                                          <  2:1651340/131‑1 (MQ=255)
                                                           cAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAgg                                                                                                          >  1:1651340/1‑131 (MQ=255)
                                                           cAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtga                                                                                       >  1:937841/1‑150 (MQ=255)
                                                            aaGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAAGCGCGCCATTACTGgcgcgcctgg                                                                                                                          >  1:1281980/1‑110 (MQ=255)
                                                            aaGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAAGCGCGCCATTACTGgcgcgcctgg                                                                                                                          <  2:1281980/114‑5 (MQ=255)
                                                                gCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtgaaatga                                                                                  <  2:256752/150‑1 (MQ=255)
                                                                  aaaGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtgaaatgaaa                                                                                <  1:1211692/150‑1 (MQ=255)
                                                                  aaaGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtgaaatgaaa                                                                                <  2:937841/150‑1 (MQ=255)
                                                                  aaaGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCtgaaatgaaa                                                                                <  2:1239260/150‑1 (MQ=255)
                                                                          ggTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTg                                                                                                    <  2:953532/122‑1 (MQ=255)
                                                                          ggTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTg                                                                                                    >  1:953532/1‑122 (MQ=255)
                                                                             gACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTg                                                                     <  2:98112/150‑1 (MQ=255)
                                                                                 gTGNTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCAt                                                                 >  2:1191068/1‑150 (MQ=255)
                                                                                  tGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCAtt                                                                >  1:1286766/1‑150 (MQ=255)
                                                                                   ggTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCAttt                                                               <  1:186739/150‑1 (MQ=255)
                                                                                        gCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTaaa                                                                           <  2:744091/133‑1 (MQ=255)
                                                                                        gCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTaaa                                                                           >  1:744091/1‑133 (MQ=255)
                                                                                          gAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATcc                                                        >  2:1684341/1‑150 (MQ=255)
                                                                                             aTTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATc                                                     >  1:1770766/1‑150 (MQ=255)
                                                                                              ttCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCg                                                    <  2:1286766/150‑1 (MQ=255)
                                                                                                  gCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTCAAGCATATGCAGTTGCAATCCATCGGCaa                                                >  2:903721/1‑150 (MQ=255)
                                                                                                       gTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAAtttt                                           >  2:1392522/1‑150 (MQ=255)
                                                                                                          cGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCta                                        <  2:624618/150‑1 (MQ=255)
                                                                                                             cAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCAAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAg                                     <  2:548653/150‑1 (MQ=255)
                                                                                                                 aCGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCa                                 <  2:1232858/150‑1 (MQ=255)
                                                                                                                   gcgTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAgg                               >  2:1737846/1‑150 (MQ=255)
                                                                                                                   gcgTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAgg                               <  1:1191068/150‑1 (MQ=255)
                                                                                                                   gcgTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAgg                               <  1:1089535/150‑1 (MQ=255)
                                                                                                                        gCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTg                          >  2:1728231/1‑150 (MQ=255)
                                                                                                                        gCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTg                          <  1:1234373/150‑1 (MQ=255)
                                                                                                                          ccTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGaa                        >  1:1112840/1‑150 (MQ=255)
                                                                                                                                 aTGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACTACTCAAATGAAATTAAAGCATTGTTATTTGCAATCCATCGGCAAATTTTCTATAGGGCAGGCTCTGAAGAGCCTg                 >  2:1004754/1‑150 (MQ=255)
                                                                                                                                 aTGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTg                 <  1:918725/150‑1 (MQ=255)
                                                                                                                                   gcgcCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAAGGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGc               <  1:1612076/150‑1 (MQ=255)
                                                                                                                                        cccAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCtttt          >  2:1387284/1‑150 (MQ=255)
                                                                                                                                         ccAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTc                                          >  2:1097589/1‑117 (MQ=255)
                                                                                                                                         ccAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTc                                          <  1:1097589/117‑1 (MQ=255)
                                                                                                                                                gCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTtg  <  2:1770766/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
ACCGACATGCGCACCCAGCTAGGCCCGACCTCGGCAGTCTCGCGTCGCCTGAATACACCCAAGGGCAAAGTGATGGTGACCGTGGTTGGCGAAATTCCGCTGGGTACGGCAGAACGCGTCGCCCTGTCGATGCGCCCCCAGGATGCCCAGGCGCGCCAGTAATGGCGCGCTTTTGCGTGCCGAGCCAAGGCAGGTGGGCATACGCCTGAAATGAAATTAAAGCATTGTCATTTGCAATCCATCGGCAAATTTTCTATAGGTCAGGCTTTGAAGAGCCTGGCCTTTTCTGCGTTG  >  NC_002947/1629337‑1629630

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