Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,735,206 T→C 6.0% Y258H (TAT→CAT)  mepH → peptidoglycan DD‑endopeptidase MepH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,735,2060TC 6.0% 402.7 / 11.4 132Y258H (TAT→CAT) mepHpeptidoglycan DD‑endopeptidase MepH
Reads supporting (aligned to +/- strand):  ref base T (66/58);  new base C (6/2);  total (72/60)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.16e-01

Only 100 of 133 total aligned reads displayed.

CGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGCTATTGTACGTATGCAAA  >  NC_000913/1735058‑1735355
                                                                                                                                                    |                                                                                                                                                     
cGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTa                                                                                                                                                      <  1:334622/150‑1 (MQ=255)
 ggCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTAt                                                                                                                                                     <  1:41012/150‑1 (MQ=255)
  gCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATg                                                                                                                                                    <  1:353166/150‑1 (MQ=255)
      ccTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTgg                                                                                                                                                >  2:1190604/1‑150 (MQ=255)
      ccTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTgg                                                                                                                                                <  2:2334840/150‑1 (MQ=255)
         ggTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGgcgc                                                                                                                                             >  1:1058931/1‑150 (MQ=255)
           tCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCtc                                                                                                                                           >  1:443633/1‑150 (MQ=255)
            cTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCtcg                                                                                                                                          <  1:1210379/150‑1 (MQ=255)
            cTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCtcg                                                                                                                                          >  2:1145297/1‑150 (MQ=255)
            cTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCAACATAGTGGCGCtgc                                                                                                                                          <  1:2158396/150‑3 (MQ=255)
             tttttCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCtcgt                                                                                                                                         >  2:2419467/1‑150 (MQ=255)
              ttttCCGTACTCAGGGACTCGGCACAGCCGATCATGTCGTCGTGTATGTCGCCAACGGCAAATTTATTGAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCtcgtc                                                                                                                                        <  1:1197650/150‑1 (MQ=255)
               tttCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCtcgtcg                                                                                                                                       >  1:721673/1‑150 (MQ=255)
                ttCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCATTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCgg                                                                                                                                      <  2:1048278/150‑1 (MQ=255)
                ttCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCgg                                                                                                                                      <  2:2255934/150‑1 (MQ=255)
                 tCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCggg                                                                                                                                     <  2:1924103/150‑1 (MQ=255)
                   cGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTa                                                                                                                                   >  2:1109248/1‑150 (MQ=255)
                    gTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTaa                                                                                                                                  >  2:2193946/1‑150 (MQ=255)
                      aCTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATg                                                                                                                                <  1:1232289/150‑1 (MQ=255)
                      aCTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATg                                                                                                                                <  2:1423140/150‑1 (MQ=255)
                      aCTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATg                                                                                                                                >  2:1996976/1‑150 (MQ=255)
                         cAGGGACGCGGCACAGCCGATCATGTCGTCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGAcc                                                                                                                             >  2:452183/1‑150 (MQ=255)
                         cAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCTGGTAATGAcc                                                                                                                             <  2:2058055/150‑1 (MQ=255)
                             gACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCaa                                                                                                                         <  2:1185922/150‑1 (MQ=255)
                             gACGCAGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCACTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCaa                                                                                                                         >  2:1316452/1‑150 (MQ=255)
                              aCGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCaaa                                                                                                                        >  1:1890002/1‑150 (MQ=255)
                               cgcgGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCaaaa                                                                                                                       <  2:243352/150‑1 (MQ=255)
                                 cgGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAac                                                                                                                     >  1:792454/1‑150 (MQ=255)
                                  ggCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAaca                                                                                                                    <  2:488441/150‑1 (MQ=255)
                                        ggcgtgcATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGc                                                                                                              >  2:216433/7‑150 (MQ=255)
                                        gCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGc                                                                                                              <  1:2479092/150‑1 (MQ=255)
                                          cGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTa                                                                                                            >  1:1746294/1‑150 (MQ=255)
                                            aTCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTaaa                                                                                                          <  2:2268241/150‑1 (MQ=255)
                                             tCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTaaaa                                                                                                         >  1:523067/1‑150 (MQ=255)
                                              cATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAAc                                                                                                        <  2:1639624/150‑1 (MQ=255)
                                                tGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACtt                                                                                                      >  2:1982315/1‑150 (MQ=255)
                                                tGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACtt                                                                                                      >  1:1657943/1‑150 (MQ=255)
                                                 gTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACttt                                                                                                     >  2:1857703/1‑150 (MQ=255)
                                                   cGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTAc                                                                                                   <  2:2185639/150‑1 (MQ=255)
                                                   cGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTAc                                                                                                   <  2:2425517/150‑1 (MQ=255)
                                                    ggCGTTTATGTCGGAAAATGCAAATTTATTCAGTCACAGCGCACAGGTCAGGAAATTGAAATCACTTCTCTCATTGAAGACTACTTTCAGCGCCACTATGTTGGCGCTGGTTGGTTAATGACCCCAAAAACACTTCGCTAAAACTTTAcc                                                                                                  <  1:1505980/150‑1 (MQ=255)
                                                    ggCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTAcc                                                                                                  >  2:1098767/1‑150 (MQ=255)
                                                       gtgtATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCtg                                                                                               >  1:961569/1‑150 (MQ=255)
                                                          tATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCtgttg                                                                                            <  1:1885765/150‑1 (MQ=255)
                                                            tGTCGGAAACGAAAAATTTATTCACTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTgttgtt                                                                                          >  2:1808824/1‑150 (MQ=255)
                                                             gTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGCCCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTa                                                                                         >  1:2104762/1‑150 (MQ=255)
                                                                      ggcaaATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAg                                                                                <  2:1468344/150‑1 (MQ=255)
                                                                       gcaaATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAgg                                                                               >  2:1069003/1‑150 (MQ=255)
                                                                        caaATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAggg                                                                              >  1:2200722/1‑150 (MQ=255)
                                                                        caaATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAggg                                                                              <  1:2082374/150‑1 (MQ=255)
                                                                        caaATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAggg                                                                              <  2:1058931/150‑1 (MQ=255)
                                                                           aTTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTaa                                                                           >  1:1062599/1‑150 (MQ=255)
                                                                                      cACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCtttt                                                                <  1:1190604/150‑1 (MQ=255)
                                                                                       aCCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTg                                                               <  2:2371956/150‑1 (MQ=255)
                                                                                       aCCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTg                                                               <  1:807039/150‑1 (MQ=255)
                                                                                         cgcgCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGtc                                                             <  2:2267129/150‑1 (MQ=255)
                                                                                         cgcgCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGtc                                                             >  2:618141/1‑150 (MQ=255)
                                                                                         cgcgCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGtc                                                             <  2:1287285/150‑1 (MQ=255)
                                                                                           cgcACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGtctc                                                           >  1:2213972/1‑150 (MQ=255)
                                                                                                  gTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCtttt                                                    >  2:1245769/1‑150 (MQ=255)
                                                                                                  gTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCtttt                                                    <  1:1787606/150‑1 (MQ=255)
                                                                                                    cAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTaa                                                  >  1:1081017/1‑150 (MQ=255)
                                                                                                      ggAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAAtt                                                <  1:685623/150‑1 (MQ=255)
                                                                                                      ggAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAAtt                                                >  2:2480883/1‑150 (MQ=255)
                                                                                                      ggAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTAATCTTTTGTATCAACTTTTAAtt                                                >  2:2103328/1‑150 (MQ=255)
                                                                                                        aaaTTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGACACTATGTTGTCCCTCGTCGGGTAATGACCCCAAAAACACTTAGCTAAAACTTTACCCTGTTGTTACGGCAACAGCGTAAGTTCATCTTTTCTCTCACCTTTTAATTTg                                              >  1:1385555/1‑150 (MQ=255)
                                                                                                         aaTTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGc                                             <  1:1534702/150‑1 (MQ=255)
                                                                                                         aaTTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGACTCACCTTTTAATTTGc                                             >  2:1320279/1‑150 (MQ=255)
                                                                                                              aaaTCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTAccc                                        <  2:1824888/150‑1 (MQ=255)
                                                                                                              aaaTCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTAccc                                        <  1:1369975/150‑1 (MQ=255)
                                                                                                               aaTCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTATCACCTTTTAATTTGCTACCCt                                       >  1:799183/1‑150 (MQ=255)
                                                                                                                   aCTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCATGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATcc                                   <  1:1645556/150‑1 (MQ=255)
                                                                                                                    cTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCa                                  >  2:1639634/1‑150 (MQ=255)
                                                                                                                          tcAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGca                            >  1:1953833/1‑150 (MQ=255)
                                                                                                                          tcAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGca                            >  1:1144757/1‑150 (MQ=255)
                                                                                                                           cAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGCCCCCAAAAACCCTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGcac                           >  1:1332183/1‑150 (MQ=255)
                                                                                                                           cAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGcac                           >  1:2046964/1‑150 (MQ=255)
                                                                                                                             gTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACaa                         >  1:967380/1‑150 (MQ=255)
                                                                                                                             gTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACaa                         >  2:689728/1‑150 (MQ=255)
                                                                                                                              tGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAAt                        <  1:1808824/150‑1 (MQ=255)
                                                                                                                               gAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCCATCCATACGCACAATa                       >  2:1717547/1‑150 (MQ=255)
                                                                                                                                  gACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAgg                    >  2:2054732/1‑150 (MQ=255)
                                                                                                                                   actactGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGc                   >  1:1379942/1‑150 (MQ=255)
                                                                                                                                   actactGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGc                   <  1:2193946/150‑1 (MQ=255)
                                                                                                                                   actactGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGc                   <  2:289415/150‑1 (MQ=255)
                                                                                                                                   actactGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGc                   <  1:1316452/150‑1 (MQ=255)
                                                                                                                                       cgacgaGCGCCAACATAGTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGCTAtt               >  2:455818/6‑150 (MQ=255)
                                                                                                                                        tGGCAGCGCGACTATGTTGGCTCGCGGCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTTCCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGGTATTg              <  2:701736/150‑1 (MQ=255)
                                                                                                                                         ggCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTGGTTACGGCAACAGGGTAAGTTCATCTGTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAAGAAGGCTATTGt             >  1:1819609/1‑150 (MQ=255)
                                                                                                                                         acgaGCGCCAACATAGTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGCTATTGt             >  2:421401/4‑150 (MQ=255)
                                                                                                                                          gCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGCTATTGTa            <  2:523067/150‑1 (MQ=255)
                                                                                                                                          cgaGCGCCAACATAGTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGCTATTGTa            >  2:1490518/3‑150 (MQ=255)
                                                                                                                                          cgaGCGCCAACATAGTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGCTATTGTa            >  2:69421/3‑150 (MQ=255)
                                                                                                                                            aGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGCTATTGTACg          <  1:2152951/150‑1 (MQ=255)
                                                                                                                                            aGCGCCAACATAGTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGCTATTGTACg          >  1:1543771/1‑150 (MQ=255)
                                                                                                                                             gcgcCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGCTATTGTACGt         >  1:699329/1‑150 (MQ=255)
                                                                                                                                              cgcCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGCTATTGTACGTa        >  2:256448/1‑150 (MQ=255)
                                                                                                                                                ccACTATCTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGATACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAAATTGCTACCCTATCCAAACGCACAATAAGGCTATTGTACGTATg      >  2:23110/1‑150 (MQ=255)
                                                                                                                                                 cACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGCTATTGTACGTATGc     <  2:1081017/150‑1 (MQ=255)
                                                                                                                                                    tATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGTACAATAGCCTTATTGTGCGTATggata  <  1:205177/150‑5 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGGCGACCTGGTCTTTTTCCGTACTCAGGGACGCGGCACAGCCGATCATGTCGGCGTGTATGTCGGCAACGGCAAATTTATTCAGTCACCGCGCACAGGTCAGGAAATTCAAATCACTTCTCTCAGTGAAGACTACTGGCAGCGCCACTATGTTGGCGCTCGTCGGGTAATGACCCCAAAAACACTTCGCTAAAACTTTACCCTGTTGTTACGGCAACAGGGTAAGTTCATCTTTTGTCTCACCTTTTAATTTGCTACCCTATCCATACGCACAATAAGGCTATTGTACGTATGCAAA  >  NC_000913/1735058‑1735355

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: