Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 236,541 C→T 100% R159* (CGA→TGA)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913236,5410CT99.3% 493.7 / ‑4.8 147R159* (CGA→TGA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (84/62);  total (85/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 152 total aligned reads displayed.

GTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAG  >  NC_000913/236445‑236637
                                                                                                |                                                                                                
gTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaa                                                                                               >  2:250675/1‑100 (MQ=255)
gTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaa                                                                                               >  1:2372187/1‑100 (MQ=255)
  cACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaacg                                                                                             <  2:1823702/100‑1 (MQ=255)
  cACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAAAAGCCTCGATGGGTTATGTGCTGGCTACGAAATAGATAACAGTAaatgaacg                                                                                             >  1:3909119/1‑100 (MQ=255)
   agcGATTGCGTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGGAACAGCCTCGATGCGTTATGTGCTCTCTACGAAATAGATAACAGTAAATGAACGc                                                                                            <  1:63497/98‑1 (MQ=255)
      gATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGc                                                                                         >  1:243704/1‑100 (MQ=255)
      gATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGc                                                                                         <  1:32296/100‑1 (MQ=255)
          gCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACgg                                                                                     >  1:1266969/1‑100 (MQ=255)
          gCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACgg                                                                                     >  1:2917201/1‑100 (MQ=255)
          gCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACgg                                                                                     >  2:1213926/1‑100 (MQ=255)
           cctaccGGTAGAGAGCAAACTGTTCACCGCTAAGCGCAACCGACTCGACGAGCTCCGCGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCat                                                                                       <  1:352186/92‑2 (MQ=255)
           ccTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACggg                                                                                    >  1:1961469/1‑100 (MQ=255)
            cTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACgggg                                                                                   <  1:1514961/100‑1 (MQ=255)
              tGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGc                                                                                  <  2:732528/99‑1 (MQ=255)
                 ggTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTa                                                                              >  1:258316/1‑100 (MQ=255)
                 ggTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTa                                                                              <  1:2746690/100‑1 (MQ=255)
                 ggTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACACTGCACGGGGCATTa                                                                              >  2:4217025/1‑100 (MQ=255)
                      cGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGa                                                                         >  1:3164382/1‑100 (MQ=255)
                         ngAAAATGTTTCCCGGTAAGCGCAACAGCCNCNANGCGNTNTGNGCNCGCTNCGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                      >  2:1970139/2‑100 (MQ=255)
                         ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                      >  1:1512963/1‑100 (MQ=255)
                         ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                      >  2:722027/1‑100 (MQ=255)
                         ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                      <  1:1019269/100‑1 (MQ=255)
                         ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                      <  1:1122639/100‑1 (MQ=255)
                         ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                      <  1:1019397/100‑1 (MQ=255)
                          gAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGcc                                                                     <  2:1091735/100‑1 (MQ=255)
                           aaaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGccc                                                                    >  2:293602/1‑100 (MQ=255)
                               tGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGAt                                                                <  2:2198096/100‑1 (MQ=255)
                                gTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATc                                                               <  2:4018750/100‑1 (MQ=255)
                                 tttCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATcc                                                              >  2:1721701/1‑100 (MQ=255)
                                  ttCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCt                                                             >  2:826778/1‑100 (MQ=255)
                                  ttCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCt                                                             >  2:2241663/1‑100 (MQ=255)
                                  ttCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCt                                                             >  2:3726419/1‑100 (MQ=255)
                                  ttCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCt                                                             >  1:2423256/1‑100 (MQ=255)
                                    cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTGGCTACGAAGTAGATAACAGTAAATGAACGCTGCACGGGGCATGACTCGAGGCGCAGATCGTTg                                                           >  1:787546/1‑100 (MQ=255)
                                    cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTg                                                           >  2:1763700/1‑100 (MQ=255)
                                    cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTg                                                           >  1:2182460/1‑100 (MQ=255)
                                    cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTg                                                           >  2:2455429/1‑100 (MQ=255)
                                     ccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTCCTCGATGCCCAGATCCTTGc                                                          <  1:2603479/100‑1 (MQ=255)
                                     ccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGc                                                          >  1:2056854/1‑100 (MQ=255)
                                     ccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGc                                                          <  2:176399/100‑1 (MQ=255)
                                       ggTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCgg                                                        <  2:1961473/100‑1 (MQ=255)
                                       ggTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCgg                                                        >  2:3693741/1‑100 (MQ=255)
                                          aaGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGTCATTACTCGATGCCCAGATCCTTGCgg                                                        >  1:2075954/1‑97 (MQ=255)
                                          aaGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAg                                                     <  2:1815956/100‑1 (MQ=255)
                                          aaGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAg                                                     <  1:1453231/100‑1 (MQ=255)
                                           aGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGt                                                    >  1:3499707/1‑100 (MQ=255)
                                            gcgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGtt                                                   >  2:3970993/1‑100 (MQ=255)
                                            gcgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGtt                                                   <  1:2209678/100‑1 (MQ=255)
                                             cgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGCTCCTTGCGGAAttt                                                   <  1:1467086/99‑1 (MQ=255)
                                             cgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGAGGCATTACTAGATGCCCAGATCCTTGCGGAAGttt                                                  >  1:2119463/1‑100 (MQ=255)
                                                 aCAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCt                                              <  2:2454258/100‑1 (MQ=255)
                                                   aGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTgg                                            >  1:3731283/1‑100 (MQ=255)
                                                    gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           >  2:348447/1‑100 (MQ=255)
                                                    gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           >  2:1992200/1‑100 (MQ=255)
                                                    gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           <  2:2406600/100‑1 (MQ=255)
                                                     ccTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCg                                          >  2:3285649/1‑100 (MQ=255)
                                                      cTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGa                                         >  2:2086351/1‑100 (MQ=255)
                                                      cTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGa                                         <  1:1771268/100‑1 (MQ=255)
                                                        cGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATg                                       <  2:3891439/100‑1 (MQ=255)
                                                          aTGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGAc                                     <  2:158480/100‑1 (MQ=255)
                                                            gCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCg                                   >  1:2145313/1‑100 (MQ=255)
                                                               ttATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtg                                >  1:1359505/1‑100 (MQ=255)
                                                                tATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACGCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtgg                               <  1:81633/100‑1 (MQ=255)
                                                                 aTGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtggt                              <  2:275958/100‑1 (MQ=255)
                                                                  tgtgCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTc                             <  2:1129735/100‑1 (MQ=255)
                                                                      cTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAAc                         <  1:1763387/100‑1 (MQ=255)
                                                                        cGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGt                       >  1:404660/1‑100 (MQ=255)
                                                                         gCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGCAGTTTATCTGGCGATGACCGGTGGTCAAACGTc                      >  1:2699042/1‑100 (MQ=255)
                                                                         gCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTc                      <  1:3301739/100‑1 (MQ=255)
                                                                         gCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACCCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTc                      <  1:1337601/100‑1 (MQ=255)
                                                                          cTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCg                     >  1:181307/1‑100 (MQ=255)
                                                                           tACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGa                    >  2:3953092/1‑100 (MQ=255)
                                                                            aCGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGAt                   >  2:2047077/1‑100 (MQ=255)
                                                                             cGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                  >  1:1139747/1‑100 (MQ=255)
                                                                             cGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                  >  1:1959337/1‑100 (MQ=255)
                                                                             cGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                  >  2:1086944/1‑100 (MQ=255)
                                                                             cGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                  >  1:677710/1‑100 (MQ=255)
                                                                                 aTAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCtt              <  1:324566/100‑1 (MQ=255)
                                                                                  tAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCttt             >  1:125993/1‑100 (MQ=255)
                                                                                   aGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCtttt            >  2:2610907/1‑100 (MQ=255)
                                                                                    gATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTg           >  2:583974/1‑100 (MQ=255)
                                                                                     aTAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGc          <  1:2553824/100‑1 (MQ=255)
                                                                                     aTAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGc          <  1:3294390/100‑1 (MQ=255)
                                                                                     aTAACAGTAAATGAACGCTGCACGGGGCAGTACTCGATGCCCAGATCCTTGCGGAAGGTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGc          <  1:507005/100‑1 (MQ=255)
                                                                                      tAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTGAAACGTCGATGGCTTTTGCg         <  1:2615161/100‑1 (MQ=255)
                                                                                      tAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCg         <  1:702828/100‑1 (MQ=255)
                                                                                      tAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCg         >  1:3318577/1‑100 (MQ=255)
                                                                                      tAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCg         >  2:2531433/1‑100 (MQ=255)
                                                                                      tAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCg         >  2:59913/1‑100 (MQ=255)
                                                                                         cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg      >  2:3484020/1‑100 (MQ=255)
                                                                                         cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg      >  2:3090906/1‑100 (MQ=255)
                                                                                         cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg      <  2:4265113/100‑1 (MQ=255)
                                                                                         cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg      <  2:215090/100‑1 (MQ=255)
                                                                                         cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTGCGATgg     >  1:180859/1‑100 (MQ=255)
                                                                                          aGTAAATGAACGCTGCACGGGGCATTTCTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATgg     >  1:3752681/1‑100 (MQ=255)
                                                                                          aGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATgg     <  1:3403898/100‑1 (MQ=255)
                                                                                          aGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTCCCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATgg     >  2:1084829/1‑100 (MQ=255)
                                                                                           gTAAATGACCGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGa    <  2:3272951/100‑1 (MQ=255)
                                                                                            tAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGaa   <  2:1200792/100‑1 (MQ=255)
                                                                                             aaaTGAACGCTGCACGGGGCATTACTCGATGCCCAGATACTTGCGGCAGTTTATCTGGCGATGACCGGCGGTCAAACGTCGATGGCTTTTGCGATGGAAg  >  1:459983/1‑100 (MQ=255)
                                                                                                |                                                                                                
GTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAG  >  NC_000913/236445‑236637

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