Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,526,900 G→C 5.0% R122P (CGC→CCC) ‡ PP_4013 → MutT/nudix family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,526,9000GC 5.0% 670.0 / 18.1 199R122P (CGC→CCC) ‡PP_4013MutT/nudix family protein
Reads supporting (aligned to +/- strand):  ref base G (70/119);  new base C (5/5);  total (75/124)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 209 total aligned reads displayed.

CACCGCCCCCAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGACCAGCCCAG  >  NC_002947/4526751‑4527048
                                                                                                                                                     |                                                                                                                                                    
cacCGCCCCCAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGgcg                                                                                                                                                      >  2:4567901/1‑150 (MQ=255)
 acCGCCCCCAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGgcgc                                                                                                                                                     >  2:3320157/1‑150 (MQ=255)
 acCGCCCCCAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGgcgc                                                                                                                                                     >  1:2217014/1‑150 (MQ=255)
 acCGCCCCCAGCAAAGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGgcgc                                                                                                                                                     >  2:205026/1‑150 (MQ=255)
    gCCCCCAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGc                                                                                                                                                  <  1:4808987/150‑1 (MQ=255)
    gCCCCCAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGc                                                                                                                                                  <  2:5262860/150‑1 (MQ=255)
      ccccAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGa                                                                                                                                                >  2:456325/1‑150 (MQ=255)
       cccAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGCGCCAGCGAg                                                                                                                                               >  1:1122044/1‑150 (MQ=255)
         cAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAgtg                                                                                                                                             <  2:814183/150‑3 (MQ=255)
         cAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCt                                                                                                                                             <  2:4546326/150‑1 (MQ=255)
         cAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGCGCCAGCGAgtg                                                                                                                                             <  2:1122044/150‑3 (MQ=255)
          aGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAgtgg                                                                                                                                            <  2:3867770/150‑4 (MQ=255)
          aGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTg                                                                                                                                            >  1:5006733/1‑150 (MQ=255)
          aGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGCGCCAGCGAgtgg                                                                                                                                            <  1:1042578/150‑4 (MQ=255)
                ggCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCcg                                                                                                                                      <  1:3051803/150‑1 (MQ=255)
                  cGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCcgcg                                                                                                                                    >  2:572492/1‑150 (MQ=255)
                  cGTGACCTACCAGCGCGTCTGCTGCGCCGCCCGCCCCGTGCGCCATCACGCCGCCCGGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCcgcg                                                                                                                                    <  2:4410307/150‑1 (MQ=255)
                        cTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCt                                                                                                                              <  1:695105/150‑1 (MQ=255)
                         tACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTg                                                                                                                             <  1:4343476/150‑1 (MQ=255)
                          aCCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTgg                                                                                                                            <  2:3660139/150‑1 (MQ=255)
                           ccAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTgg                                                                                                                            >  2:1404281/1‑149 (MQ=255)
                           ccAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTgg                                                                                                                            <  1:1404281/149‑1 (MQ=255)
                           ccAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGa                                                                                                                           <  1:2576959/150‑1 (MQ=255)
                            cAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGAt                                                                                                                          >  2:901040/1‑150 (MQ=255)
                            cAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGAt                                                                                                                          <  1:1503457/150‑1 (MQ=255)
                             aGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGc                                                                                                                                                  >  2:5323692/1‑125 (MQ=255)
                             aGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGc                                                                                                                                                  <  1:5323692/125‑1 (MQ=255)
                               cgcgTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGAc                                                                                                                       <  2:5209385/150‑1 (MQ=255)
                               cgcgTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCAGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGTCCCGCGACTAACTGCTTGCCTACCCCACTCGCTGGCGCAGCGTGCTGGTGCCGCGCTGCCTGGATGAc                                                                                                                       <  1:2342165/150‑1 (MQ=255)
                                 cgTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTa                                                                                                                     >  1:4717415/1‑150 (MQ=255)
                                      gCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCa                                                                                                                >  1:2235108/1‑150 (MQ=255)
                                      gCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCa                                                                                                                >  2:4353733/1‑150 (MQ=255)
                                        ttCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCTCTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCaaa                                                                                                              <  1:1244593/150‑1 (MQ=255)
                                        ttCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCaaa                                                                                                              <  2:5006733/150‑1 (MQ=255)
                                        ttCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCaaa                                                                                                              <  2:3855923/150‑1 (MQ=255)
                                           gccgccCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGc                                                                                                           >  2:3153560/1‑150 (MQ=255)
                                              gcccgcccCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCg                                                                                                        <  2:3345046/150‑1 (MQ=255)
                                                        tGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGcc                                                                                              >  2:4071435/1‑150 (MQ=255)
                                                          cgcCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTc                                                                                            <  1:5369808/150‑1 (MQ=255)
                                                             cATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGTGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGAc                                                                                         >  2:2048269/1‑150 (MQ=255)
                                                             cATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGAc                                                                                         <  1:1654534/150‑1 (MQ=255)
                                                             cATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGAc                                                                                         <  1:3837476/150‑1 (MQ=255)
                                                               tCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCt                                                                                       >  2:1167519/1‑150 (MQ=255)
                                                                       aCCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTAcgcg                                                                               >  1:4439165/1‑150 (MQ=255)
                                                                       aCCTGACACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTCcgcg                                                                               >  2:3176189/1‑150 (MQ=255)
                                                                         cTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGAc                                                                             <  2:2885186/150‑1 (MQ=255)
                                                                            gCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCaaa                                                                                                              >  1:1077699/1‑114 (MQ=255)
                                                                            gCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGa                                                                          <  2:1994499/150‑1 (MQ=255)
                                                                                  gACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGAcgcgtc                                                                    <  2:379034/150‑3 (MQ=255)
                                                                                   aCAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACgcggcg                                                                   >  1:2846195/1‑150 (MQ=255)
                                                                                     aGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGcgcc                                                                 >  1:3906930/1‑150 (MQ=255)
                                                                                       cGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGcgccgc                                                               <  2:196211/150‑1 (MQ=255)
                                                                                       cGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGcgccgc                                                               <  2:5138549/150‑1 (MQ=255)
                                                                                          tatCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCggg                                                            <  1:2473007/150‑1 (MQ=255)
                                                                                              gTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTg                                                        <  1:4906524/150‑1 (MQ=255)
                                                                                               tCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGa                                                       <  1:4479479/150‑1 (MQ=255)
                                                                                               tCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGa                                                       <  1:1324193/150‑1 (MQ=255)
                                                                                                 cgcgCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATa                                                     <  2:410405/150‑1 (MQ=255)
                                                                                                    gcCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGaa                                                  >  1:4936141/1‑150 (MQ=255)
                                                                                                    gcCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGaa                                                  <  2:1508728/150‑1 (MQ=255)
                                                                                                      cGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATc                                                <  1:4031200/150‑1 (MQ=255)
                                                                                                         gtgGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGc                                             <  2:3906930/150‑1 (MQ=255)
                                                                                                         gtgGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGc                                             <  1:3988660/150‑1 (MQ=255)
                                                                                                            gCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGtt                                          >  1:1800908/1‑150 (MQ=255)
                                                                                                              tGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGtttt                                        <  2:635472/150‑1 (MQ=255)
                                                                                                               gACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGttt                                         <  2:5409755/148‑1 (MQ=255)
                                                                                                               gACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGttt                                         >  1:5409755/1‑148 (MQ=255)
                                                                                                                 cccGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTcc                                     <  2:2004306/150‑1 (MQ=255)
                                                                                                                 cccGCGACGAACTGCTGGCCGACCCCACCCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTcc                                     <  1:3481970/150‑1 (MQ=255)
                                                                                                                  ccGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTg                                                                           <  2:184362/111‑1 (MQ=255)
                                                                                                                  ccGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTg                                                                           >  1:184362/1‑111 (MQ=255)
                                                                                                                  ccGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGa                                                       >  1:2024200/1‑131 (MQ=255)
                                                                                                                  ccGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTccc                                    <  1:4298259/150‑1 (MQ=255)
                                                                                                                  ccGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTccc                                    >  1:4636809/1‑150 (MQ=255)
                                                                                                                  ccGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTccc                                    >  2:2917705/1‑150 (MQ=255)
                                                                                                                  ccGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTccc                                    >  2:28510/1‑150 (MQ=255)
                                                                                                                   cgcgACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTcccc                                   <  2:4889928/150‑1 (MQ=255)
                                                                                                                   cgcgACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCCGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTcccc                                   <  2:2154763/150‑1 (MQ=255)
                                                                                                                    gcgACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTccccc                                  <  1:1149726/150‑1 (MQ=255)
                                                                                                                        cgaACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGTCCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGa                              <  1:3153560/150‑1 (MQ=255)
                                                                                                                        cgaACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGa                              <  1:3747640/150‑1 (MQ=255)
                                                                                                                         gaACTGCTGGCCGACCCAACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGAt                             >  1:3393324/1‑150 (MQ=255)
                                                                                                                          aaCTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATa                            >  2:4157717/1‑150 (MQ=255)
                                                                                                                              gctgGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATacaca                        <  1:4165315/150‑1 (MQ=255)
                                                                                                                                   ccGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCg                                                                                                        <  1:2457233/65‑1 (MQ=255)
                                                                                                                                   ccGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCg                                                                                                        >  2:2457233/1‑65 (MQ=255)
                                                                                                                                     gACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTAc                 <  2:5393599/150‑1 (MQ=255)
                                                                                                                                      aCCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTAcc                >  1:4513577/1‑150 (MQ=255)
                                                                                                                                        cccACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCa              >  1:3302244/1‑150 (MQ=255)
                                                                                                                                         cagcTCGCTGCGCCAGCGAGTTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCAt             >  2:2496263/4‑150 (MQ=255)
                                                                                                                                         cagcTCGCTGCGCCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCAt             >  2:2287026/4‑150 (MQ=255)
                                                                                                                                          agcTCGCTGCGCCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACt                                                                            <  2:345632/84‑1 (MQ=255)
                                                                                                                                           aCTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGa           <  1:456325/150‑1 (MQ=255)
                                                                                                                                           aCTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGa           <  2:176515/150‑1 (MQ=255)
                                                                                                                                            cTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGAc          <  1:3554473/150‑1 (MQ=255)
                                                                                                                                              cGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGAcca        <  2:2235014/150‑1 (MQ=255)
                                                                                                                                                cTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGAccagc      <  2:3769543/150‑1 (MQ=255)
                                                                                                                                                 tGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGAccagcc     >  2:4726263/1‑150 (MQ=255)
                                                                                                                                                   gcgcAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGAccagccca   <  2:4214097/150‑1 (MQ=255)
                                                                                                                                                   gcgcAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGAccagccca   <  2:2382200/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CACCGCCCCCAGCAACGGCGTGACCTACCAGCGCGTCTGCTTCGCCGCCCGCCCCGTGCGCCATCACGCCGACCTGGCACTGGACAGCGATATCGTCCGCGCCGTGTGGCTGACCCGCGACGAACTGCTGGCCGACCCCACTCGCTGGCGCAGCGAGCTGGTGCCGCGCTGCCTGGATGACTACCTCAAAGGCCCGCTGCACAGCCTCGACCTGCTACGCGACTGACGCGGCGCCGCGGGTTTGATAGAATCGGCGTTTTTCCCCCTTGATACACACCGGTACCCATGACCAGCCCAG  >  NC_002947/4526751‑4527048

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: