Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,770,025 G→C 6.1% P252A (CCT→GCT) ‡ 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,0250GC 6.1% 542.1 / 18.9 163P252A (CCT→GCT) ‡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 = 9.72e-01

Only 100 of 167 total aligned reads displayed.

CTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCC  >  NC_002947/4769876‑4770174
                                                                                                                                                     |                                                                                                                                                     
ctcGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAgg                                                                                                                                                       >  2:3583351/1‑150 (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)
         cGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTACACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGtt                                                                                                                                              <  2:3243038/150‑1 (MQ=255)
          gAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTg                                                                                                                                             <  2:678584/150‑1 (MQ=255)
          gAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCACCTGCACTGATg                                                                                                                                             <  2:3628810/150‑1 (MQ=255)
          gAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGACATGAAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTg                                                                                                                                             >  2:2336209/1‑150 (MQ=255)
           aaCCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCAACAGTGCACCTGCACTGATGc                                                                                                                                            <  1:4231155/150‑1 (MQ=255)
                aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCg                                                                                                                                       >  1:3923037/1‑150 (MQ=255)
                aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCg                                                                                                                                       <  2:2387815/150‑1 (MQ=255)
                aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAACCAGCAACAGTGCACCTGCACTGATGCTgatt                                                                                                                                       <  2:2787677/150‑4 (MQ=255)
                aTACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCACGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGATGCTGGCg                                                                                                                                       >  2:392996/1‑150 (MQ=255)
                 tACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGa                                                                                                                                      >  1:1519573/1‑150 (MQ=255)
                 tACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGa                                                                                                                                      >  2:671626/1‑150 (MQ=255)
                  aCGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGAc                                                                                                                                     <  1:830677/150‑1 (MQ=255)
                  aCGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGAc                                                                                                                                     <  2:2237161/150‑1 (MQ=255)
                    gCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCACTCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGc                                                                                                                                   >  2:4186977/1‑150 (MQ=255)
                    gCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGtt                                                                                                                                              >  2:5034123/1‑139 (MQ=255)
                    gCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGtt                                                                                                                                              <  1:5034123/139‑1 (MQ=255)
                    gCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGc                                                                                                                                   >  2:4847379/1‑150 (MQ=255)
                     cATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCt                                                                                                                                  >  2:5161801/1‑150 (MQ=255)
                           gTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAg                                                                                                                            <  2:337747/150‑1 (MQ=255)
                            tCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGc                                                                                                                           >  1:3158383/1‑150 (MQ=255)
                            tCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGc                                                                                                                           >  1:691652/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)
                                 ggTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACgg                                                                                                                      >  2:1298323/1‑150 (MQ=255)
                                   tCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTg                                                                                                                    <  2:4647798/150‑1 (MQ=255)
                                   tCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTg                                                                                                                    >  2:1094503/1‑150 (MQ=255)
                                       tCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGaaaa                                                                                                                >  1:2317779/1‑150 (MQ=255)
                                       tCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGCGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGaaaa                                                                                                                >  2:1968426/1‑150 (MQ=255)
                                       tCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGCGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGaaaa                                                                                                                >  2:4707755/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)
                                            ggATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAAcgcg                                                                                                           >  1:4401744/1‑150 (MQ=255)
                                            ggATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGAGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAACAAcgcg                                                                                                           >  1:4401737/1‑150 (MQ=255)
                                               tCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGAt                                                                                                        >  1:4915709/1‑150 (MQ=255)
                                               tCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGAt                                                                                                        >  2:2199626/1‑150 (MQ=255)
                                                cACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATc                                                                                                       <  2:4439273/150‑1 (MQ=255)
                                                  ccGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGt                                                                                                     >  1:902699/1‑150 (MQ=255)
                                                   cGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTg                                                                                                    <  1:3845083/150‑1 (MQ=255)
                                                     cAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGa                                                                                                         >  2:623681/1‑143 (MQ=255)
                                                     cAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGa                                                                                                         <  1:623681/143‑1 (MQ=255)
                                                           gCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTg                                                                                            >  1:600967/1‑150 (MQ=255)
                                                             aCCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGcc                                                                                          <  1:5013438/150‑1 (MQ=255)
                                                                 ccACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCa                                                                                      <  2:1056325/150‑1 (MQ=255)
                                                                    ccccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCg                                                                                   >  1:687863/1‑150 (MQ=255)
                                                                    ccccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCg                                                                                   >  2:3152642/1‑150 (MQ=255)
                                                                     cccGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGt                                                                                  >  2:1482930/1‑150 (MQ=255)
                                                                               tCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGc                                                                        >  1:4709824/1‑150 (MQ=255)
                                                                                  aaGTGCCGTGCATCAGCGACCCCCCGGACAGGCAGCGGGTCGCACTCGGCACTCAGCCTCAGTGCAGGTGCACTGTTGCTGGCGCCGCGGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCggg                                                                     <  2:2791035/150‑1 (MQ=255)
                                                                                  aaGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                                     <  1:825682/150‑1 (MQ=255)
                                                                                  aaGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTgg                                                                     <  1:4438228/150‑1 (MQ=255)
                                                                                         gTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGAGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAAcc                                                              >  1:701251/1‑150 (MQ=255)
                                                                                              tAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTcc                                                         >  1:3749091/1‑150 (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:4332452/150‑1 (MQ=255)
                                                                                                   gACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGt                                                    <  2:4332632/150‑1 (MQ=255)
                                                                                                    acacCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGt                                                                             >  2:3165903/1‑124 (MQ=255)
                                                                                                    acacCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGt                                                                             <  1:3165903/124‑1 (MQ=255)
                                                                                                    acacCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTc                                                   <  2:701240/150‑1 (MQ=255)
                                                                                                      acCCTGGTCATGCAGCGGGTAGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTg                                                 <  2:2877169/150‑1 (MQ=255)
                                                                                                       cccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgc                                                <  1:3583351/150‑1 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcg                                               >  1:4173397/1‑150 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcg                                               >  1:4173421/1‑150 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcg                                               >  2:412354/1‑150 (MQ=255)
                                                                                                        ccTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCATGGCGTATTGCTGGTAGAACCACTCCAGGGTCTgcg                                               >  1:4791339/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:1872120/150‑1 (MQ=255)
                                                                                                                  cAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAg                                     <  1:671626/150‑1 (MQ=255)
                                                                                                                  cAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAg                                     <  1:2072114/150‑1 (MQ=255)
                                                                                                                  cAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAg                                     >  1:234587/1‑150 (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                           >  1:2795244/1‑147 (MQ=255)
                                                                                                                               tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCa                        >  1:3554182/1‑150 (MQ=255)
                                                                                                                               tCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCa                        >  2:867783/1‑150 (MQ=255)
                                                                                                                                 ggCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATa                      <  1:2442599/150‑1 (MQ=255)
                                                                                                                                 ggCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATa                      <  1:1094503/150‑1 (MQ=255)
                                                                                                                                   cAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAgc                    >  1:1788399/1‑150 (MQ=255)
                                                                                                                                    aaTCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATcgcg                   >  2:5302303/1‑150 (MQ=255)
                                                                                                                                        aGCAACAGTGCACCTGCACTGATGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGc               >  1:3371595/1‑150 (MQ=255)
                                                                                                                                              aGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGaa         <  2:3749091/150‑1 (MQ=255)
                                                                                                                                              aGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGaa         <  1:1756793/150‑1 (MQ=255)
                                                                                                                                               gTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAg        <  1:2789591/150‑1 (MQ=255)
                                                                                                                                               gTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAg        <  1:3292219/150‑1 (MQ=255)
                                                                                                                                                   aGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGgccg    <  2:691652/150‑1 (MQ=255)
                                                                                                                                                    ggTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAgcg                   <  2:5448867/134‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCGATGATCGAACCCATACGCATGAAGTCATGGGTCACTCCCGGGATCACCGCAAGCTGCACCGCCACCCCGGCTTGCTCCAAGTGCCGTGCATAAGCGACACCCTGGTCATGCAGCGGGTCGCACTCGGCAATCAGCATCAGTGCAGGTGCACTGTTGCTGGCGACGCTGCCCAGCAACGGTGAAAAACGCGGATCGTGACGGTCTGCCGGCACCGTGGCGTATTGCTGGTAGAACCACTCCAGGGTCTGCGCTTCAAGCAGGTAACCGCTGCCATAGCGCTGCACCGAAGGCCGCC  >  NC_002947/4769876‑4770174

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: