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,5410CT98.0% 325.9 / ‑2.8 101R159* (CGA→TGA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (51/47);  total (53/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 107 total aligned reads displayed.

GTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAG  >  NC_000913/236445‑236637
                                                                                                |                                                                                                
gTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaa                                                                                               >  1:110602/1‑100 (MQ=255)
gTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaa                                                                                               >  1:1724330/1‑100 (MQ=255)
gTCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaa                                                                                               >  2:1780343/1‑100 (MQ=255)
 tCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaa                                                                                               <  1:61070/99‑1 (MQ=255)
 tCACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaacgaac                                                                                              >  2:1263777/1‑100 (MQ=255)
  cACCGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAaatgaacg                                                                                             >  2:28963/1‑100 (MQ=255)
    ccGATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCt                                                                                           >  1:200579/1‑100 (MQ=255)
      gATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGc                                                                                         >  2:520180/1‑100 (MQ=255)
      gATAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGc                                                                                         >  2:4575/1‑100 (MQ=255)
       aTAGCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCa                                                                                        <  2:1727808/100‑1 (MQ=255)
          gCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACgg                                                                                     >  1:2339064/1‑100 (MQ=255)
          gCCTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCCCGCTACGAAATAGATAACAGTAAACGAACGCTGCACgg                                                                                     >  1:1675667/1‑100 (MQ=255)
           ccTTGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGACATAGATAACAGTAAATGAACGCTGCACggg                                                                                    >  2:837635/1‑100 (MQ=255)
             ttGCGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGc                                                                                  <  1:712591/100‑1 (MQ=255)
                cGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCAtt                                                                               <  2:431483/100‑1 (MQ=255)
                cGGTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCAtt                                                                               <  1:28964/100‑1 (MQ=255)
                  gTGGCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTAc                                                                             >  2:160140/1‑100 (MQ=255)
                     gCGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATTCGTTATTTTCTCTCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCg                                                                          <  1:1180570/100‑1 (MQ=255)
                      cGAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGa                                                                         >  1:231268/1‑100 (MQ=255)
                       gAGGAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGAt                                                                        <  2:1219612/100‑1 (MQ=255)
                         ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                      <  2:1248075/100‑1 (MQ=255)
                         ggAAAATGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGc                                                                      >  2:1051562/1‑100 (MQ=255)
                           aaaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGccc                                                                    <  2:2343801/100‑1 (MQ=255)
                           aaaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGccc                                                                    >  1:2118886/1‑100 (MQ=255)
                           aaaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACAGGGCATTACTCGATGccc                                                                    >  2:912341/1‑100 (MQ=255)
                            aaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCa                                                                   >  2:1313170/1‑100 (MQ=255)
                            aaaTGTTTCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCa                                                                   >  2:1290163/1‑100 (MQ=255)
                              aTGTTTCCCAGTAAGAGAAACAGCCTGGATGCGTGTTGTGCTCGCTTCGAAATAGATAACAGTATATGAACGCTGCACGGGGCATTACTCGATGTCCAg                                                                  <  2:523209/99‑1 (MQ=255)
                                  ttCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCTCAGATcc                                                              >  2:448977/1‑99 (MQ=255)
                                  ttCCCGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCt                                                             <  1:2253350/100‑1 (MQ=255)
                                    cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTg                                                           >  2:1321699/1‑100 (MQ=255)
                                    cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTg                                                           <  2:1810222/100‑1 (MQ=255)
                                    cccGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTg                                                           <  1:746621/100‑1 (MQ=255)
                                      cGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCg                                                         <  2:1412875/100‑1 (MQ=255)
                                      cGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCg                                                         >  1:1071853/1‑100 (MQ=255)
                                      cGGTAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCg                                                         >  1:1209246/1‑100 (MQ=255)
                                       ggTAAGCGCAACAGCCTCAATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCg                                                         >  2:2170773/1‑99 (MQ=255)
                                         tAAGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGaa                                                      <  2:249031/100‑1 (MQ=255)
                                          aaGCGCAACATCCTTGATGCGTTATGTGCTCGCTACGAAATAGATAACATTAAATCAACGCTGTACTGGGCATTACTCGATGCCCAGATCCTCGAGGaa                                                      <  2:2097671/99‑1 (MQ=255)
                                          aaGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAg                                                     <  1:1342747/100‑1 (MQ=255)
                                          aaGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAg                                                     <  1:1326880/100‑1 (MQ=255)
                                           aGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGt                                                    <  2:2006660/100‑1 (MQ=255)
                                           aGCGCAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGt                                                    <  1:1882649/100‑1 (MQ=255)
                                            gcgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCAGGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGt                                                    >  1:279305/1‑99 (MQ=255)
                                            gcgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGtt                                                   <  1:2196695/100‑1 (MQ=255)
                                            gcgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGtt                                                   <  2:2376974/100‑1 (MQ=255)
                                            gcgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGtt                                                   <  2:1881918/100‑1 (MQ=255)
                                             cgcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGACAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGt                                                    <  2:64584/98‑1 (MQ=255)
                                              gcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTa                                                 <  2:1963257/100‑1 (MQ=255)
                                              gcAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTa                                                 >  1:1821509/1‑100 (MQ=255)
                                               cAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTAt                                                >  1:1139016/1‑100 (MQ=255)
                                               cAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTAt                                                >  2:2012274/1‑100 (MQ=255)
                                               cAACAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAAAAGTAAATGAACGCTGCACGGGGAATTACTCGATGCCCAGATCCTTGCGGAAGTTTAt                                                >  2:1883841/1‑100 (MQ=255)
                                                 aCAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCt                                              >  1:756523/1‑100 (MQ=255)
                                                 aCAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCAAGGGGCATTACTCGATTCCCAGATCCTTGCGGAAGTTTATCt                                              >  2:241774/1‑100 (MQ=255)
                                                  cAGCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCt                                              >  1:488782/1‑99 (MQ=255)
                                                    gtctgggatcGTTATGTGCTCGCTACGAAATAGAGAACAGTAAATGGATGCTGCACGTGGAGTTACTCGATGGCCAGATCCTTGCGGAAGTTGATGTGGc                                           <  2:1883046/91‑1 (MQ=255)
                                                    gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTgg                                            <  2:1858580/99‑1 (MQ=255)
                                                    gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTg                                             >  1:1376816/1‑98 (MQ=255)
                                                    gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           >  2:1837164/1‑100 (MQ=255)
                                                    gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           <  2:2222022/100‑1 (MQ=255)
                                                    gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           >  1:359652/1‑100 (MQ=255)
                                                    gCCTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGc                                           >  2:583247/1‑100 (MQ=255)
                                                     ccTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAAAGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCg                                          >  2:1056194/1‑100 (MQ=255)
                                                      cTCGATGCGTTATGTGCTCGCTACGAAATAGATAGCAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCg                                          <  1:357388/99‑1 (MQ=255)
                                                      cTCGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGa                                         <  1:1268928/100‑1 (MQ=255)
                                                        cGATGCGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATg                                       >  1:1101962/1‑100 (MQ=255)
                                                             cGTTATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCgg                                  <  2:150984/100‑1 (MQ=255)
                                                                tATGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtgg                               <  1:1567381/100‑1 (MQ=255)
                                                                 aTGTGCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCggtggt                              <  2:1122027/100‑1 (MQ=255)
                                                                  tgtgCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTc                             >  1:1503060/1‑100 (MQ=255)
                                                                  tgtgCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTc                             >  1:1918946/1‑100 (MQ=255)
                                                                     gCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCaaa                          >  2:1417208/1‑100 (MQ=255)
                                                                     gCTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCaaa                          >  2:771472/1‑100 (MQ=255)
                                                                      cTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAAc                         <  1:1417223/100‑1 (MQ=255)
                                                                      cTCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAAc                         <  2:2366263/100‑1 (MQ=255)
                                                                       tCGCTACGAAATAGATAACAGTAAATGAACGCTGCACGTGGCATTACTCGAGGCCCAGATCCTTGCGGAAGTTTATCTGGTGATGACCGGTGGGCAAACg                        <  2:1134767/100‑1 (MQ=255)
                                                                        cGCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGt                       >  2:279642/1‑100 (MQ=255)
                                                                         gCTACGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTc                      <  2:353478/100‑1 (MQ=255)
                                                                          cTACGAAATAGATAACAGTAAACGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCg                     <  1:1263790/100‑1 (MQ=255)
                                                                             cGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                  <  2:1613953/100‑1 (MQ=255)
                                                                             cGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                  <  2:1590093/100‑1 (MQ=255)
                                                                             cGAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATg                  <  1:753610/100‑1 (MQ=255)
                                                                              gAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATgg                 >  2:1197585/1‑100 (MQ=255)
                                                                              gAAATAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATgg                 >  2:473028/1‑100 (MQ=255)
                                                                                  tAGATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCttt             <  1:2109346/100‑1 (MQ=255)
                                                                                    gATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCTCAGATCCTTGCGGAAGTTTATCTGGCGATGCAAGGTGGTCAAACATCGATGGCtttt            >  1:52186/1‑99 (MQ=255)
                                                                                    gATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCACACGTCGATGGCTTTTg           >  2:2119243/1‑100 (MQ=255)
                                                                                    gATAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTg           <  2:2118862/100‑1 (MQ=255)
                                                                                      tAACAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTATg           >  1:758174/1‑98 (MQ=255)
                                                                                        aCAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGAt       <  1:279645/100‑1 (MQ=255)
                                                                                         cAGTAACTGAACGCTACACGGGGAATTACTCCATGCCCAGATCCTTGCAGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTATGCGATg      >  2:1654366/1‑100 (MQ=255)
                                                                                         cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg      >  1:1004395/1‑100 (MQ=255)
                                                                                         cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATg      >  2:1526769/1‑100 (MQ=255)
                                                                                         cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGGCAAACGTCGATGGCTTTTGCGATg      >  1:1671009/1‑100 (MQ=255)
                                                                                         cAGTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATACTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGGtttt            >  2:1989046/1‑94 (MQ=255)
                                                                                           gTAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGa    >  1:2055826/1‑100 (MQ=255)
                                                                                            tAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGaa   >  2:1294604/1‑100 (MQ=255)
                                                                                            tAAATGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGaa   >  2:1296616/1‑100 (MQ=255)
                                                                                             aaaTGAACGCTGCACGGGGCATTACTCGATGCCCAGATCCTTGCGGAAGTTTATCTGGCGATGACCGGTGGTCAAACGTCGATGGCTTTTGCGATGGAAg  <  1:661853/100‑1 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: