Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,759,091 A→C 7.5% L39R (CTG→CGG)  PP_2412 ← ParA family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,759,0910AC 7.5% 217.8 / 7.4 75L39R (CTG→CGG) PP_2412ParA family protein
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base C (2/4);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GGTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCGGGGTCAATCCTAT  >  NC_002947/2758946‑2759239
                                                                                                                                                 |                                                                                                                                                     
ggTGATTGCCACGCGGTGCTTCTTGCCGGCGACGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCaa                                                                                                                                                   >  2:1917513/1‑150 (MQ=255)
 gTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAg                                                                                                                                                  >  2:2287852/1‑150 (MQ=255)
    aTTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGtt                                                                                                                                               <  2:2241054/150‑1 (MQ=255)
     ttGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTc                                                                                                                                              <  1:952235/150‑1 (MQ=255)
      tGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCTGGTCGACCTCGACCCGCaggccaat                                                                                                                                             <  1:1744097/150‑8 (MQ=255)
        ccACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCTGGTCGACCTCGACCCGCag                                                                                                                                                   <  1:903613/142‑2 (MQ=255)
        ccACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCTGGTCGACCTCGACCCGCag                                                                                                                                                   >  2:903613/1‑141 (MQ=255)
        ccACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGt                                                                                                                                           <  2:185499/150‑1 (MQ=255)
         cACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTa                                                                                                                                          <  2:924516/150‑1 (MQ=255)
         cACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTa                                                                                                                                          <  2:2189277/150‑1 (MQ=255)
         cACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTa                                                                                                                                          >  2:1548519/1‑150 (MQ=255)
          aCGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAg                                                                                                                                         >  2:1787451/1‑150 (MQ=255)
               gTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCAGGGTCGCATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCCTTGGCCTGCGGGTCGCGGTCGACCAGCAAGGTTCGGTAGCCGTc                                                                                                                                    >  2:469114/1‑150 (MQ=255)
                tGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCCGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTcg                                                                                                                                   >  2:1556103/1‑150 (MQ=255)
                  cttcttGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAg                                                                                                                                                  <  2:604346/133‑1 (MQ=255)
                  cttcttGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAg                                                                                                                                                  >  1:604346/1‑133 (MQ=255)
                                 tgTACCCGGCGCCAAGGTCTGGCGGAAGAAGTCGGCGATGCCTGGTTGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGTCCTGCGGGTCGACGTCGACCCGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAgcggc                                                                                                                  <  2:797528/149‑1 (MQ=255)
                                   tAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGCGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAACGGCGa                                                                                                                <  2:1648660/150‑1 (MQ=255)
                                   tAACCGGCGACAAGGTCCGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGa                                                                                                                <  1:735498/150‑1 (MQ=255)
                                     aCCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAgg                                                                                                              <  2:1838641/150‑1 (MQ=255)
                                      ccGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGa                                                                                                                <  2:1820971/147‑1 (MQ=255)
                                      ccGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGa                                                                                                                >  1:1820971/1‑147 (MQ=255)
                                       cGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGtt                                                                                                            >  1:1561215/1‑150 (MQ=255)
                                       cGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGATCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGtt                                                                                                            >  1:1450655/1‑150 (MQ=255)
                                       cGGCGACAAGGTATGGCGGAAGAAGGCGGCGATGCCGGGTGGAATGGTGTCGTGGATCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGCCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGACGtt                                                                                                            <  2:540941/150‑1 (MQ=255)
                                           ggcAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAg                                                                                                        <  2:1395232/148‑1 (MQ=255)
                                           gACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAg                                                                                                        >  1:2355506/1‑150 (MQ=255)
                                              aaGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCg                                                                                                     >  2:1837802/1‑150 (MQ=255)
                                                    tGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTg                                                                                               >  1:12615/1‑150 (MQ=255)
                                                        ggAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGATAGCCCTCGGCGGCGCTGGCAGCGGCGCGGTTGCAGGCGACGCTGGAtt                                                                                           >  2:1167393/1‑150 (MQ=255)
                                                         gaagaaGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGAGAGGTAGAAGGTGGCATTGGCATGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGAGGCGAGGTTGCAGGCGGTGCAGGAttt                                                                                          >  1:222701/1‑150 (MQ=255)
                                                                   gCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCa                                                                                <  1:469114/150‑1 (MQ=255)
                                                                    cGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCAc                                                                               <  1:1917513/150‑1 (MQ=255)
                                                                     gATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCAcc                                                                              >  1:2157719/1‑150 (MQ=255)
                                                                      aTGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTg                                                                                         <  1:1947512/138‑1 (MQ=255)
                                                                      aTGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTg                                                                                         >  2:1947512/1‑138 (MQ=255)
                                                                      aTGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCt                                                                             <  2:1561215/150‑1 (MQ=255)
                                                                          cGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCt                                                                         <  1:278416/150‑1 (MQ=255)
                                                                              tGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGAtt                                                                     <  1:92851/150‑1 (MQ=255)
                                                                              tGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGAtt                                                                     <  1:1556103/150‑1 (MQ=255)
                                                                                 aaTGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGaa                                                                  >  2:1929011/1‑150 (MQ=255)
                                                                                 aaTGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTCGCATTAGCCTGAGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGaa                                                                  >  1:1024578/1‑150 (MQ=255)
                                                                                    ggTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAAc                                                               <  2:18762/150‑1 (MQ=255)
                                                                                      tgtCGTTGACCATCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACca                                                             >  1:2384892/1‑150 (MQ=255)
                                                                                       gtCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACcac                                                            <  1:1929011/150‑1 (MQ=255)
                                                                                        tCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCACCCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACcaca                                                           >  1:55947/1‑150 (MQ=255)
                                                                                                cAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCTGGTCGACCTCGACCCGCag                                                                                                                                                   <  2:244311/54‑2 (MQ=18)
                                                                                                cAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCTGGTCGACCTCGACCCGCag                                                                                                                                                   >  1:244311/1‑53 (MQ=18)
                                                                                                        tGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTa                                                                                                                                          <  2:410/55‑1 (MQ=255)
                                                                                                        tGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTa                                                                                                                                          >  1:410/1‑55 (MQ=255)
                                                                                                        tGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCAcc                                                                              >  1:1060653/1‑115 (MQ=255)
                                                                                                        tGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCAcc                                                                              <  2:1060653/115‑1 (MQ=255)
                                                                                                         gAGGTAGGTGGTGGCATTCGTCTGCGGAGCGTGGTCGACCAGCCAGGTTCGGTCGCCCTCGGCGGCGCTGGCAGCGGCGCGGTTGCAGGCGATGCTGGCTTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTc                                         <  2:55947/150‑1 (MQ=37)
                                                                                                            gTAGAAGGTGGTATTGGCCTGAGGGGCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCATGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATc                                      <  2:222701/150‑1 (MQ=255)
                                                                                                             tAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATcc                                     >  1:483068/1‑150 (MQ=255)
                                                                                                                  ggTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGaa                                <  2:12615/150‑1 (MQ=255)
                                                                                                                       cATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTaaa                           >  2:1295863/1‑150 (MQ=255)
                                                                                                                        aTTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGc                                                                                        >  2:1909599/1‑89 (MQ=255)
                                                                                                                        aTTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGc                                                                                        <  1:1909599/89‑1 (MQ=255)
                                                                                                                        aTTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTaaaa                          >  2:595932/1‑150 (MQ=255)
                                                                                                                            gCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGt                      >  2:184930/1‑150 (MQ=255)
                                                                                                                              cTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCa                    <  2:2157719/150‑1 (MQ=255)
                                                                                                                                gCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCAAGGAAGCAGGCGATACTGGAAGTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAgg                  >  1:901551/1‑150 (MQ=255)
                                                                                                                                   ggtcgaggtcgaCCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGaa                                >  1:778913/1‑133 (MQ=255)
                                                                                                                                   ggtcgaggtcgaCCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGaa                                <  2:778913/133‑1 (MQ=255)
                                                                                                                                     tcgaggtcgaCCAGCAAGGTTCGGTAGCCCTCGGCAGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCggg             <  2:1292346/150‑1 (MQ=255)
                                                                                                                                       gaggtcgaCCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCGGGGt           <  1:2287852/150‑1 (MQ=255)
                                                                                                                                            tgagcAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCGGGGTCAATc      <  1:70690/146‑1 (MQ=255)
                                                                                                                                             gaCCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCGGGGTCAATcc     <  1:1787451/150‑1 (MQ=255)
                                                                                                                                              aCCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGaaaa                                                                >  1:939586/1‑91 (MQ=255)
                                                                                                                                              aCCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGaaaa                                                                <  2:939586/91‑1 (MQ=255)
                                                                                                                                               ccAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCGGGGTCAATCCTa   <  1:1548519/150‑1 (MQ=255)
                                                                                                                                                cAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGTCCACGTTTGCGCTCATCCATGAAGCAAAAATGTCAGGTCGAGGTCAATCCTAt  <  1:674290/150‑1 (MQ=38)
                                                                                                                                                cAGCAAGGCTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCgg              >  1:810864/1‑138 (MQ=255)
                                                                                                                                                cAGCAAGGCTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCgg              <  2:810864/138‑1 (MQ=255)
                                                                                                                                                 aGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCGGGGTCAATCCTa   >  1:2009410/1‑148 (MQ=255)
                                                                                                                                                 aGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCGGGGTCAATCCTa   <  2:2009410/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GGTGATTGCCACGCGGTGCTTCTTGCCGGCGGCGGTAACCGGCGACAAGGTCTGGCGGAAGAAGTCGGCGATGCCGGGTGGAATGGTGTCGTTGACCAGCCCGGTGAGGTAGAAGGTGGCATTGGCCTGCGGGTCGAGGTCGACCAGCAAGGTTCGGTAGCCCTCGGCGGCGCTGGCAGCGGCGAGGTTGCAGGCGATGCTGGATTTGCCGACGCCACCTTTCTGATTGAAAACCACACGACGCATGAGACGCTCATCCATGAAGTAAAAATGTCAGGTCGGGGTCAATCCTAT  >  NC_002947/2758946‑2759239

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