Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,713,362 A→T 5.1% K91* (AAA→TAA)  tag → 3‑methyl‑adenine DNA glycosylase I, constitutive

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,713,3620AT 5.1% 258.1 / 3.0 79K91* (AAA→TAA) tag3‑methyl‑adenine DNA glycosylase I, constitutive
Reads supporting (aligned to +/- strand):  ref base A (37/38);  new base T (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

AGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTC  >  NC_000913/3713213‑3713510
                                                                                                                                                     |                                                                                                                                                    
agGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGa                                                                                                                                                      >  1:1656858/1‑150 (MQ=255)
 ggCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGaa                                                                                                                                                     >  2:1318145/1‑150 (MQ=255)
 ggCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGACGGGATTATCCGCCATAGAGGGaa                                                                                                                                                     >  1:449221/1‑150 (MQ=255)
  gCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGaaa                                                                                                                                                    <  1:408542/150‑1 (MQ=255)
  gCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGaaa                                                                                                                                                    <  1:198805/150‑1 (MQ=255)
    tGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAAt                                                                                                                                                  <  1:1031289/150‑1 (MQ=255)
         tATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAgg                                                                                                                                             <  1:831095/150‑1 (MQ=255)
         tATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAgg                                                                                                                                             <  1:1442386/150‑1 (MQ=255)
           tCGTGGATCAGCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCa                                                                                                                                           <  1:180881/150‑1 (MQ=255)
            cGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCaa                                                                                                                                          >  1:809036/1‑150 (MQ=255)
            cGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCaa                                                                                                                                          >  2:721566/1‑150 (MQ=255)
              tGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGTAAATTCAGGCAatt                                                                                                                                        >  2:161344/1‑150 (MQ=255)
              tGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGTAAATTCAGGCAatt                                                                                                                                        >  2:155374/1‑150 (MQ=255)
               ggATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAatta                                                                                                                                       >  1:268671/1‑150 (MQ=255)
                gATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAattat                                                                                                                                      >  1:875063/1‑150 (MQ=255)
                    aCCGTCATGAAAATACGCGATAACTATCGCGCCTGCTTTCTTGAGTTCGATCCGGTGAAGGTCGATGCATTGCAGTAAGAGTATGTCGAAAGACTGTTACAGGACGCCGGGATTATCCGCGATCGAGGGAAAATTCAGTCAATTATTGGt                                                                                                                                  <  2:449221/150‑1 (MQ=255)
                             aaaaaaCGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCgg                                                                                                                         <  1:398377/150‑1 (MQ=255)
                             aaaaaaCGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCgg                                                                                                                         >  1:756242/1‑150 (MQ=255)
                              aaaaaCGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCggg                                                                                                                        >  2:710016/1‑150 (MQ=255)
                              aaaaaCGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCggg                                                                                                                        >  2:713391/1‑150 (MQ=255)
                               aaaaCGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGc                                                                                                                       <  2:696926/150‑1 (MQ=255)
                                aaaCGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCg                                                                                                                      <  2:801002/150‑1 (MQ=255)
                                       aaaaCTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTAACCGCCATCGAGGGAAAATTCAGCCAATTATTGGTAATGCGCGGGCGTACCTGc                                                                                                               >  1:1671096/1‑150 (MQ=255)
                                              tCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGa                                                                                                        <  1:1099315/150‑1 (MQ=255)
                                                gcgcCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAAc                                                                                                      <  1:893666/150‑1 (MQ=255)
                                                 cgcCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACa                                                                                                     <  1:759505/150‑1 (MQ=255)
                                                       cTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCATACCTGCAAATGGAACAGAACgg                                                                                               >  1:1512803/1‑150 (MQ=255)
                                                         ttCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCg                                                                                             >  2:610018/1‑150 (MQ=255)
                                                          tcatcaGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGa                                                                                            >  2:1189562/1‑150 (MQ=255)
                                                            atcaGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAAc                                                                                          >  1:457762/1‑150 (MQ=255)
                                                              caGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCg                                                                                        >  2:415033/1‑150 (MQ=255)
                                                               aGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACGGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGt                                                                                       >  1:1417632/1‑150 (MQ=255)
                                                                       ccGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGAc                                                                               <  1:333088/150‑1 (MQ=255)
                                                                           tGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACGTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTg                                                                           <  2:268671/150‑1 (MQ=255)
                                                                           tGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTg                                                                           >  2:1376788/1‑150 (MQ=255)
                                                                             aaGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACTGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTc                                                                         <  2:1417632/150‑1 (MQ=255)
                                                                             aaGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTc                                                                         <  1:710016/150‑1 (MQ=255)
                                                                             aaGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTc                                                                         <  1:713391/150‑1 (MQ=255)
                                                                              aGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCt                                                                        >  1:951102/1‑150 (MQ=255)
                                                                               ggTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTg                                                                       <  1:415033/150‑1 (MQ=255)
                                                                                 gcGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGt                                                                     <  2:504752/149‑1 (MQ=255)
                                                                                    cagcaATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGt                                                                  >  1:1529288/1‑150 (MQ=255)
                                                                                       caATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTg                                                               <  2:1220001/150‑1 (MQ=255)
                                                                                         aTGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTa                                                             <  1:13783/150‑1 (MQ=255)
                                                                                             aGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACAGTTTGTCGACTTTCTCTGGTCGTTTGTAAatc                                                         >  1:1264974/1‑150 (MQ=255)
                                                                                              ggAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGTCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAatca                                                        <  2:1321181/150‑1 (MQ=255)
                                                                                              ggAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAatca                                                        >  1:208359/1‑150 (MQ=255)
                                                                                                aaGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAatcata                                                      >  1:291539/1‑149 (MQ=255)
                                                                                                    ggATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGcc                                                  <  2:809036/150‑1 (MQ=255)
                                                                                                      aTGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCcac                                                <  2:1656858/150‑1 (MQ=255)
                                                                                                        gTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAg                                              >  1:14555/1‑150 (MQ=255)
                                                                                                         tCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAgg                                             <  1:121561/150‑1 (MQ=255)
                                                                                                          cGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGt                                            >  2:441275/1‑150 (MQ=255)
                                                                                                           gAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTg                                           <  1:1318145/150‑1 (MQ=255)
                                                                                                           gAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTg                                           <  2:250938/150‑1 (MQ=255)
                                                                                                           gAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTAACTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTAGTCTGGTCGTTTGTAAATCATCAGCCACAGGTg                                           >  2:1505248/1‑150 (MQ=255)
                                                                                                               gACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGacac                                       >  1:817935/1‑150 (MQ=255)
                                                                                                                 cTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACaa                                     <  1:721566/150‑1 (MQ=255)
                                                                                                                 cTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACaa                                     >  1:594961/1‑150 (MQ=255)
                                                                                                                  tGGTACAGGACGCCGGGATTATCCGCCATCGAGGGTAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAg                                    <  1:155374/150‑1 (MQ=255)
                                                                                                                  tGGTACAGGACGCCGGGATTATCCGCCATCGAGGGTAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAg                                    <  1:161344/150‑1 (MQ=255)
                                                                                                                  tGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAg                                    <  2:756242/150‑1 (MQ=255)
                                                                                                                  tGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAg                                    >  1:987694/1‑150 (MQ=255)
                                                                                                                  tGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCATACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAg                                    <  2:1512803/150‑1 (MQ=255)
                                                                                                                   ggTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGc                                   >  1:1218596/1‑150 (MQ=255)
                                                                                                                   ggTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGc                                   >  1:837284/1‑150 (MQ=255)
                                                                                                                   ggTACAGGACGCCGGGATTATCAGCCATCGAGGGAAAATTCAGTCAATTATTGGTAATGCGCGGGAGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGc                                   >  2:1088082/1‑150 (MQ=255)
                                                                                                                   ggTACACGACGACTAGATTATCCGACATCCACGGAAAATTAAGGCAATTATTGCTAATCATATGTCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGc                                   >  2:750030/1‑150 (MQ=255)
                                                                                                                      aCAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCaac                                >  1:75737/1‑150 (MQ=255)
                                                                                                                         ggACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAAc                             <  1:979953/150‑1 (MQ=255)
                                                                                                                           aCGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGt                           >  1:1115724/1‑150 (MQ=255)
                                                                                                                            cGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGtt                          >  2:412062/1‑150 (MQ=255)
                                                                                                                                gggATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGc                      <  1:339751/150‑1 (MQ=255)
                                                                                                                                     tATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAAt                 <  1:441275/150‑1 (MQ=255)
                                                                                                                                         ctccATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGTTCGTTTGTAATTCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTccc             <  2:14137/148‑1 (MQ=255)
                                                                                                                                          gCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCca            >  1:836880/1‑150 (MQ=255)
                                                                                                                                           ccATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGGGACACAAGCCACAACGTTGAGCGAAATTCCcac           <  2:20074/150‑1 (MQ=255)
                                                                                                                                               cGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCt       <  2:14555/150‑1 (MQ=255)
                                                                                                                                               cGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCt       <  1:1189562/150‑1 (MQ=255)
                                                                                                                                                gAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTa      <  2:594961/150‑1 (MQ=255)
                                                                                                                                                  gggAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACg    >  2:343357/1‑150 (MQ=255)
                                                                                                                                                   ggAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGt   >  2:526211/1‑150 (MQ=255)
                                                                                                                                                   ggAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGAGAAATTCCCACATCTACGt   >  1:1290699/1‑150 (MQ=255)
                                                                                                                                                    gAAAATTCAGTCAATTATTGGTAATGCGCGGGCGTACCTGCATATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTTACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTc  <  2:1671096/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGGCTGGATTATCGTGGATCACCGTCCTCAAAAAACGCGAAAACTATCGCGCCTGCTTTCATCAGTTCGATCCGGTGAAGGTCGCAGCAATGCAGGAAGAGGATGTCGAAAGACTGGTACAGGACGCCGGGATTATCCGCCATCGAGGGAAAATTCAGGCAATTATTGGTAATGCGCGGGCGTACCTGCAAATGGAACAGAACGGCGAACCGTTTGTCGACTTTGTCTGGTCGTTTGTAAATCATCAGCCACAGGTGACACAAGCCACAACGTTGAGCGAAATTCCCACATCTACGTC  >  NC_000913/3713213‑3713510

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: ‑