Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,522,780 C→G 7.0% intergenic (‑29/‑91) hofM ← / → mrcA DNA utilization protein HofM/peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,522,7800CG 7.0% 330.4 / 12.2 100intergenic (‑29/‑91)hofM/mrcADNA utilization protein HofM/peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcA
Reads supporting (aligned to +/- strand):  ref base C (55/38);  new base G (4/3);  total (59/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

Only 100 of 103 total aligned reads displayed.

CAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGG  >  NC_000913/3522635‑3522929
                                                                                                                                                 |                                                                                                                                                     
cAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGAc                                                                                                                                                   <  1:2019/150‑1 (MQ=255)
cAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGAc                                                                                                                                                   >  1:1605673/1‑150 (MQ=255)
   cGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACggg                                                                                                                                                <  1:325587/150‑1 (MQ=255)
     gCAACCGCCACCAGCGTTGCAAAAATCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTTCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGta                                                                                                                                              <  1:427306/150‑1 (MQ=255)
     gCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCAGAGTGACGGGtg                                                                                                                                              <  2:1226102/150‑2 (MQ=255)
     gCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGta                                                                                                                                              >  1:1588649/1‑150 (MQ=255)
      cAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCAGAGTGACGGGtga                                                                                                                                             <  2:651562/150‑3 (MQ=255)
         ccGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTatata                                                                                                                                          >  2:12287/1‑149 (MQ=255)
         ccGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATc                                                                                                                                          >  1:805182/1‑150 (MQ=255)
           gCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTagatgga                                                                                                                                        >  1:703941/1‑147 (MQ=255)
             caccaGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATg                                                                                                                                      <  2:1614650/150‑1 (MQ=255)
                       gCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCtt                                                                                                                            >  2:1313125/1‑150 (MQ=255)
                         aaaaaGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGc                                                                                                                          >  2:1042729/1‑150 (MQ=255)
                              gCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTta                                                                                                                     <  2:206313/150‑1 (MQ=255)
                               cATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTtat                                                                                                                    >  2:175696/1‑150 (MQ=255)
                               cATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTtat                                                                                                                    >  2:1205585/1‑150 (MQ=255)
                                aTTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTtata                                                                                                                   <  1:1548908/150‑1 (MQ=255)
                                   cTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTAtacta                                                                                                                >  2:524592/1‑150 (MQ=255)
                                    ttttGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCGCTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTAtactac                                                                                                               >  1:869368/1‑150 (MQ=255)
                                      ttGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGCTCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACcg                                                                                                             >  2:354750/1‑150 (MQ=255)
                                      ttGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACcg                                                                                                             <  1:1117811/150‑1 (MQ=255)
                                       tGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGAGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACcgc                                                                                                            >  2:1555420/1‑150 (MQ=255)
                                           ccccGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgt                                                                                                        <  2:495498/150‑1 (MQ=255)
                                           ccccGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgt                                                                                                        >  2:72793/1‑150 (MQ=255)
                                              cGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttg                                                                                                     <  1:190925/150‑1 (MQ=255)
                                                cACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgtt                                                                                                   <  1:1587424/150‑1 (MQ=255)
                                                 aCGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgttt                                                                                                  <  2:1620373/150‑1 (MQ=255)
                                                 aCGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgttt                                                                                                  <  1:785689/150‑1 (MQ=255)
                                                  ccaacacAACCCCTACCCCATCATGCTATTGTAAATGCAAACCAATTTGCCAGAACATAAATGCCATAGTGATGATCACCTTATCACCCGTCACTCTGACGGGAATATCAATGCGTCTGGCATGCCTTTATACTACCGCGCGTTTGTTta                                                                                                 >  2:875909/7‑150 (MQ=255)
                                                  cGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTta                                                                                                 >  1:22905/1‑150 (MQ=255)
                                                  cGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTta                                                                                                 >  1:876012/1‑150 (MQ=255)
                                                  cGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTta                                                                                                 >  2:695459/1‑150 (MQ=255)
                                                   gATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTtat                                                                                                >  2:1406607/1‑150 (MQ=255)
                                                    aTCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTtata                                                                                               <  1:818441/150‑1 (MQ=255)
                                                       gCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAAc                                                                                            >  1:1467900/1‑150 (MQ=255)
                                                       gCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAAc                                                                                            >  2:667697/1‑150 (MQ=255)
                                                       gCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGAGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAAc                                                                                            >  2:302879/1‑150 (MQ=255)
                                                        cAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACt                                                                                           >  2:466525/1‑150 (MQ=255)
                                                            cgctaccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGccc                                                                                       <  2:300358/150‑1 (MQ=255)
                                                            cgctaccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGccc                                                                                       >  2:25628/1‑150 (MQ=255)
                                                            cgctaccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGccc                                                                                       <  2:896201/150‑1 (MQ=255)
                                                              ctaccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGACTTTATACTACCGCGCGTTTGTTTATAAACTGCCCaa                                                                                     >  2:1140836/1‑150 (MQ=255)
                                                                 ccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATg                                                                                  <  2:587854/150‑1 (MQ=255)
                                                                 ccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATg                                                                                  >  1:1200501/1‑150 (MQ=255)
                                                                  cgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGa                                                                                 >  2:553510/1‑150 (MQ=255)
                                                                   gctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGaa                                                                                <  1:557487/150‑1 (MQ=255)
                                                                    cttcttGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGTCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGaaa                                                                               <  1:524592/150‑1 (MQ=255)
                                                                         tGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTaaa                                                                          <  1:553510/150‑1 (MQ=255)
                                                                           cTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATg                                                                        <  1:46791/150‑1 (MQ=255)
                                                                               tgtAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGaa                                                                    >  2:412375/1‑150 (MQ=255)
                                                                               tgtAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGaa                                                                    >  1:1430599/1‑150 (MQ=255)
                                                                                 tAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAAt                                                                  >  1:584280/1‑150 (MQ=255)
                                                                                 tAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAAt                                                                  >  1:576129/1‑150 (MQ=255)
                                                                                   aaTGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAAttt                                                                >  2:1621070/1‑150 (MQ=255)
                                                                                       cAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATAACCCGACACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCat                                                            >  2:1443252/1‑149 (MQ=255)
                                                                                         aaCCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTg                                                          >  1:1223621/1‑150 (MQ=255)
                                                                                          aCCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGa                                                         <  1:1453841/150‑1 (MQ=255)
                                                                                            cAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAg                                                       <  2:1381158/150‑1 (MQ=255)
                                                                                             aaTTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGt                                                      >  2:1431612/1‑150 (MQ=255)
                                                                                              aTTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGtt                                                     <  1:175696/150‑1 (MQ=255)
                                                                                                  gCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTa                                                 >  1:880700/1‑150 (MQ=255)
                                                                                                  gCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTa                                                 <  1:1313125/150‑1 (MQ=255)
                                                                                                   ccAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTaa                                                >  1:158019/1‑150 (MQ=255)
                                                                                                   ccAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTaa                                                >  2:556547/1‑150 (MQ=255)
                                                                                                      gATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGt                                             <  2:1111178/150‑1 (MQ=255)
                                                                                                        tCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTAt                                           >  1:1256369/1‑150 (MQ=255)
                                                                                                                  cATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATcc                                 <  2:1588649/150‑1 (MQ=255)
                                                                                                                      gtgATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCg                                                   >  1:25121/1‑128 (MQ=255)
                                                                                                                      gtgATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCg                                                   <  2:25121/128‑1 (MQ=255)
                                                                                                                      gtgATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGc                             <  2:1467900/150‑1 (MQ=255)
                                                                                                                      gtgATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGc                             >  1:59811/1‑150 (MQ=255)
                                                                                                                      gtgATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGc                             >  1:19072/1‑150 (MQ=255)
                                                                                                                      gtgATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGc                             >  2:680242/1‑150 (MQ=255)
                                                                                                                           gatCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCt                        >  2:701175/1‑150 (MQ=255)
                                                                                                                             tctcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGt                      >  1:117436/1‑150 (MQ=255)
                                                                                                                              ctcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGc                             <  1:3023/142‑1 (MQ=255)
                                                                                                                              ctcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGc                             >  2:3023/1‑142 (MQ=255)
                                                                                                                              ctcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGtt                     >  2:892177/1‑150 (MQ=255)
                                                                                                                               tcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTg                    <  2:805182/150‑1 (MQ=255)
                                                                                                                               tcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTg                    >  1:303318/1‑150 (MQ=255)
                                                                                                                                ccTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGc                   <  2:30337/150‑1 (MQ=255)
                                                                                                                                ccTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGc                   >  1:1030081/1‑150 (MQ=255)
                                                                                                                                ccTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGc                   >  1:1127882/1‑150 (MQ=255)
                                                                                                                                  ttATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCAt                 <  1:1242173/150‑1 (MQ=255)
                                                                                                                                   tATCACCCGGCACTCTGACGGGTATATCAATGCGGCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGGCTGTTGCAtt                <  2:555661/150‑1 (MQ=255)
                                                                                                                                    tataCCCGTCAGAGTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTc               >  1:954920/4‑150 (MQ=255)
                                                                                                                                    aTCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTc               <  1:1042729/150‑1 (MQ=255)
                                                                                                                                     tCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTct              >  2:585164/1‑150 (MQ=255)
                                                                                                                                     ataCCCGTCAGAGTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTct              >  2:313208/3‑150 (MQ=255)
                                                                                                                                     ataCCCGTCAGAGTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTct              >  2:1513433/3‑150 (MQ=255)
                                                                                                                                     ataCCCGTCAGAGTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTct              >  1:1270975/3‑150 (MQ=255)
                                                                                                                                       aCCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTctgc            >  1:49127/1‑150 (MQ=255)
                                                                                                                                         ccGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCATTGCAGTCTGTTGCATTctgctg          >  2:450674/1‑150 (MQ=255)
                                                                                                                                          cGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTgg         <  1:1080687/150‑1 (MQ=255)
                                                                                                                                           gTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTggg        >  1:756840/1‑150 (MQ=255)
                                                                                                                                             cACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAg      <  1:861974/150‑1 (MQ=255)
                                                                                                                                              aCTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGc     >  1:1391637/1‑150 (MQ=255)
                                                                                                                                               ctctGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCa    >  1:98458/1‑150 (MQ=255)
                                                                                                                                               ctctGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCa    <  1:531628/150‑1 (MQ=255)
                                                                                                                                                 ctGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAgg  <  1:1431612/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGCAGG  >  NC_000913/3522635‑3522929

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: