Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,535,612 A→C 10.4% I86S (ATC→AGC)  envZ ← sensory histidine kinase EnvZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,535,6120AC10.4% 248.8 / 17.2 77I86S (ATC→AGC) envZsensory histidine kinase EnvZ
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base C (3/5);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.65e-01

GGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACATACGCACT  >  NC_000913/3535514‑3535710
                                                                                                  |                                                                                                  
ggCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAgagagaga                                                                                                    >  1:515022/1‑99 (MQ=255)
 gccgccCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGAt                                                                                                   <  2:1602803/99‑1 (MQ=255)
 gccgccCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGAt                                                                                                   <  2:1318397/99‑1 (MQ=255)
 gccgccCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGAt                                                                                                   <  2:822579/99‑1 (MQ=255)
    gccCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGCTccc                                                                                                >  2:219628/1‑99 (MQ=255)
     cccAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATcccc                                                                                               >  2:1152017/1‑99 (MQ=255)
      ccAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCa                                                                                              <  2:1458721/99‑1 (MQ=255)
       cAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGCTCCCCAg                                                                                             >  1:289645/1‑99 (MQ=255)
       cAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAg                                                                                             <  1:2366022/99‑1 (MQ=255)
        aGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGc                                                                                            >  1:313690/1‑99 (MQ=255)
            gctgcGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCAc                                                                                        <  2:91655/99‑1 (MQ=255)
             ctgcGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACg                                                                                       >  1:1218390/1‑99 (MQ=255)
             ctgcGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACg                                                                                       >  1:911543/1‑99 (MQ=255)
             ctgcGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACg                                                                                       >  1:1945014/1‑99 (MQ=255)
              tgcGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACgg                                                                                      <  2:1881513/99‑1 (MQ=255)
              tgcGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACgg                                                                                      <  1:1212046/99‑1 (MQ=255)
                   cATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGCTCCCCAGCTCACGGTAGAt                                                                                 <  1:1846014/99‑1 (MQ=255)
                   cATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGAt                                                                                 <  1:2180724/99‑1 (MQ=255)
                    aTCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGAtc                                                                                <  1:755385/99‑1 (MQ=255)
                    aTCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGAtc                                                                                <  1:1258582/99‑1 (MQ=255)
                       tGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTcc                                                                             >  1:210748/1‑99 (MQ=255)
                       tGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTcc                                                                             >  2:1027385/1‑99 (MQ=255)
                         aTGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCg                                                                           >  2:1322725/1‑99 (MQ=255)
                         aTGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCg                                                                           >  1:22804/1‑99 (MQ=255)
                          tGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGa                                                                          >  1:2392328/1‑99 (MQ=255)
                          tGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGa                                                                          >  2:337097/1‑99 (MQ=255)
                          tGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGa                                                                          >  2:427082/1‑99 (MQ=255)
                           gtcTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGAc                                                                         <  2:471998/97‑1 (MQ=255)
                            gCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACg                                                                        >  2:198456/1‑99 (MQ=255)
                             cTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACgg                                                                       >  2:1518471/1‑99 (MQ=255)
                             cTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACgg                                                                       >  2:1732763/1‑99 (MQ=255)
                             cTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACgg                                                                       >  1:986660/1‑99 (MQ=255)
                                aaGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGaaa                                                                    <  1:2169866/99‑1 (MQ=255)
                                     ttCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCggg                                                               >  1:1017381/1‑99 (MQ=255)
                                      tCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGCTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGa                                                              <  2:185178/99‑1 (MQ=255)
                                       cATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAg                                                             >  2:306054/1‑99 (MQ=255)
                                       cATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAg                                                             <  2:2124471/99‑1 (MQ=255)
                                        aTAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAgg                                                            <  2:628936/99‑1 (MQ=255)
                                        aTAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAgg                                                            <  1:1670986/99‑1 (MQ=255)
                                            tgtTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGGCGGAAAGCGGGAGGcaca                                                        >  1:2197032/1‑99 (MQ=255)
                                               tGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCAcaacc                                                     <  2:1068330/99‑1 (MQ=255)
                                                 cgcCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCAcaaccaa                                                   <  2:405068/99‑1 (MQ=255)
                                                 cgcCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCAcaaccaa                                                   <  2:165878/99‑1 (MQ=255)
                                                 cgcCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCAcaaccaa                                                   <  2:1522590/99‑1 (MQ=255)
                                                  gcCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCAcaaccaac                                                  >  2:1117455/1‑99 (MQ=255)
                                                      aaCGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGCTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCg                                              <  2:342592/99‑1 (MQ=255)
                                                       aCGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGt                                             >  1:1677416/1‑99 (MQ=255)
                                                       aCGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGt                                             <  2:1239496/99‑1 (MQ=255)
                                                         gCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGc                                           <  1:1837440/99‑1 (MQ=255)
                                                           aGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGCTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCg                                         >  1:244824/1‑99 (MQ=255)
                                                             aCCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGtc                                       >  1:480889/1‑99 (MQ=255)
                                                             aCCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGtc                                       >  1:2472175/1‑99 (MQ=255)
                                                               cTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGtcct                                     >  2:847795/1‑99 (MQ=255)
                                                                tGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGtcctc                                    >  1:2439047/1‑99 (MQ=255)
                                                                   ctctTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAg                                 <  2:312806/99‑1 (MQ=255)
                                                                   ctctTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAg                                 <  1:1520858/99‑1 (MQ=255)
                                                                   ctctTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAg                                 >  2:427693/1‑99 (MQ=255)
                                                                    tctTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGt                                <  2:676224/99‑1 (MQ=255)
                                                                    tctTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGt                                >  2:2243804/1‑99 (MQ=255)
                                                                    tctTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGt                                >  2:2311645/1‑99 (MQ=255)
                                                                    tctTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGt                                >  2:960801/1‑99 (MQ=255)
                                                                    tctTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGt                                >  1:849453/1‑99 (MQ=255)
                                                                     ctTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGtt                               <  1:1066080/99‑1 (MQ=255)
                                                                      ttCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTg                              >  1:1985409/1‑99 (MQ=255)
                                                                        cGGCAGCCTCGTTGGAGTAGAGAGAGCTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCa                            <  1:87679/99‑1 (MQ=255)
                                                                         ggCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAg                           <  2:1234970/99‑1 (MQ=255)
                                                                          gCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGt                          <  2:851147/99‑1 (MQ=255)
                                                                           cAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGtt                         <  1:524004/99‑1 (MQ=255)
                                                                             gCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTg                       >  2:1932073/1‑99 (MQ=255)
                                                                                tCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCg                    >  2:480667/1‑99 (MQ=255)
                                                                                  gTTGGAGTAGAGAGAGCTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGAGCCGTCCTCCAGTTGCAGTTTGTCGGt                  <  1:2256722/99‑1 (MQ=255)
                                                                                    tGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGtca                <  2:2331602/99‑1 (MQ=255)
                                                                                    tGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGtca                <  2:1374693/99‑1 (MQ=255)
                                                                                      gAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGtcatc              >  1:1134606/1‑99 (MQ=255)
                                                                                      gAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGtcatc              <  1:899666/99‑1 (MQ=255)
                                                                                      gAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGtcatc              >  2:2239399/1‑99 (MQ=255)
                                                                                          agagagagATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACa          <  1:1967486/99‑1 (MQ=255)
                                                                                          agagagagATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACa          <  2:457389/99‑1 (MQ=255)
                                                                                            agagagaTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACATa        <  1:236369/99‑1 (MQ=255)
                                                                                            agagagaTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACATa        >  1:1259331/1‑99 (MQ=255)
                                                                                            agagagaTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACATa        <  1:1391190/99‑1 (MQ=255)
                                                                                            agagagaTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACATa        <  1:548373/99‑1 (MQ=255)
                                                                                              agagaTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACATACg      >  1:2077497/1‑99 (MQ=255)
                                                                                              agagaTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACATACg      <  1:1836815/99‑1 (MQ=255)
                                                                                               gagaTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACATACGc     <  2:2254433/99‑1 (MQ=255)
                                                                                                 gaTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACATACGCAc   >  2:230330/1‑99 (MQ=255)
                                                                                                  aTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACATACGCACt  >  2:846882/1‑99 (MQ=255)
                                                                                                  aTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACATACGCACt  <  1:751413/99‑1 (MQ=255)
                                                                                                  aTCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACATACGCACt  >  2:1842734/1‑99 (MQ=255)
                                                                                                  |                                                                                                  
GGCCGCCCAGTTGCTGCGCCATCTGATGGCTTAAGAATTCATAGTGTTGCGCCCAACGCAGACCTGCCTCTTCGGCAGCCTCGTTGGAGTAGAGAGAGATCCCCAGCTCACGGTAGATCTCCCGACGGAAAGCGGGAGGCACAACCAACTGCGTGCCGTCCTCCAGTTGCAGTTTGTCGGTCATCAACATACGCACT  >  NC_000913/3535514‑3535710

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