Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,770,024 G→C 6.1% P252R (CCT→CGT) ‡ PP_4218 ← lipase/esterase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,770,0240GC 6.1% 541.9 / 18.7 163P252R (CCT→CGT) ‡PP_4218lipase/esterase family protein
Reads supporting (aligned to +/- strand):  ref base G (93/60);  new base C (2/8);  total (95/68)
Fisher's exact test for biased strand distribution p-value = 1.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

CTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCG  >  NC_002947/4769876‑4770172
                                                                                                                                                    |                                                                                                                                                    
ctcGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAgg                                                                                                                                                     >  2:3583351/1‑150 (MQ=255)
 tcGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGt                                                                                                                                                    >  1:3741143/1‑150 (MQ=255)
 tcGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGt                                                                                                                                                    <  1:481781/150‑1 (MQ=255)
         cGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCACCTGCACTgat                                                                                                                                            <  1:1679703/150‑3 (MQ=255)
         cGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCATGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGtt                                                                                                                                            >  2:3659650/1‑150 (MQ=255)
          gAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTg                                                                                                                                           <  2:678584/150‑1 (MQ=255)
          gAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCACCTGCACTGATg                                                                                                                                           <  2:3628810/150‑1 (MQ=255)
           aaCCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTGATGc                                                                                                                                          <  1:4231155/150‑1 (MQ=255)
           aaCCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTGATGc                                                                                                                                          <  2:1131422/150‑1 (MQ=255)
              ccATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTgg                                                                                                                                       >  2:2028723/1‑150 (MQ=255)
                aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCg                                                                                                                                     >  2:3292219/1‑150 (MQ=255)
                aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAACCAGCAACAGTGCACCTGCACTGATGCTgatt                                                                                                                                     <  2:2787677/150‑4 (MQ=255)
                aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCACGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGATGCTGGCg                                                                                                                                     >  2:392996/1‑150 (MQ=255)
                 tACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGa                                                                                                                                    >  2:671626/1‑150 (MQ=255)
                 tACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGa                                                                                                                                    >  1:1519573/1‑150 (MQ=255)
                    gCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCACTCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGc                                                                                                                                 >  2:4186977/1‑150 (MQ=255)
                    gCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGtt                                                                                                                                            >  2:5034123/1‑139 (MQ=255)
                    gCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGc                                                                                                                                 >  2:4847379/1‑150 (MQ=255)
                     cATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCt                                                                                                                                >  2:5161801/1‑150 (MQ=255)
                            tCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGc                                                                                                                         >  1:3158383/1‑150 (MQ=255)
                              aTGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCaa                                                                                                                       >  1:2521463/1‑150 (MQ=255)
                              aTGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCaa                                                                                                                       >  1:2066677/1‑150 (MQ=255)
                               tGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAg                                                                                                                          <  1:4605642/146‑1 (MQ=255)
                                   tCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTg                                                                                                                  <  2:4647798/150‑1 (MQ=255)
                                       tCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGaaaa                                                                                                              >  2:3111616/1‑150 (MQ=255)
                                       tCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGaaaa                                                                                                              >  1:2317779/1‑150 (MQ=255)
                                       tCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGCGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGaaaa                                                                                                              >  2:4707755/1‑150 (MQ=255)
                                         ccGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAc                                                                                                            >  2:234446/1‑150 (MQ=255)
                                         ccGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAc                                                                                                            >  2:234508/1‑150 (MQ=255)
                                         ccGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAc                                                                                                            >  1:54928/1‑150 (MQ=255)
                                           gggATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAcgc                                                                                                          <  1:1079661/150‑1 (MQ=255)
                                            ggATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGAGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAACAAcgcg                                                                                                         >  1:4401737/1‑150 (MQ=255)
                                               tCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGAt                                                                                                      >  2:2199626/1‑150 (MQ=255)
                                               tCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGAt                                                                                                      >  1:4915709/1‑150 (MQ=255)
                                                cACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATc                                                                                                     <  2:4439273/150‑1 (MQ=255)
                                                  ccGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGt                                                                                                   >  1:902959/1‑150 (MQ=255)
                                                   cGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTg                                                                                                  <  1:3845083/150‑1 (MQ=255)
                                                    gCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGa                                                                                                 <  1:652533/150‑1 (MQ=255)
                                                     cAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGa                                                                                                       >  2:623681/1‑143 (MQ=255)
                                                      aaGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACg                                                                                               <  1:886380/150‑1 (MQ=255)
                                                          tGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCt                                                                                           >  1:4405596/1‑150 (MQ=255)
                                                           gCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTg                                                                                          >  1:600967/1‑150 (MQ=255)
                                                             aCCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGcc                                                                                        <  1:5013438/150‑1 (MQ=255)
                                                                 ccACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCa                                                                                    <  2:1056325/150‑1 (MQ=255)
                                                                    ccccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCg                                                                                 >  2:3152642/1‑150 (MQ=255)
                                                                    ccccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCg                                                                                 >  1:687863/1‑150 (MQ=255)
                                                                     cccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGt                                                                                >  2:1482930/1‑150 (MQ=255)
                                                                                  aaGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                                   <  1:4438228/150‑1 (MQ=255)
                                                                                  aaGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                                   <  1:825682/150‑1 (MQ=255)
                                                                                       ccGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAg                                                                                                                                                      >  2:356931/1‑62 (MQ=255)
                                                                                                aGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGt                                                                                                   >  2:1419321/1‑104 (MQ=255)
                                                                                                 gCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGTCGTATTGCTGGTAGAACCACTCCAgg                                                    <  1:3248483/150‑1 (MQ=255)
                                                                                                 gCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAgg                                                    >  1:4712253/1‑150 (MQ=255)
                                                                                                 gCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAgg                                                    >  1:3472668/1‑150 (MQ=255)
                                                                                                   gACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGt                                                  <  2:4332632/150‑1 (MQ=255)
                                                                                                   gACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGt                                                  <  2:4332452/150‑1 (MQ=255)
                                                                                                    acacCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGt                                                                           <  1:3165903/124‑1 (MQ=255)
                                                                                                    acacCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTc                                                 <  2:701240/150‑1 (MQ=255)
                                                                                                      acCCTGGTCATGCAGCGGGTCTCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTg                                               <  2:1965703/150‑1 (MQ=255)
                                                                                                      acCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTg                                               >  2:240469/1‑150 (MQ=255)
                                                                                                       cccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgc                                              <  1:3583351/150‑1 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcg                                             >  2:412354/1‑150 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcg                                             >  1:4173397/1‑150 (MQ=255)
                                                                                                         cTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcgc                                            >  1:3868059/1‑150 (MQ=255)
                                                                                                           ggTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCtt                                          <  1:2827014/150‑1 (MQ=255)
                                                                                                            gTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTc                                         >  1:1277222/1‑150 (MQ=255)
                                                                                                               aTGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCCCCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAg                                      >  2:933217/1‑150 (MQ=255)
                                                                                                               aTGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAg                                      <  1:1948484/150‑1 (MQ=255)
                                                                                                                tGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGc                                     <  1:3689106/150‑1 (MQ=255)
                                                                                                                tGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGc                                     <  1:1872120/150‑1 (MQ=255)
                                                                                                                  cAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAg                                   <  1:671626/150‑1 (MQ=255)
                                                                                                                  cAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAg                                   >  1:420269/1‑150 (MQ=255)
                                                                                                                  cAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAg                                   >  1:234587/1‑150 (MQ=255)
                                                                                                                  cAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAg                                   <  1:2072114/150‑1 (MQ=255)
                                                                                                                   aGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAgg                                  >  1:4740717/1‑150 (MQ=255)
                                                                                                                        gTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAAcc                             <  2:4933160/150‑1 (MQ=255)
                                                                                                                         tCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCg                            >  2:3359247/1‑150 (MQ=255)
                                                                                                                             aCTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGc                        >  2:5322201/1‑150 (MQ=255)
                                                                                                                               tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTg                         <  2:2795244/147‑1 (MQ=255)
                                                                                                                               tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTg                         >  1:2795244/1‑147 (MQ=255)
                                                                                                                               tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCa                      >  1:3554182/1‑150 (MQ=255)
                                                                                                                               tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCa                      >  2:867783/1‑150 (MQ=255)
                                                                                                                               tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCa                      >  2:2771236/1‑150 (MQ=255)
                                                                                                                                 ggCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATa                    <  1:1094503/150‑1 (MQ=255)
                                                                                                                                 ggCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATa                    <  1:2442599/150‑1 (MQ=255)
                                                                                                                                   cAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAgc                  >  1:1788399/1‑150 (MQ=255)
                                                                                                                                    aaTCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATcgcg                 >  2:5302303/1‑150 (MQ=255)
                                                                                                                                      tCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCt               >  1:692100/1‑150 (MQ=255)
                                                                                                                                        aGCAACAGTGCACCTGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGc             >  1:3371595/1‑150 (MQ=255)
                                                                                                                                         gCAACAGTGCACCTGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAgg                                  >  1:4710094/1‑128 (MQ=255)
                                                                                                                                           aTCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCAcc          >  2:3219418/1‑150 (MQ=255)
                                                                                                                                              aGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGaa       <  2:3749091/150‑1 (MQ=255)
                                                                                                                                              aGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGaa       <  2:2292471/150‑1 (MQ=255)
                                                                                                                                              aGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGaa       <  1:1756793/150‑1 (MQ=255)
                                                                                                                                               gTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAg      <  1:3292219/150‑1 (MQ=255)
                                                                                                                                               gTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAg      <  1:2789591/150‑1 (MQ=255)
                                                                                                                                                  cAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGgcc   <  1:711596/150‑1 (MQ=255)
                                                                                                                                                   aGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGgccg  <  2:691652/150‑1 (MQ=255)
                                                                                                                                                    ggTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAgcg                 >  1:5448867/1‑134 (MQ=255)
                                                                                                                                                    ggTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAgcg                 <  2:5448867/134‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCG  >  NC_002947/4769876‑4770172

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: