Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,367,931 G→A 46.3% P199S (CCC→TCC)  ftsH ← integral membrane ATP‑dependent zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,367,9310GA46.3% 44.4 / 284.9 203P199S (CCC→TCC) ftsHintegral membrane ATP‑dependent zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (64/45);  new base A (48/46);  total (112/91)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 205 total aligned reads displayed.

GTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTCAGCCAGCGTGGTCTTGACCTGGTC  >  NC_002947/5367782‑5368078
                                                                                                                                                     |                                                                                                                                                   
  gCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATg                                                                                                                                                   <  1:4191208/150‑1 (MQ=255)
   cttcttGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTgg                                                                                                                                                  <  2:602571/150‑1 (MQ=255)
   cttcttGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATgg                                                                                                                                                  <  2:845604/150‑1 (MQ=255)
     tcttGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCAATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGc                                                                                                                                                <  2:2427771/150‑1 (MQ=255)
        tGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGa                                                                                                                                             <  1:1721568/150‑1 (MQ=255)
               cTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAg                                                                                                                                      >  1:142484/1‑150 (MQ=255)
               cTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAg                                                                                                                                      >  1:3340605/1‑150 (MQ=255)
                tCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGc                                                                                                                                     >  1:4275593/1‑150 (MQ=255)
                tCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGc                                                                                                                                     >  1:1140300/1‑150 (MQ=255)
                   aaCATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGcaca                                                                                                                                  >  2:262078/1‑150 (MQ=255)
                    aCATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCacac                                                                                                                                 <  1:3366277/150‑1 (MQ=255)
                     cATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACAcc                                                                                                                                <  1:371424/150‑1 (MQ=255)
                       tGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACcgc                                                                                                                              >  1:2043253/1‑150 (MQ=255)
                         tCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACcgcgc                                                                                                                            >  2:1749443/1‑150 (MQ=255)
                         tCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACcgcgc                                                                                                                            <  2:4172757/150‑1 (MQ=255)
                             gcACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGa                                                                                                                        <  1:4444666/150‑1 (MQ=255)
                                cacGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATAc                                                                                                                     <  1:4521074/150‑1 (MQ=255)
                                 acGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACg                                                                                                                    >  1:3874333/1‑150 (MQ=255)
                                  cGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACgg                                                                                                                   <  2:2957435/150‑1 (MQ=255)
                                   ggCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCATCACACCGCGCGGGATACGGc                                                                                                                  >  2:1022750/1‑150 (MQ=255)
                                    gCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGcc                                                                                                                 <  2:2735490/150‑1 (MQ=255)
                                      tGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCAc                                                                                                               >  2:3829602/1‑150 (MQ=255)
                                        gCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCAccc                                                                                                             >  1:182496/1‑150 (MQ=255)
                                         cACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCa                                                                                                            <  2:3738463/150‑1 (MQ=255)
                                          aCCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAg                                                                                                           >  2:232127/1‑150 (MQ=255)
                                               cacCCACGAACATGTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCt                                                                                                      <  1:2997797/150‑1 (MQ=255)
                                                 cccACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTgg                                                                                                    <  1:697387/150‑1 (MQ=255)
                                                 cccACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTgg                                                                                                    <  1:2172967/150‑1 (MQ=255)
                                                 cccACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTgg                                                                                                    <  1:1480819/150‑1 (MQ=255)
                                                  ccACGAACATTTCCACGAAGTCCGAACCGGAACTGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGa                                                                                                   >  2:3389661/1‑150 (MQ=255)
                                                  ccACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGa                                                                                                   >  2:4241861/1‑150 (MQ=255)
                                                  ccACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGa                                                                                                   >  1:3051197/1‑150 (MQ=255)
                                                   cACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGaa                                                                                                  >  2:490481/1‑150 (MQ=255)
                                                   cACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGaa                                                                                                  <  1:1110364/150‑1 (MQ=255)
                                                    aCGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAAc                                                                                                 >  2:2766320/1‑150 (MQ=255)
                                                    aCGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAAc                                                                                                 >  2:2516705/1‑150 (MQ=255)
                                                      gAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACtt                                                                                               <  2:361593/150‑1 (MQ=255)
                                                       aaCATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTg                                                                                              <  2:3503160/150‑1 (MQ=255)
                                                         cATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGcc                                                                                            <  2:3938659/150‑1 (MQ=255)
                                                         cATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGcc                                                                                            <  2:4257853/150‑1 (MQ=255)
                                                         cATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGcc                                                                                            <  1:1747273/150‑1 (MQ=255)
                                                          aTTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGccc                                                                                           <  1:967329/150‑1 (MQ=255)
                                                           tttCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCg                                                                                          >  1:2561948/1‑150 (MQ=255)
                                                              ccACGAAGTCCGAACCGGAAATGGGGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGAt                                                                                       >  2:766679/1‑150 (MQ=255)
                                                               cACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATc                                                                                      >  1:2212168/1‑150 (MQ=255)
                                                               cACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATc                                                                                      <  1:3138566/150‑1 (MQ=255)
                                                               cACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGCTc                                                                                      >  2:3283754/1‑150 (MQ=255)
                                                                aCGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATcg                                                                                     >  2:2468004/1‑150 (MQ=255)
                                                                aCGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATcg                                                                                     >  1:260393/1‑150 (MQ=255)
                                                                aCGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATcg                                                                                     >  2:3248037/1‑150 (MQ=255)
                                                                 cGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATcgc                                                                                    >  2:2145522/1‑150 (MQ=255)
                                                                  gAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATcgcg                                                                                   >  1:3913437/1‑150 (MQ=255)
                                                                  gAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATcgcg                                                                                   >  1:1779534/1‑150 (MQ=255)
                                                                   aaGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCgcgc                                                                                  >  2:2402614/1‑150 (MQ=255)
                                                                    aGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCa                                                                                 <  2:4320840/150‑1 (MQ=255)
                                                                    aGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCa                                                                                 <  1:2988646/150‑1 (MQ=255)
                                                                    aGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCa                                                                                 <  2:628086/150‑1 (MQ=255)
                                                                      tCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAgg                                                                               >  1:399759/1‑150 (MQ=255)
                                                                       ccGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGa                                                                              >  1:859198/1‑150 (MQ=255)
                                                                       ccGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAg                                                                                                           <  1:800383/121‑1 (MQ=255)
                                                                         gacccGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAAc                                                                            >  2:2002180/4‑150 (MQ=255)
                                                                         gAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGGCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAAc                                                                            >  1:2543220/1‑150 (MQ=255)
                                                                         gAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAAc                                                                            >  1:3637691/1‑150 (MQ=255)
                                                                              ggAAATGGTGAAGAATGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGAc                                                                       <  1:4253793/150‑1 (MQ=255)
                                                                              ggAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGAc                                                                       <  2:3391206/150‑1 (MQ=255)
                                                                               gAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGAcc                                                                      >  1:1258808/1‑150 (MQ=255)
                                                                                      tGAAGAACGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGc                                                               >  2:4108026/1‑150 (MQ=255)
                                                                                       gaagaaCGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGcc                                                              >  1:4146290/1‑150 (MQ=255)
                                                                                        aagaaCGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGACCTCGACCAGCTCGCCa                                                             >  1:748682/1‑150 (MQ=255)
                                                                                         agaaCGGTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCaa                                                            <  1:2226627/150‑1 (MQ=255)
                                                                                               gTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCtt                                                      >  2:1491603/1‑150 (MQ=255)
                                                                                               gTACTTTCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCtt                                                      <  2:3253592/150‑1 (MQ=255)
                                                                                                   tttCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCtt                                                  <  1:1032139/150‑1 (MQ=255)
                                                                                                   tttCGCTTCGCCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCtt                                                  >  1:3191069/1‑150 (MQ=255)
                                                                                                    ttcgcttcgcCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTg                                                 <  1:4112232/150‑1 (MQ=255)
                                                                                                        cttcgcCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCt                                             >  2:514522/1‑150 (MQ=255)
                                                                                                         ttcgcCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCtc                                            <  1:3256404/150‑1 (MQ=255)
                                                                                                         ttcgcCCACGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCtc                                            >  2:3754278/1‑150 (MQ=255)
                                                                                                          tcgcCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCtcg                                           >  1:2486777/1‑150 (MQ=255)
                                                                                                           cgcCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCtcgt                                          <  1:3574516/150‑1 (MQ=255)
                                                                                                               cgcgATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCa                                      >  1:4235848/1‑150 (MQ=255)
                                                                                                                gcgATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAg                                     >  2:2368672/1‑150 (MQ=255)
                                                                                                                       ccTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCaa                              >  2:4437688/1‑150 (MQ=255)
                                                                                                                        cTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAAc                             >  1:171449/1‑150 (MQ=255)
                                                                                                                        cTTGGCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAAc                             <  1:4365047/150‑1 (MQ=255)
                                                                                                                         ttGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACg                            <  1:3220120/150‑1 (MQ=255)
                                                                                                                            gCCAGCAGGGTCTTACCGGTACCGGATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTca                         <  2:4383781/150‑1 (MQ=255)
                                                                                                                             ccAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTcag                        <  2:2275541/150‑1 (MQ=255)
                                                                                                                             ccAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTcag                        >  1:3955619/1‑150 (MQ=255)
                                                                                                                              cagcagGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTcagc                       >  1:396232/1‑150 (MQ=255)
                                                                                                                                 cagGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCCGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTcagccag                    >  2:3417028/1‑150 (MQ=255)
                                                                                                                                    ggTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTCAGCCAGCGt                 >  1:480908/1‑150 (MQ=255)
                                                                                                                                     gTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTCAGCCAGCGTg                <  2:2043253/150‑1 (MQ=255)
                                                                                                                                        ttACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTCAGCCAGCGTGGTc             >  2:4350245/1‑150 (MQ=255)
                                                                                                                                           ccggtaccggATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTCAGCCAGCGTGGTCTTg          >  2:3628182/1‑150 (MQ=255)
                                                                                                                                              gtaccggATGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTCAGCCAGCGTGGTCTTGAcc       >  1:1896086/1‑150 (MQ=255)
                                                                                                                                               taccggGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTCAGCCAGCGTGGTCTTGACCt      <  2:1815556/150‑1 (MQ=255)
                                                                                                                                                accggGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTCAGCCAGCGTGGTCTTGACCTg     <  2:2881305/150‑1 (MQ=255)
                                                                                                                                                  cggGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTCAGCCAGCGTGGTCTTGACCTGGt   <  2:1140300/150‑1 (MQ=255)
                                                                                                                                                   gggTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTCAGCCAGCGTGGTCTTGACCTGGTc  <  2:3051197/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGGGTGGGCCGACCATCAGCACACCGCGCGGGATACGGCCACCCAGGCGCTGGAACTTGCCCGGATCGCGCAGGAACTCGACCAGCTCGCCAACCTCTTCCTTGGCCTCGTCGCAGCCCGCAACGTCAGCCAGCGTGGTCTTGACCTGGTC  >  NC_002947/5367782‑5368078

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