Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,646,791 A→C 5.6% Q21P (CAG→CCG)  PP_2312 → lytic transglycosylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,646,7910AC 5.6% 1413.6 / 48.2 415Q21P (CAG→CCG) PP_2312lytic transglycosylase family protein
Reads supporting (aligned to +/- strand):  ref base A (172/219);  new base C (17/7);  total (189/226)
Fisher's exact test for biased strand distribution p-value = 1.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 420 total aligned reads displayed.

CATGCTGTTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGAGT  >  NC_002947/2646642‑2646940
                                                                                                                                                     |                                                                                                                                                     
cATGCTGTTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCa                                                                                                                                                       <  1:428990/150‑1 (MQ=255)
   gCTGTTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGgc                                                                                                                                                    >  2:2283091/1‑150 (MQ=255)
   gCTGTTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGgc                                                                                                                                                    <  2:3706327/150‑1 (MQ=255)
   gCTGTTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGgc                                                                                                                                                    <  2:7521420/150‑1 (MQ=255)
         gCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTc                                                                                                                                              <  2:361662/150‑1 (MQ=255)
          ccGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCg                                                                                                                                             >  2:8223475/1‑150 (MQ=255)
             tGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCAAGGATGCCGCTCACCGACCACGGACCCCTTTTCGCATGCAGCGATGCAGGCTCATTTTTATGCTGACCCGGGGGGTGGCCTTTGCCCGTCCCCCCCACGCGCGCTCGCCg                                                                                                                                          >  1:2443206/1‑150 (MQ=255)
                 cTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCCGGCGAGCGCGCCTGGGc                                                                                                                                      <  1:303890/150‑1 (MQ=255)
                 cTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGc                                                                                                                                      <  1:5198312/150‑1 (MQ=255)
                  ttttACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGcc                                                                                                                                     >  2:3451942/1‑150 (MQ=255)
                       ccGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCGGCCCCGGCGAGCGCGCCTGGGcgggacc                                                                                                                                <  1:6853422/150‑7 (MQ=255)
                       ccGTCGTTTATCAAATGGCGGTGCAACCCATGGATGCCGCACAACCACCACGGACAACTTTTAGCATGCACCGATGCTGGCTCATTTTTCTGCAGACGCTGGGGTTGGCCTTTGCCCGTCCCGCCCAGGCGCGCTCGCCGGGGCCgcagc                                                                                                                                >  2:1734995/1‑150 (MQ=255)
                        cgtcgtTTATCACTTGGCGGTGCACCCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTGTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCgcagca                                                                                                                               <  1:4397592/150‑1 (MQ=255)
                          tcgtTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCaaa                                                                                                                             <  1:3444672/150‑1 (MQ=255)
                          tcgtTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCaaa                                                                                                                             <  1:7884214/150‑1 (MQ=255)
                                 tCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATAcc                                                                                                                      >  1:470213/1‑150 (MQ=255)
                                   aCTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGc                                                                                                                    >  1:4792444/1‑150 (MQ=255)
                                            tGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAgg                                                                                                           <  1:915593/150‑1 (MQ=255)
                                                aCCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCg                                                                                                       >  2:2187561/1‑150 (MQ=255)
                                                 cccATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGt                                                                                                      <  2:3310899/150‑1 (MQ=255)
                                                  ccATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTg                                                                                                     >  2:1817903/1‑150 (MQ=255)
                                                   cATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGa                                                                                                    >  2:3058109/1‑150 (MQ=255)
                                                   cATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGa                                                                                                    <  1:6606654/150‑1 (MQ=255)
                                                    aTGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGAc                                                                                                   <  1:8718327/150‑1 (MQ=255)
                                                      ggATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCt                                                                                                 <  1:8301965/150‑1 (MQ=255)
                                                        aTGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTgg                                                                                               >  2:1986001/1‑150 (MQ=255)
                                                        aTGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTgg                                                                                               >  1:1011389/1‑150 (MQ=255)
                                                            cGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCa                                                                                           <  1:2046164/150‑1 (MQ=255)
                                                            cGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCa                                                                                           <  2:1249458/150‑1 (MQ=255)
                                                                 aCCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTc                                                                                      <  2:1157052/150‑1 (MQ=255)
                                                                   cGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTAGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGc                                                                                    >  1:3777584/1‑150 (MQ=255)
                                                                    gACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATGTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCa                                                                                   <  2:8288532/150‑1 (MQ=255)
                                                                     aCCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCAc                                                                                  <  2:8451873/150‑1 (MQ=255)
                                                                       cACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTGTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCa                                                                                <  2:7257343/150‑1 (MQ=255)
                                                                          ggACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCa                                                                             >  2:871453/1‑150 (MQ=255)
                                                                          ggACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCa                                                                             >  2:8543092/1‑150 (MQ=255)
                                                                          ggACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCa                                                                             >  2:6408309/1‑150 (MQ=255)
                                                                           gACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCaa                                                                            <  1:2187561/150‑1 (MQ=255)
                                                                            acacGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCTCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCaaa                                                                           <  1:546230/150‑1 (MQ=255)
                                                                              acGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGt                                                                         >  1:4867347/1‑150 (MQ=255)
                                                                                gTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCa                                                                                           >  2:1685195/1‑130 (MQ=255)
                                                                                gTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTAc                                                                       >  1:395114/1‑150 (MQ=255)
                                                                                  tttCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCAcc                                                                                 >  1:6088308/1‑138 (MQ=255)
                                                                                  tttCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCAcc                                                                                 <  2:6088308/138‑1 (MQ=255)
                                                                                   ttCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGc                                                                    <  2:4073443/150‑1 (MQ=255)
                                                                                     cGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGa                                                                  <  1:8401652/150‑1 (MQ=255)
                                                                                      gCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAg                                                                 <  2:2408786/150‑1 (MQ=255)
                                                                                            aGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTgg                                                           >  2:8395909/1‑150 (MQ=255)
                                                                                            aGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTgg                                                           >  2:5914640/1‑150 (MQ=255)
                                                                                            aGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTgg                                                           <  1:1296495/150‑1 (MQ=255)
                                                                                            aGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTgg                                                           <  1:5296855/150‑1 (MQ=255)
                                                                                            aGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTgg                                                           <  1:5717819/150‑1 (MQ=255)
                                                                                                 tGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAAc                                                      >  2:2051074/1‑150 (MQ=255)
                                                                                                     ggCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCaga                                                  <  2:8807403/150‑1 (MQ=255)
                                                                                                     ggCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCaga                                                  <  1:335068/150‑1 (MQ=255)
                                                                                                     ggCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCaga                                                  <  2:7202174/150‑1 (MQ=255)
                                                                                                      gCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCagag                                                 >  2:158854/1‑150 (MQ=255)
                                                                                                      gCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCagag                                                 <  2:3942729/150‑1 (MQ=255)
                                                                                                       cTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGc                                                >  2:521767/1‑150 (MQ=255)
                                                                                                        tCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGcc                                               >  2:6529597/1‑150 (MQ=255)
                                                                                                        tCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGcc                                               >  1:7387562/1‑150 (MQ=255)
                                                                                                        tCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGcc                                               >  2:4327162/1‑150 (MQ=255)
                                                                                                        tCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGcc                                               >  1:803284/1‑150 (MQ=255)
                                                                                                        tCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGcc                                               >  2:7772863/1‑150 (MQ=255)
                                                                                                         cATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCg                                              <  2:448011/150‑1 (MQ=255)
                                                                                                                   ctgACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCt                                    <  1:8548819/150‑1 (MQ=255)
                                                                                                                    tgACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTg                                   <  1:6475822/150‑1 (MQ=255)
                                                                                                                    tgACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTg                                   <  2:8066421/150‑1 (MQ=255)
                                                                                                                     gACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTgg                                  <  2:5032661/150‑1 (MQ=255)
                                                                                                                     gACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTgg                                  <  1:8636985/150‑1 (MQ=255)
                                                                                                                      aCGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCTCAGATTCGCACCAGCAAAGTACTTCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggt                                 <  1:5418686/150‑1 (MQ=255)
                                                                                                                      aCGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggt                                 <  2:7407361/150‑1 (MQ=255)
                                                                                                                      aCGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggt                                 <  2:6883712/150‑1 (MQ=255)
                                                                                                                        gCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtga                               <  1:670211/150‑1 (MQ=255)
                                                                                                                        gCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtga                               <  2:6693590/150‑1 (MQ=255)
                                                                                                                        gCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtga                               <  1:4050856/150‑1 (MQ=255)
                                                                                                                        gCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtga                               <  2:1057084/150‑1 (MQ=255)
                                                                                                                         cTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtgag                              <  1:7155390/150‑1 (MQ=255)
                                                                                                                          tGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTggtgagg                             <  2:5496815/150‑1 (MQ=255)
                                                                                                                          tGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCGGGTCAACCAGAGCCGCAACAGTTCTggtgagg                             <  2:2889979/150‑1 (MQ=255)
                                                                                                                              gTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGaa                         <  1:157918/150‑1 (MQ=255)
                                                                                                                                   ccTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTg                    >  1:8275871/1‑150 (MQ=255)
                                                                                                                                    cTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGa                   <  2:5710395/150‑1 (MQ=255)
                                                                                                                                    cTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGa                   <  1:3545739/150‑1 (MQ=255)
                                                                                                                                     tttGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAg                  <  2:2696035/150‑1 (MQ=255)
                                                                                                                                     tttGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCCGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAg                  >  2:8609742/1‑150 (MQ=255)
                                                                                                                                     attGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAg                  <  2:4794212/149‑1 (MQ=255)
                                                                                                                                      ttGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGc                 <  1:6320772/150‑1 (MQ=255)
                                                                                                                                           ggTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTc            <  1:2576024/150‑1 (MQ=255)
                                                                                                                                           ggTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTc            <  1:3279529/150‑1 (MQ=255)
                                                                                                                                           ggTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTc            >  2:7894184/1‑150 (MQ=255)
                                                                                                                                            ctgcggCCCCGGCGAGCGCGCCTGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCg           >  1:7248035/6‑150 (MQ=255)
                                                                                                                                            ctgcggCCCCGGCGAGCGCGCCTGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCg           >  2:5744091/6‑150 (MQ=255)
                                                                                                                                              cccGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGc         <  2:6617511/150‑1 (MQ=255)
                                                                                                                                               ccGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCa        >  2:4908944/1‑150 (MQ=255)
                                                                                                                                                 gCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCGGACAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATc      <  1:1862704/150‑1 (MQ=255)
                                                                                                                                                  cccAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCg     <  1:2861732/150‑1 (MQ=255)
                                                                                                                                                   ccAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGa    >  1:8142951/1‑150 (MQ=255)
                                                                                                                                                   ccAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGa    <  1:1713225/150‑1 (MQ=255)
                                                                                                                                                     aGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGAGt  >  1:5926863/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATGCTGTTGCCGTGGTCTTTTACCGTCGTTTATCACTTGGCGGTGCAACCCATGGATGCCGCTCACCGACCACGGACACGTTTTCGCATGCAGCGATGCTGGCTCATTTTTCTGCTGACGCTGGGGTTGGCCTTTGCCGGTCCCGCCCAGGCGCGCTCGCCGGGGCCGCAGCAAAACATACCGCCTGCCAAGGTGCGTGACCTGGCGCAGATTCGCACCAGCAAAGTACTGCGAGTGCTGGTCAACCAGAGCCGCAACAGTTCTGGTGAGGTGAAGGGTGAGCCAGTCGGCATCGAGT  >  NC_002947/2646642‑2646940

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: