Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,770,016 T→A 5.0% M255L (ATG→TTG)  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,0160TA 5.0% 540.9 / 11.8 159M255L (ATG→TTG) PP_4218lipase/esterase family protein
Reads supporting (aligned to +/- strand):  ref base T (94/57);  new base A (2/6);  total (96/63)
Fisher's exact test for biased strand distribution p-value = 5.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 159 total aligned reads displayed.

TCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACC  >  NC_002947/4769874‑4770164
                                                                                                                                              |                                                                                                                                                    
tcctcGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCa                                                                                                                                               <  1:1590269/150‑1 (MQ=255)
 cctcGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCTGCAACAGTGCac                                                                                                                                              <  1:2677515/150‑2 (MQ=255)
  ctcGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAgg                                                                                                                                             >  2:3583351/1‑150 (MQ=255)
   tcGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGt                                                                                                                                            <  1:481781/150‑1 (MQ=255)
   tcGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGt                                                                                                                                            >  1:3740748/1‑150 (MQ=255)
      atgatCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCa                                                                                                                                         >  2:1756793/1‑150 (MQ=255)
           cGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGtt                                                                                                                                    <  2:2714986/150‑1 (MQ=255)
           cGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCACCTGCACTgat                                                                                                                                    <  1:1679703/150‑3 (MQ=255)
           cGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCATGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGtt                                                                                                                                    >  2:3659650/1‑150 (MQ=255)
            gAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCACGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTGATg                                                                                                                                   <  2:4142677/150‑1 (MQ=255)
            gAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTg                                                                                                                                   <  2:681469/150‑1 (MQ=255)
            gAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTg                                                                                                                                   <  2:678584/150‑1 (MQ=255)
            gAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTg                                                                                                                                   >  2:711596/1‑150 (MQ=255)
            gAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTg                                                                                                                                   >  1:4439273/1‑150 (MQ=255)
            gAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCACCTGCACTGATg                                                                                                                                   <  2:3628773/150‑1 (MQ=255)
            gAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGACATGAAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTg                                                                                                                                   >  2:2336209/1‑150 (MQ=255)
                ccATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTgg                                                                                                                               >  2:2028723/1‑150 (MQ=255)
                  aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCg                                                                                                                             >  2:3292219/1‑150 (MQ=255)
                  aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCg                                                                                                                             >  1:3923037/1‑150 (MQ=255)
                  aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCg                                                                                                                             <  2:2387815/150‑1 (MQ=255)
                  aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCACGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGATGCTGGCg                                                                                                                             >  2:392996/1‑150 (MQ=255)
                   tACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGa                                                                                                                            >  2:671626/1‑150 (MQ=255)
                    aCGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGAc                                                                                                                           <  1:830677/150‑1 (MQ=255)
                    aCGAATGAAGTCATGGGTCACTCCCGTGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGAc                                                                                                                           <  1:2307129/150‑1 (MQ=255)
                      gCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCACTCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGc                                                                                                                         >  2:4186977/1‑150 (MQ=255)
                      gCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGtt                                                                                                                                    <  1:5034123/139‑1 (MQ=255)
                             gTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAg                                                                                                                  <  2:337747/150‑1 (MQ=255)
                              tCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGc                                                                                                                 >  1:3158383/1‑150 (MQ=255)
                                aTGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCaa                                                                                                               >  1:2066677/1‑150 (MQ=255)
                                aTGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCaa                                                                                                               >  1:2521463/1‑150 (MQ=255)
                                 tGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAg                                                                                                                  >  2:4605642/1‑146 (MQ=255)
                                 tGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAg                                                                                                                  <  1:4605642/146‑1 (MQ=255)
                                   ggTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACgg                                                                                                            >  2:1298323/1‑150 (MQ=255)
                                     tCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTg                                                                                                          >  2:1094503/1‑150 (MQ=255)
                                         tCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGaaaa                                                                                                      >  2:3111616/1‑150 (MQ=255)
                                         tCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGCGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGaaaa                                                                                                      >  2:4707755/1‑150 (MQ=255)
                                           ccGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAc                                                                                                    >  1:54964/1‑150 (MQ=255)
                                           ccGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAc                                                                                                    >  2:234508/1‑150 (MQ=255)
                                             gggATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAcgc                                                                                                  <  1:1079661/150‑1 (MQ=255)
                                              ggATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAcgcg                                                                                                 >  1:4401744/1‑150 (MQ=255)
                                                 tCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGAt                                                                                              >  2:2199626/1‑150 (MQ=255)
                                                 tCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGAt                                                                                              >  1:4915709/1‑150 (MQ=255)
                                                    ccGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGt                                                                                           >  1:902699/1‑150 (MQ=255)
                                                    ccGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGt                                                                                           >  1:902959/1‑150 (MQ=255)
                                                      gCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGa                                                                                         <  1:652533/150‑1 (MQ=255)
                                                       cAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGa                                                                                               <  1:623681/143‑1 (MQ=255)
                                                            tGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCt                                                                                   >  1:4405596/1‑150 (MQ=255)
                                                               aCCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGcc                                                                                <  1:5013438/150‑1 (MQ=255)
                                                               aCCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGcc                                                                                <  1:4721283/150‑1 (MQ=255)
                                                                  gCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGc                                                                             >  1:4369323/1‑150 (MQ=255)
                                                                      ccccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCGTCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCg                                                                         >  1:687849/1‑150 (MQ=255)
                                                                      ccccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCg                                                                         >  1:687863/1‑150 (MQ=255)
                                                                          ggCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGc                                                                     <  2:5299876/150‑1 (MQ=255)
                                                                             ttGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTa                                                                  <  1:3278563/150‑1 (MQ=255)
                                                                             ttGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTa                                                                  >  1:1763435/1‑150 (MQ=255)
                                                                                    aaGTGCCGTGCATCAGCGACCCCCCGGACAGGCAGCGGGTCGCACTCGGCACTCAGCCTCAGTGCAGGTGCACTGTTGCTGGCGCCGCGGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCggg                                                           <  2:2791035/150‑1 (MQ=255)
                                                                                    aaGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                           <  1:825682/150‑1 (MQ=255)
                                                                                         ccGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAg                                                                                                                                              >  2:356931/1‑62 (MQ=255)
                                                                                         ccGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAg                                                                                                                                              <  1:356931/62‑1 (MQ=255)
                                                                                                tAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTcc                                               >  1:3749091/1‑150 (MQ=255)
                                                                                                tAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGGGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTCTTGCTGGTAGAACCACTcc                                               >  2:549153/1‑150 (MQ=255)
                                                                                                  aGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGt                                                                                           <  1:1419321/104‑1 (MQ=255)
                                                                                                  aGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGt                                                                                           >  2:1419321/1‑104 (MQ=255)
                                                                                                  aGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAg                                             >  2:5048939/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                                                                   >  2:3165903/1‑124 (MQ=255)
                                                                                                      acacCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGt                                                                   <  1:3165903/124‑1 (MQ=255)
                                                                                                        acCCTGGTCATGCAGCGGGTCTCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTg                                       <  2:1965703/150‑1 (MQ=255)
                                                                                                        acCCTGGTCATGCAGCGGGTAGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTg                                       <  2:2877169/150‑1 (MQ=255)
                                                                                                          ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcg                                     >  2:3868281/1‑150 (MQ=255)
                                                                                                          ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcg                                     >  1:2450248/1‑150 (MQ=255)
                                                                                                          ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcg                                     >  1:4173421/1‑150 (MQ=255)
                                                                                                           cTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcgc                                    >  1:2814166/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)
                                                                                                                cATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCaa                               <  1:2336209/150‑1 (MQ=255)
                                                                                                                 aTGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCCCCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAg                              >  2:933217/1‑150 (MQ=255)
                                                                                                                  tGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGc                             <  1:1872120/150‑1 (MQ=255)
                                                                                                                  tGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGc                             <  1:3689106/150‑1 (MQ=255)
                                                                                                                  tGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGATCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGc                             <  2:5343982/150‑1 (MQ=255)
                                                                                                                   gCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCa                            >  1:2024419/1‑150 (MQ=255)
                                                                                                                    cAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAg                           >  1:420269/1‑150 (MQ=255)
                                                                                                                    cAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAg                           >  1:234587/1‑150 (MQ=255)
                                                                                                                     aGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAgg                          >  1:4740717/1‑150 (MQ=255)
                                                                                                                          gTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAAcc                     >  2:4849917/1‑150 (MQ=255)
                                                                                                                          gTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAAcc                     <  2:4933160/150‑1 (MQ=255)
                                                                                                                                 tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCa              >  2:867783/1‑150 (MQ=255)
                                                                                                                                 tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCa              >  1:3554182/1‑150 (MQ=255)
                                                                                                                                 tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCa              >  2:2771236/1‑150 (MQ=255)
                                                                                                                                   ggCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATa            <  1:1094503/150‑1 (MQ=255)
                                                                                                                                      aaTCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATcgcg         >  2:5302303/1‑150 (MQ=255)
                                                                                                                                        tCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCt       >  1:692100/1‑150 (MQ=255)
                                                                                                                                          aGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGc     >  1:1947053/1‑150 (MQ=255)
                                                                                                                                          aGCAACAGTGCACCTGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGc     >  1:3371595/1‑150 (MQ=255)
                                                                                                                                           gCAACAGTGCACCTGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAgg                          <  2:4710094/128‑1 (MQ=255)
                                                                                                                                           gCAACAGTGCACCTGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAgg                          >  1:4710094/1‑128 (MQ=255)
                                                                                                                                             aTCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCAcc  >  2:3219418/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TCCTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACC  >  NC_002947/4769874‑4770164

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: