Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 165,039 G→C 6.1% A104P (GCG→CCG)  mrcB → peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913165,0390GC 6.1% 264.6 / 4.9 78A104P (GCG→CCG) mrcBpeptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base C (2/3);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGAG  >  NC_000913/164894‑165188
                                                                                                                                                 |                                                                                                                                                     
aggGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGc                                                                                                                                                   >  1:1264435/1‑150 (MQ=255)
      tAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTa                                                                                                                                             >  1:600394/1‑150 (MQ=255)
       aaGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTAt                                                                                                                                            >  2:1373690/1‑150 (MQ=255)
       aaGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCGCAGGCAGttgcc                                                                                                                                            <  2:1936292/150‑4 (MQ=255)
        aGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCGCAGGCAGttgcca                                                                                                                                           <  1:774998/150‑5 (MQ=255)
         gCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCGCAGGCAGttgccag                                                                                                                                          <  2:1381555/150‑6 (MQ=255)
           cTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGcc                                                                                                                                        >  2:1759904/1‑150 (MQ=255)
           cTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGcc                                                                                                                                        <  1:1062972/150‑1 (MQ=255)
               tGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAAt                                                                                                                                    <  1:1999364/150‑1 (MQ=255)
               tGGCAAACGAGGCTGGCTATGGCTACTGCTAAAACTTGATATCGTTTTTGCCGTGCTGATCAACATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTAATGGCCGAAt                                                                                                                                    >  2:1876722/1‑150 (MQ=255)
                  cAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGt                                                                                                                                 >  2:855330/1‑150 (MQ=255)
                    aaCGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCa                                                                                                                               >  1:1061979/1‑150 (MQ=255)
                     aCGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCaa                                                                                                                              <  1:1163311/150‑1 (MQ=255)
                      cgcgGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAAt                                                                                                                             <  1:1531865/150‑1 (MQ=255)
                        cgGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCt                                                                                                                           >  1:879046/1‑150 (MQ=255)
                            tggctatggctaCTGCTAAAACTGGCTATCTTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCACTCTTGAg                                                                                                                       <  1:1917731/150‑1 (MQ=255)
                            tggctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCACTCTTGAg                                                                                                                       <  1:1910562/150‑1 (MQ=255)
                            tggctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAg                                                                                                                       <  1:1373690/150‑1 (MQ=255)
                            tggctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGGTTATGGCCGAATGGTCAATCTTGAg                                                                                                                       <  2:1995950/150‑1 (MQ=255)
                             ggctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGc                                                                                                                      <  1:954337/150‑1 (MQ=255)
                             ggctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGc                                                                                                                      <  1:316878/150‑1 (MQ=255)
                               ctatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCa                                                                                                                    <  1:1617277/150‑1 (MQ=255)
                                tatggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAg                                                                                                                   <  2:640942/150‑1 (MQ=255)
                                   ggctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACa                                                                                                                >  2:475549/1‑150 (MQ=255)
                                    gctaCTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACAt                                                                                                               >  1:1055618/1‑150 (MQ=255)
                                         tGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCa                                                                                                          <  2:1264435/150‑1 (MQ=255)
                                          gCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCAt                                                                                                         >  2:1502297/1‑150 (MQ=255)
                                          gCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCAt                                                                                                         >  2:1203593/1‑150 (MQ=255)
                                              aaaCTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGc                                                                                                     >  1:1332466/1‑150 (MQ=255)
                                                   ggCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGaa                                                                                                >  1:1486825/1‑150 (MQ=255)
                                                      tATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACga                                                                                             <  2:966303/150‑1 (MQ=255)
                                                          gTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATg                                                                                         >  1:1923936/1‑150 (MQ=255)
                                                           tttttGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGACAGACATGACCATCAGCAAGAACGAGATgg                                                                                        <  2:706682/150‑1 (MQ=255)
                                                             tttGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTg                                                                                      <  2:1676106/150‑1 (MQ=255)
                                                               tGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGaa                                                                                    >  2:1334041/1‑150 (MQ=255)
                                                                  cGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAgct                                                                                 <  2:1022039/150‑1 (MQ=255)
                                                                       tGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTgg                                                                            <  1:1355419/150‑1 (MQ=255)
                                                                         aTCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAg                                                                          <  2:976474/150‑1 (MQ=255)
                                                                         aTCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAg                                                                          <  2:1703617/150‑1 (MQ=255)
                                                                          tCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAgg                                                                         <  1:855330/150‑1 (MQ=255)
                                                                           cGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGc                                                                        <  2:1923936/150‑1 (MQ=255)
                                                                           cGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGc                                                                        <  2:1776857/150‑1 (MQ=255)
                                                                           cGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGc                                                                        >  2:479408/1‑150 (MQ=255)
                                                                             ccATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGAc                                                                                                            >  1:630285/1‑112 (MQ=255)
                                                                             ccATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGAc                                                                                                            <  2:630285/112‑1 (MQ=255)
                                                                                tttACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGAccc                                                                   <  2:879046/150‑1 (MQ=255)
                                                                                 ttACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCa                                                                  <  1:353837/150‑1 (MQ=255)
                                                                                  tACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAg                                                                 >  1:1964617/1‑150 (MQ=255)
                                                                                    cGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTa                                                               >  1:354795/1‑150 (MQ=255)
                                                                                      gCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCGGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATc                                                             <  2:632674/150‑1 (MQ=255)
                                                                                          ttATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCa                                                         >  2:55717/1‑150 (MQ=255)
                                                                                            aTCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAgg                                                       <  1:273936/150‑1 (MQ=255)
                                                                                            aTCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAgg                                                       >  2:1321707/1‑150 (MQ=255)
                                                                                            aTCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAgg                                                       <  1:1759904/150‑1 (MQ=255)
                                                                                               tcGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGgtgt                                                    >  1:1159688/1‑150 (MQ=255)
                                                                                               tcGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGACTGGCAACTGCCTCCGGCAGTTAATGGCCGAATGGTCAAACTGGCGCCAGACATGACCATCAGCAAGAACGAGATGGAGAAGCTGCTGCAGGCGACCCAGTATAGTCAGgtgt                                                    >  1:1291540/1‑150 (MQ=255)
                                                                                                cGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTc                                                   <  1:1321707/150‑1 (MQ=255)
                                                                                                     aaaaaaTTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGCGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGCCCCAGTATCGTCAGGTGTCGaaaa                                              >  2:426448/1‑150 (MQ=255)
                                                                                                      aaaaaTTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAAt                                             >  1:398728/1‑150 (MQ=255)
                                                                                                        aaaTTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGa                                           <  1:1334041/150‑1 (MQ=255)
                                                                                                          aTTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGAcc                                         <  1:1605975/150‑1 (MQ=255)
                                                                                                          aTTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGAcc                                         <  2:1950243/150‑1 (MQ=255)
                                                                                                               tAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTcc                                    <  2:845223/150‑1 (MQ=255)
                                                                                                                 gCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCAAGTATCGTCAGGTGTCGAAAATGACCCGTCCTg                                  >  2:1935424/1‑150 (MQ=255)
                                                                                                                   cGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGc                                >  1:70123/1‑150 (MQ=255)
                                                                                                                   cGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGc                                >  1:1205442/1‑150 (MQ=255)
                                                                                                                    gTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCg                               <  2:1055618/150‑1 (MQ=255)
                                                                                                                       ttGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAAt                            <  2:1061979/150‑1 (MQ=255)
                                                                                                                        tGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAAtt                           >  1:1773053/1‑150 (MQ=255)
                                                                                                                             gCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCg                      <  2:1773053/150‑1 (MQ=255)
                                                                                                                             gCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCg                      >  2:686677/1‑150 (MQ=255)
                                                                                                                                   tCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAgg                <  2:1964617/150‑1 (MQ=255)
                                                                                                                                   tCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAgg                >  2:1973754/1‑150 (MQ=255)
                                                                                                                                    ccataaaCTGCCGCAGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGc               >  2:1110062/6‑150 (MQ=255)
                                                                                                                                     tGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGcc              >  1:911853/1‑150 (MQ=255)
                                                                                                                                     tGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGcc              >  1:1353463/1‑150 (MQ=255)
                                                                                                                                      ggCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCa             >  2:1669641/1‑150 (MQ=255)
                                                                                                                                         aaCTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGTCCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACa          <  1:32407/150‑1 (MQ=255)
                                                                                                                                           ctgcctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGc        >  1:347190/1‑150 (MQ=255)
                                                                                                                                              cctgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCAtt     >  2:1831854/1‑150 (MQ=255)
                                                                                                                                                tgcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTga   <  2:354795/150‑1 (MQ=255)
                                                                                                                                                 gcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTgag  >  1:1735624/1‑150 (MQ=255)
                                                                                                                                                 gcGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTgag  >  2:842315/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
AGGGCGTAAGCCTCGTGGCAAACGCGGCTGGCTATGGCTACTGCTAAAACTGGCTATCGTTTTTGCCGTGCTGATCGCCATTTACGGCGTTTATCTCGATCAAAAAATTCGTAGCCGTATTGATGGCAAGGTCTGGCAACTGCCTGCGGCAGTTTATGGCCGAATGGTCAATCTTGAGCCAGACATGACCATCAGCAAGAACGAGATGGTGAAGCTGCTGGAGGCGACCCAGTATCGTCAGGTGTCGAAAATGACCCGTCCTGGCGAATTTACCGTGCAGGCCAACAGCATTGAG  >  NC_000913/164894‑165188

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: