Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,522,778 C→G 6.1% intergenic (‑27/‑93) hofM ← / → mrcA DNA catabolic putative pilus assembly protein/penicillin‑binding protein 1a, murein transglycosylase and transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,522,7780CG 6.1% 424.1 / 12.6 128intergenic (‑27/‑93)hofM/mrcADNA catabolic putative pilus assembly protein/penicillin‑binding protein 1a, murein transglycosylase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (50/69);  new base G (3/5);  total (54/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 131 total aligned reads displayed.

TTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGC  >  NC_000913/3522630‑3522926
                                                                                                                                                    |                                                                                                                                                    
ttCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCAct                                                                                                                                                     >  1:1627134/1‑150 (MQ=255)
   tcCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCCACCGCTCCCTCTTCGTGCTGTTGTTAATGCACACCCATTTGCCAGCACTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTg                                                                                                                                                  >  1:277163/1‑150 (MQ=255)
   tcCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTg                                                                                                                                                  >  1:207173/1‑150 (MQ=255)
    ccAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGa                                                                                                                                                 <  1:940587/150‑1 (MQ=255)
    ccAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTGGCTGTTGTAAATGCACACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGa                                                                                                                                                 >  1:1234602/1‑150 (MQ=255)
        cGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACggg                                                                                                                                             <  2:1704952/150‑1 (MQ=255)
         ggCAACCGCCACCAGCGTTGCAAAAAGCAATCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGt                                                                                                                                            <  2:1464385/150‑1 (MQ=255)
            aaCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCAGAGTGACGGGtgat                                                                                                                                         <  1:376718/150‑4 (MQ=255)
            aaCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCAGAGTGACGGGtgat                                                                                                                                         <  2:742759/150‑4 (MQ=255)
            aaCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTAACACCCGTCACTCTGACGGGtata                                                                                                                                         <  1:1631732/150‑1 (MQ=255)
             aCCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTatat                                                                                                                                        >  2:774219/1‑150 (MQ=255)
              ccGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACTGGTatata                                                                                                                                       >  1:253791/1‑149 (MQ=255)
                 ccaccaGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAAt                                                                                                                                    >  2:1383368/1‑150 (MQ=255)
                  caccaGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAAtt                                                                                                                                   <  2:477659/150‑2 (MQ=255)
                  caccaGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATg                                                                                                                                   <  2:359204/150‑1 (MQ=255)
                    ccaGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCg                                                                                                                                 <  1:112501/150‑1 (MQ=255)
                      aGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTc                                                                                                                               >  1:1702475/1‑150 (MQ=255)
                       gCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCt                                                                                                                              <  2:630834/150‑1 (MQ=255)
                       gCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGAAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCt                                                                                                                              <  2:1554256/150‑1 (MQ=255)
                        cGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTg                                                                                                                             <  1:1144058/150‑1 (MQ=255)
                          ttGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCATTTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCGATGCGTCTGGc                                                                                                                           <  2:98529/150‑1 (MQ=255)
                           tGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCt                                                                                                                          <  1:1090731/150‑1 (MQ=255)
                            gCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCtt                                                                                                                         <  1:1480271/150‑1 (MQ=255)
                             cAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCCTCTGGCTTg                                                                                                                        >  2:1140655/1‑150 (MQ=255)
                               aaaaGCATTCTTTTGCGCCCCGCACGATCGCAACCCCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGcc                                                                                                                      <  2:1638718/150‑1 (MQ=255)
                                  aGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCttt                                                                                                                   <  2:1147314/150‑1 (MQ=255)
                                    cATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTtat                                                                                                                 >  1:1589841/1‑150 (MQ=255)
                                    cATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTtat                                                                                                                 <  2:458685/150‑1 (MQ=255)
                                     aTTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGcc                                                                                                                      <  1:235626/144‑1 (MQ=255)
                                                 cccGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtt                                                                                                    >  2:282132/1‑150 (MQ=255)
                                                 cccGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtt                                                                                                    >  1:1258254/1‑150 (MQ=255)
                                                 cccGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtt                                                                                                    >  1:110293/1‑150 (MQ=255)
                                                    gCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgt                                                                                                 <  1:756814/150‑1 (MQ=255)
                                                      aCGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgttt                                                                                               <  2:360222/150‑1 (MQ=255)
                                                      aCGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgttt                                                                                               <  2:513522/150‑1 (MQ=255)
                                                       cGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTta                                                                                              <  1:1005171/150‑1 (MQ=255)
                                                         aTCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTtata                                                                                            <  2:1433345/150‑1 (MQ=255)
                                                            gCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAAc                                                                                         >  1:1525375/1‑150 (MQ=255)
                                                            gCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAAc                                                                                         >  2:1020772/1‑150 (MQ=255)
                                                            gCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAAc                                                                                         <  1:1169718/150‑1 (MQ=255)
                                                            gCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAAc                                                                                         >  1:758394/1‑150 (MQ=255)
                                                             cAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGAGGGGTATATCAATGCGTCTGGCTTGCCTTTATACTAACGCGCGTTTGTTTATAAACt                                                                                        >  2:68176/1‑150 (MQ=255)
                                                               accgctaccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGc                                                                                      <  1:411457/150‑1 (MQ=255)
                                                               accgctaccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGc                                                                                      >  1:1477264/1‑150 (MQ=255)
                                                               accgctaccgctTCTTGCTGTTGTAAATGCAAACCAATTCGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGc                                                                                      >  1:1733759/1‑150 (MQ=255)
                                                                ccgctaccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGcc                                                                                     >  2:799156/1‑150 (MQ=255)
                                                                 cgctaccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGccc                                                                                    >  1:634299/1‑150 (MQ=255)
                                                                  gctaccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATGGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCa                                                                                   <  2:1239201/150‑1 (MQ=255)
                                                                   ctaccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCgtttgtt                                                                                                <  2:1476628/136‑1 (MQ=255)
                                                                     accgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAAt                                                                                >  1:208583/1‑150 (MQ=255)
                                                                      ccgctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATg                                                                               <  2:487349/150‑1 (MQ=255)
                                                                        gctTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGaa                                                                             >  1:1376832/1‑150 (MQ=255)
                                                                         cttcttGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGaaa                                                                            <  1:204447/150‑1 (MQ=255)
                                                                          ttcttGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAATTGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAAc                                                                           <  2:492326/150‑1 (MQ=255)
                                                                          ttcttGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTTTCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAAc                                                                           <  1:417983/150‑1 (MQ=255)
                                                                           tcttGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACt                                                                          <  2:207173/150‑1 (MQ=255)
                                                                                cTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATg                                                                     >  1:531487/1‑150 (MQ=255)
                                                                                  gttgtAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATggg                                                                   <  2:731800/150‑1 (MQ=255)
                                                                                  gttgtAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATggg                                                                   >  1:713777/1‑150 (MQ=255)
                                                                                  gttgtAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGACTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATggg                                                                   >  1:720092/1‑150 (MQ=255)
                                                                                    tgtAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCATTTTACTACCGCGCGGTTGTTTATAAAAATCCCAAAATAAACTAAATTGGaaa                                                                >  2:1157195/1‑150 (MQ=255)
                                                                                    tgtAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTATAACAGATATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGaa                                                                 >  2:233258/1‑150 (MQ=255)
                                                                                      tAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAAt                                                               <  2:445112/150‑1 (MQ=255)
                                                                                       aaaTGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAAtt                                                              >  2:1044155/1‑150 (MQ=255)
                                                                                          tGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTcc                                                           >  1:527490/1‑150 (MQ=255)
                                                                                            cAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGTTGTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAg                                                         <  1:135147/150‑1 (MQ=255)
                                                                                            cAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAg                                                         <  2:1110265/150‑1 (MQ=255)
                                                                                             aaaCCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGt                                                        <  2:1531464/150‑1 (MQ=255)
                                                                                                    attGCCAGAGCATAACTGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTc                                                 <  2:333941/149‑1 (MQ=255)
                                                                                                     ttGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCg                                                <  1:1062718/150‑1 (MQ=255)
                                                                                                       gCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCAAGTGAAGTTCGTa                                              >  1:868663/1‑150 (MQ=255)
                                                                                                        ccAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTaa                                             >  2:246834/1‑150 (MQ=255)
                                                                                                         cAGATCTTAAATGCCAGTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATTAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTaaa                                            >  2:147836/1‑150 (MQ=255)
                                                                                                          aGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAg                                           <  1:1093396/150‑1 (MQ=255)
                                                                                                          aGATCTTAAAAGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAg                                           <  2:86777/150‑1 (MQ=255)
                                                                                                             tCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTAt                                        >  1:914266/1‑150 (MQ=255)
                                                                                                                 aaaTGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGATTGTTTATAAACTGCCCAAATGAAACTAAATGGTAAATTTCCAGTGAAGTTCGTAAAGTAttttt                                    >  2:368987/1‑150 (MQ=255)
                                                                                                                  agtgcggttGTGAGGATCTCCTTATCACCGGTCGGTGTGACGGGTATATCAATGCGTCTGGCTTGGCTTTGTACGGCCGCGGGTGTGTTTATAAGCTGCCCAAATGAGGGTAGGTGGGATATTTCGTGTGATGTTCGTATTGTAtttttt                                   >  2:439676/8‑150 (MQ=255)
                                                                                                                     gCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGAt                                <  2:1525375/150‑1 (MQ=255)
                                                                                                                     gCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGAt                                <  1:1725264/150‑1 (MQ=255)
                                                                                                                      ccATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATc                               >  1:134493/1‑150 (MQ=255)
                                                                                                                       cATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATcc                              <  1:664090/150‑1 (MQ=255)
                                                                                                                        aTTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCt                             <  2:1376832/150‑1 (MQ=255)
                                                                                                                           tttattatCTCCTTATCACCCGTCACTCTGACGTGTATATCAATGCGTCTTGCTTGCCTTTATAATACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGc                          <  1:1211171/144‑1 (MQ=255)
                                                                                                                           gtgATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTCCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGc                          <  2:956868/150‑1 (MQ=255)
                                                                                                                           gtgATGATCTCCTTATAACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGc                          <  1:1056468/150‑1 (MQ=255)
                                                                                                                             gatgatCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAAAGGGAAATTTCCAGTGAAGTTCGTATAGTATTTTTTGATCCTTGCAg                        >  2:705312/1‑150 (MQ=255)
                                                                                                                              atgatCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGt                       <  2:758394/150‑1 (MQ=255)
                                                                                                                                   ctcCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATTAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGtt                  >  2:1184542/1‑150 (MQ=255)
                                                                                                                                    tcCTTATCACCCGTCACTCTGACGGGAATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAACCTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTg                 <  2:1274875/150‑1 (MQ=255)
                                                                                                                                    accTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCGTTTATACGACCGCGCGTTTGTTTATTAACTGCCCAAATGAAACTAAATGGGAAAGTTCCAGTGAAGTTCGTAATGTATTTTTTGATCCTTGCAGTCTGTTg                 >  1:611144/2‑150 (MQ=255)
                                                                                                                                       ttATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCAt              <  1:799694/150‑1 (MQ=255)
                                                                                                                                         aTCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTc            >  1:1045752/1‑150 (MQ=255)
                                                                                                                                           taCCCGTCAGAGTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTctg          >  1:93404/2‑150 (MQ=255)
                                                                                                                                           cACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTctg          >  2:498911/1‑150 (MQ=255)
                                                                                                                                           cACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTctg          >  1:1466051/1‑150 (MQ=255)
                                                                                                                                            aCCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTctgc         <  1:793234/150‑1 (MQ=255)
                                                                                                                                            aCCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTctgc         <  1:537114/150‑1 (MQ=255)
                                                                                                                                              ccGTCACTCTGACGGGTATATCAATGAGTCTGGATTTCCTTTATTCTACCGCGCGTTTTTTTATAAACTGCCCAAATGAAACTAAATGTGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGTTCCTTGCAGTCTGTTGCATTctgctg       <  1:308115/150‑1 (MQ=255)
                                                                                                                                                  cACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAg   <  2:1271571/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCTCCAGCGGCAACCGCCACCAGCGTTGCAAAAAGCATTCTTTTGCGCCCCGCACGATCGCAACCGCTACCGCTTCTTGCTGTTGTAAATGCAAACCAATTTGCCAGATCTTAAATGCCATTGTGATGATCTCCTTATCACCCGTCACTCTGACGGGTATATCAATGCGTCTGGCTTGCCTTTATACTACCGCGCGTTTGTTTATAAACTGCCCAAATGAAACTAAATGGGAAATTTCCAGTGAAGTTCGTAAAGTATTTTTTGATCCTTGCAGTCTGTTGCATTCTGCTGGGAGC  >  NC_000913/3522630‑3522926

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