Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,301,538 C→A 6.1% intergenic (‑139/‑192) sbcD ← / → PP_5497 exonuclease subunit SbcD/rhs repeat‑associated core domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,301,5380CA 6.1% 265.9 / 3.3 81intergenic (‑139/‑192)sbcD/PP_5497exonuclease subunit SbcD/rhs repeat‑associated core domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (38/38);  new base A (2/3);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGTTGATTGCGGCGGAAAAATTGT  >  NC_002947/2301389‑2301686
                                                                                                                                                     |                                                                                                                                                    
aaCAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGAccc                                                                                                                                                      >  1:1169188/1‑150 (MQ=255)
aaCAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGAccc                                                                                                                                                      >  2:812212/1‑150 (MQ=255)
    gACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTG                                                                                                                                                  <  1:448142‑M1/150‑1 (MQ=255)
     tcGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGCGCAACACGATCGCGGCCCGGTCGCCGCAGATATAAGTTCTCATCACTTACTGATTTAGCAGGGACCCTGTgg                                                                                                                                                 <  1:528135/149‑1 (MQ=255)
     aCGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGG                                                                                                                                                 <  2:938281‑M1/150‑1 (MQ=255)
     aCGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGG                                                                                                                                                 <  2:848573‑M1/150‑1 (MQ=255)
       gCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCGAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGA                                                                                                                                               >  2:528135‑M1/1‑150 (MQ=255)
       gCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGA                                                                                                                                               >  2:904956‑M1/1‑150 (MQ=255)
         aTGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGC                                                                                                                                             >  1:263460‑M1/1‑150 (MQ=255)
         aTGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGC                                                                                                                                             >  2:443945‑M1/1‑150 (MQ=255)
                    gTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCcg                                                                                                                                  >  2:150972/1‑150 (MQ=255)
                     ttGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACGCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGCTTTACCAGGGGGTCTGTGGGAGCGGGCGggccccc                                                                                                                                 >  2:419978‑M1/1‑143 (MQ=255)
                        gtgtCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCgaa                                                                                                                              >  2:1609170/1‑150 (MQ=255)
                          gtCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGccc                                                                                                                                   >  2:317504/1‑143 (MQ=255)
                          gtCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGccc                                                                                                                                   <  1:317504/143‑1 (MQ=255)
                           tCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCgaagaa                                                                                                                           <  1:1480020/150‑1 (MQ=255)
                           tCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCgaagaa                                                                                                                           >  1:764201/1‑150 (MQ=255)
                               gTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCa                                                                                                                       >  1:1149460/1‑150 (MQ=255)
                                tCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGcgcccgcgaacaccggcg                                                                                                                      >  2:264813‑M1/1‑132 (MQ=255)
                                    gCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCgaagaa                                                                                                                           >  2:1501125/1‑141 (MQ=255)
                                    gCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCgaagaa                                                                                                                           <  1:1501125/141‑1 (MQ=255)
                                    gCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAAcgcg                                                                                                                  <  1:23050/150‑1 (MQ=255)
                                      aaGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGt                                                                                                                >  1:797229/1‑150 (MQ=255)
                                          ggggAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCAt                                                                                                            <  2:1733382/150‑1 (MQ=255)
                                          ggggAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCAt                                                                                                            <  1:1609170/150‑1 (MQ=255)
                                           gggAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATg                                                                                                           <  2:1519658/150‑1 (MQ=255)
                                           gggAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATg                                                                                                           <  1:961746/150‑1 (MQ=255)
                                           gggAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATg                                                                                                           <  2:764201/150‑1 (MQ=255)
                                                  ggCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCGATGAAGGATCCAACGCGGTGCATGGCACCgg                                                                                                    <  1:1322583/150‑1 (MQ=255)
                                                  ggCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATc                                                                                                                         >  2:577976/1‑129 (MQ=255)
                                                  ggCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATc                                                                                                                         <  1:577976/129‑1 (MQ=255)
                                                   gCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGc                                                                                                   >  1:207121/1‑150 (MQ=255)
                                                          ggTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGcgcccgcgaacaccggcgcagccggtgccatgcaccgcgttgga                                                                                            <  1:264813‑M1/150‑45 (MQ=255)
                                                           gTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCgg                                                                                           >  1:95095/1‑150 (MQ=255)
                                                           gTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCATGTGGGAGCGG                                                                                                                                           <  2:263018‑M1/102‑1 (MQ=255)
                                                           gTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCATGTGGGAGCGG                                                                                                                                           >  1:263018‑M1/1‑102 (MQ=255)
                                                            tAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGgt                                                                                          <  2:263460/150‑1 (MQ=255)
                                                            tAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGgt                                                                                          <  1:812212/150‑1 (MQ=255)
                                                              aCCCAAGCGTGGGGCCGCTGTGTGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGCGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGgtgt                                                                                        >  1:907235/1‑150 (MQ=255)
                                                               cccAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGtt                                                                                       <  2:1130115/150‑1 (MQ=255)
                                                                  aaGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTcgc                                                                                    >  1:1168064/1‑150 (MQ=255)
                                                                   aGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTcgcg                                                                                   >  2:765921/1‑150 (MQ=255)
                                                                    gCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCgg                                                                                  <  1:544735/150‑1 (MQ=255)
                                                                    gCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCgg                                                                                  <  2:1169188/150‑1 (MQ=255)
                                                                        ggggCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCATGTGGGAGCGGGCGcgcccgcgaacaccggcgcagccggtgccatgcaccgcgttggattcttcgcgggcac                                                                              <  1:1898470‑M1/150‑59 (MQ=255)
                                                                           gCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCAGCACTTACGGATTTAGCAGGCACCCTGTTGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGcggc                                                                           <  1:1694266‑M2/150‑5 (MQ=255)
                                                                              gCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGcgcccgc                                                                        <  1:1454281‑M2/150‑8 (MQ=255)
                                                                                tgtgCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGcgcccgcgaacaccggcgcagccggtgccatgcaccgcgttggattcttcgcgggcacgcccgctc                                                                      <  2:1931715‑M1/150‑67 (MQ=255)
                                                                                   gCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGcgcccgctccca                                                                   >  1:954133‑M2/1‑138 (MQ=255)
                                                                                       cccGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGcgcccgcgaacaccggcgcagccggtgccatgcaccgcgttggattcttcgcgggcacgcccgc                                                                        >  1:188461‑M1/1‑77 (MQ=255)
                                                                                       cccGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGcgcccgcgaacaccggcgcagccggtgccatgcaccgcgttggattcttcgcgggcacgcccgc                                                                        <  2:188461‑M1/141‑65 (MQ=255)
                                                                                        ccGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCa                                                                                                             >  1:1820427/1‑103 (MQ=255)
                                                                                        ccGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCa                                                                                                             <  2:1820427/103‑1 (MQ=255)
                                                                                        ccGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATgg                                                              >  2:494155/1‑150 (MQ=255)
                                                                                         cGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGa                                                             <  1:150972/150‑1 (MQ=255)
                                                                                          gATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGAc                                                            <  2:657330/150‑1 (MQ=255)
                                                                                           aTCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGACACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGAcc                                                           >  1:1738933/1‑150 (MQ=255)
                                                                                            tCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCa                                                          >  1:1830630/1‑150 (MQ=255)
                                                                                                gggCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCt                                                      >  2:1617890/1‑150 (MQ=255)
                                                                                                gggCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCt                                                      >  2:399777/1‑150 (MQ=255)
                                                                                                       ccAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGtgcccgcgaacaccggcgcagccggtgccatgcac                                                                                                     <  2:671446‑M1/96‑36 (MQ=255)
                                                                                                       ccAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGtgcccgcgaacaccggcgcagccggtgccatgcac                                                                                                     >  1:671446‑M1/1‑61 (MQ=255)
                                                                                                                aTAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGt                                      >  1:984935/1‑150 (MQ=255)
                                                                                                                aTAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGt                                      >  1:376544/1‑150 (MQ=255)
                                                                                                                aTAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGt                                      >  1:668936/1‑150 (MQ=255)
                                                                                                                aTAGAACTGCGAATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGt                                      <  1:399777/150‑1 (MQ=255)
                                                                                                                 tAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTc                                     <  1:23174/150‑1 (MQ=255)
                                                                                                                     aCTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGtt                                 <  1:1243345/150‑1 (MQ=255)
                                                                                                                      cTGCGCATCACTTACTGATTTAGCAGGGTCCATGTGGGAGCGGGCGcgcccgcgaacaccggcgcagccggtgccatgcaccgcgttggattcttcgcgg                                                                                  >  2:507163‑M1/1‑46 (MQ=255)
                                                                                                                      cTGCGCATCACTTACTGATTTAGCAGGGTCCATGTGGGAGCGGGCGcgcccgcgaacaccggcgcagccggtgccatgcaccgcgttggattcttcgcgg                                                                                  <  1:507163‑M1/100‑55 (MQ=255)
                                                                                                                            aTCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGc                          <  1:419978/150‑1 (MQ=255)
                                                                                                                                  tACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGTTGa                    <  1:494155/150‑1 (MQ=255)
                                                                                                                                     tGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGTTGATTg                 >  1:590046/1‑150 (MQ=255)
                                                                                                                                     tGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGTTGATTg                 >  1:1338643/1‑150 (MQ=255)
                                                                                                                                       aTTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTcgcg                                                                                   >  2:1795716/1‑82 (MQ=255)
                                                                                                                                       aTTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTcgcg                                                                                   <  1:1795716/82‑1 (MQ=255)
                                                                                                                                            gCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTg                                                <  1:516680/112‑1 (MQ=255)
                                                                                                                                            gCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTg                                                >  2:516680/1‑112 (MQ=255)
                                                                                                                                            gCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGTTGATTGCGGCGGa          >  1:1449685/1‑150 (MQ=255)
                                                                                                                                              aGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGt                       <  2:1888222/135‑1 (MQ=255)
                                                                                                                                              aGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGt                       >  1:1888222/1‑135 (MQ=255)
                                                                                                                                                ggACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGTTGATTGCGGCGGaaaaa      >  1:1499056/1‑150 (MQ=255)
                                                                                                                                                  aCCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGTTGATTGCGGCGGAAAAAtt    <  2:1741200/150‑1 (MQ=255)
                                                                                                                                                    ccTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTc                                     <  2:143877/115‑1 (MQ=255)
                                                                                                                                                    ccTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTc                                     >  1:143877/1‑115 (MQ=255)
                                                                                                                                                    ccTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGTTGATTgcgg              >  1:1071541/1‑138 (MQ=255)
                                                                                                                                                    ccTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGTTGATTgcgg              <  2:1071541/138‑1 (MQ=255)
                                                                                                                                                    ccTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGTTGATTGCGGCGGAAAAATtgt  <  1:443945/150‑1 (MQ=255)
                                                                                                                                                    ccTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGTTGATTGCGGCGGAAAAATtgt  >  1:207269/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACAGACGCATGGATGACCTGTTGGTGTCGAGTCGGGCAAGCGGGGAAAGGGCGGTATGGTAACCCAAGCGTGGGGCCGCTGTGCGGCCCGATCGCGGGCAAGCCAGCTCTCATAGAACTGCGCATCACTTACTGATTTAGCAGGGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAATCCAACGCGGTGCATGGCACCGGCTGCGCCGGTGTTCGCGGGCGCGCCCGCTCCCACATGGACCACGCTTGCCTGGCATTTTTGTCAGTTGTCTCGCCGTTGATTGCGGCGGAAAAATTGT  >  NC_002947/2301389‑2301686

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: