Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,896,142 C→A 8.9% R216L (CGG→CTG)  PP_3439 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,896,1420CA 8.9% 347.9 / 14.8 112R216L (CGG→CTG) PP_3439AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (49/53);  new base A (5/5);  total (54/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

Only 100 of 112 total aligned reads displayed.

AGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATGGGC  >  NC_002947/3895994‑3896284
                                                                                                                                                    |                                                                                                                                              
aGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCg                                                                                                                                               <  1:697682/150‑1 (MQ=255)
 ggAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAGTTGCGCCAgg                                                                                                                                              <  1:1099561/150‑1 (MQ=255)
   aGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGc                                                                                                                                            <  2:1940536/150‑1 (MQ=255)
   aGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAGTTGCGCCAGGGc                                                                                                                                            >  1:2048684/1‑150 (MQ=255)
    gcgcCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGcc                                                                                                                                           <  2:672707/150‑1 (MQ=255)
       cccAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGt                                                                                                                                        >  2:97678/1‑150 (MQ=255)
       cccAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGAGGCCCTGGCGCAACTGCGCCCGGGCCAGt                                                                                                                                        >  1:1891084/1‑150 (MQ=255)
        ccAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGtt                                                                                                                                       <  2:2244133/150‑1 (MQ=255)
        ccAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGtt                                                                                                                                       >  2:1112275/1‑150 (MQ=255)
        ccAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAGTTGCGCCAGGGccgctc                                                                                                                                       <  2:2048684/150‑5 (MQ=255)
          aGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGGGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGGGCCCCTGGCGCAACTGCCCCCCCGCCAGTGGc                                                                                                                                     >  2:140484/1‑150 (MQ=255)
            gCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCAc                                                                                                                                   <  1:964069/150‑1 (MQ=255)
            gCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATCGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCAc                                                                                                                                   <  1:140484/150‑1 (MQ=255)
              tAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCAcgc                                                                                                                                 >  2:312055/1‑150 (MQ=255)
                gCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACgcgc                                                                                                                               <  1:1219642/150‑1 (MQ=255)
                  gACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCa                                                                                                                                          >  1:1254671/1‑137 (MQ=255)
                  gACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCa                                                                                                                                          <  2:1254671/137‑1 (MQ=255)
                    cGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGt                                                                                                                           <  2:525409/150‑1 (MQ=255)
                       cATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTa                                                                                                                        >  2:1175325/1‑150 (MQ=255)
                        aTTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAg                                                                                                                       >  2:927783/1‑150 (MQ=255)
                        aTTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAg                                                                                                                       >  2:228139/1‑150 (MQ=255)
                         tttCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGc                                                                                                                      >  2:2287144/1‑150 (MQ=255)
                          ttCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAGTTGCGCCAggg                                                                                                                                             >  2:1491942/1‑126 (MQ=255)
                              gTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGcc                                                                                                                                           <  2:226657/124‑1 (MQ=255)
                              gTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGcc                                                                                                                                           >  1:226657/1‑124 (MQ=255)
                               tAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCCATTAACCCCACCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAg                                                                                                                >  2:65417/1‑150 (MQ=255)
                                aaTGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGc                                                                                                               <  1:228139/150‑1 (MQ=255)
                                 aTGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGcc                                                                                                              >  1:1145038/1‑150 (MQ=255)
                                     cccGGTGTCGCCGCCAGGTGCTGTTTGAATGCCTGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAGTTGCGCCAGGGc                                                                                                                                            <  1:741474/116‑1 (MQ=255)
                                     cccGGTGTCGCCGCCAGGTGCTGTTTGAATGCCTGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAGTTGCGCCAGGGc                                                                                                                                            >  2:741474/1‑116 (MQ=255)
                                     cccGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgc                                                                                                          >  1:1310193/1‑150 (MQ=255)
                                     cccGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgc                                                                                                          <  1:2287144/150‑1 (MQ=255)
                                      ccGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgcg                                                                                                         >  2:737988/1‑150 (MQ=255)
                                              gccgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAg                                                                                                                                         >  1:2119792/1‑110 (MQ=255)
                                              gccgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAg                                                                                                                                         <  2:2119792/110‑1 (MQ=255)
                                                cgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAg                                                                                                                       <  2:630317/126‑1 (MQ=255)
                                                   cAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGaa                                                                                            <  1:1175325/150‑1 (MQ=255)
                                                    aGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAg                                                                                           <  1:1685651/150‑1 (MQ=255)
                                                        gCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGtt                                                                                       <  2:1891084/150‑1 (MQ=255)
                                                         cTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTg                                                                                      >  2:478064/1‑150 (MQ=255)
                                                           gTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgcggc                                                                                                       <  1:1545934/131‑1 (MQ=255)
                                                           gTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgcggc                                                                                                       >  2:1545934/1‑131 (MQ=255)
                                                           gTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGc                                                                                    >  1:1189858/1‑150 (MQ=255)
                                                            tttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCa                                                                                                                     >  2:1239814/1‑116 (MQ=255)
                                                                aaTGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAg                                                                               <  1:927783/150‑1 (MQ=255)
                                                                aaTGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAg                                                                               <  2:1572035/150‑1 (MQ=255)
                                                                aaTGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAg                                                                               >  2:229904/1‑150 (MQ=255)
                                                                  tGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCgg                                                                                                                                              >  1:1300768/1‑85 (MQ=255)
                                                                  tGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCgg                                                                                                                                              <  2:1300768/85‑1 (MQ=255)
                                                                  tGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGc                                                                             >  2:162473/1‑150 (MQ=255)
                                                                   gCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCa                                                                            >  2:1687053/1‑150 (MQ=255)
                                                                    cccGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCAc                                                                           >  2:481439/1‑150 (MQ=255)
                                                                    cccGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCAc                                                                           <  2:1189858/150‑1 (MQ=255)
                                                                     ccGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgcggc                                                                                                       >  2:1919733/1‑121 (MQ=255)
                                                                     ccGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgcggc                                                                                                       <  1:1919733/121‑1 (MQ=255)
                                                                       gCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGaa                                                                        >  1:729487/1‑150 (MQ=255)
                                                                        cTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAAt                                                                       <  1:97678/150‑1 (MQ=255)
                                                                         tGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATc                                                                      <  1:959335/150‑1 (MQ=255)
                                                                         tGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATc                                                                      <  1:162473/150‑1 (MQ=255)
                                                                          gCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCa                                                                     <  2:1741060/150‑1 (MQ=255)
                                                                           cAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGaa                                                                                            >  1:119629/1‑126 (MQ=255)
                                                                           cAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGaa                                                                                            <  2:119629/126‑1 (MQ=255)
                                                                              gtgtGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTa                                                                 <  1:1687053/150‑1 (MQ=255)
                                                                                  gCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCCGGAGTTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAc                                                             <  1:737988/150‑1 (MQ=255)
                                                                                   ccTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTgg                                                                                     >  1:1958374/1‑125 (MQ=255)
                                                                                   ccTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTgg                                                                                     <  2:1958374/125‑1 (MQ=255)
                                                                                   ccTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAcc                                                            >  2:4906/1‑150 (MQ=255)
                                                                                   ccTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAcc                                                            >  2:557395/1‑150 (MQ=255)
                                                                                   ccTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAcc                                                            >  2:2111258/1‑150 (MQ=255)
                                                                                     tGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAgg                                                                   <  1:1323475/141‑1 (MQ=255)
                                                                                     tGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAgg                                                                   >  2:1323475/1‑141 (MQ=255)
                                                                                     tGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCga                                                          <  2:226897/150‑1 (MQ=255)
                                                                                     tGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGAAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCga                                                          >  2:1172563/1‑150 (MQ=255)
                                                                                        ccAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGNGTCAGGCCGAAGCGTTGGCCAAAGGCACGAATCAGGTAGGACCGAGAc                                                       <  1:1644262/149‑1 (MQ=255)
                                                                                         cAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACa                                                      <  1:1112275/150‑1 (MQ=255)
                                                                                          aGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAgg                                                                                                                   >  2:165485/1‑88 (MQ=255)
                                                                                          aGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAgg                                                                                                                   <  1:165485/88‑1 (MQ=255)
                                                                                            cGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCCGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAgg                                                                   >  1:313878/1‑134 (MQ=255)
                                                                                            cGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCCGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAgg                                                                   <  2:313878/134‑1 (MQ=255)
                                                                                            cGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGc                                                   <  1:243576/150‑1 (MQ=255)
                                                                                              aaCCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCg                                                 >  1:1603456/1‑150 (MQ=255)
                                                                                                cccGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCAt                                                                                                            >  2:1305465/1‑89 (MQ=255)
                                                                                                cccGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCAt                                                                                                            <  1:1305465/89‑1 (MQ=255)
                                                                                                       tcctGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGTACCGAGACAGGCCGCAAGCCCCg                                        <  1:1172563/150‑1 (MQ=255)
                                                                                                         ctGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCa                                      <  1:2111258/150‑1 (MQ=255)
                                                                                                            tgCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGTACCGAGACAGGCCGCAAGCCCCGCAGAt                                   <  1:1920882/150‑1 (MQ=255)
                                                                                                            tgCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGAt                                   >  2:1357012/1‑150 (MQ=255)
                                                                                                                  ttCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTc                             >  1:2054596/1‑150 (MQ=255)
                                                                                                                          ttGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGa                     <  1:478064/150‑1 (MQ=255)
                                                                                                                               cGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCggg                >  2:348371/1‑150 (MQ=255)
                                                                                                                                   ccgggCGCAGTTGCGCCAGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCAc            >  2:903817/4‑150 (MQ=255)
                                                                                                                                   ccTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGAGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGcccc                                         >  2:2235835/1‑121 (MQ=255)
                                                                                                                                   ccTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGAGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGcccc                                         <  1:2235835/121‑1 (MQ=255)
                                                                                                                                    cTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGAt                                   <  1:1625666/126‑1 (MQ=255)
                                                                                                                                        cgcAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGa       <  1:481439/150‑1 (MQ=255)
                                                                                                                                         gcAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGAt      <  1:1357012/150‑1 (MQ=255)
                                                                                                                                           aaCTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAgg                                                    >  1:1394005/1‑102 (MQ=255)
                                                                                                                                           aaCTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAgg                                                    <  2:1394005/102‑1 (MQ=255)
                                                                                                                                            gttGCGCCAGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATggg   >  1:1010522/3‑150 (MQ=255)
                                                                                                                                             cTGCGCCCGGGCCAGTGGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATGGgc  <  1:2036817/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
AGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATGGGC  >  NC_002947/3895994‑3896284

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: