Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,549,769 Δ1 bp 5.6% coding (302/2085 nt) malQ ← 4‑alpha‑glucanotransferase (amylomaltase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,549,7690G. 5.6% 227.7 / 7.6 72coding (302/2085 nt)malQ4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base . (2/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

AACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGCGT  >  NC_000913/3549627‑3549917
                                                                                                                                              |                                                                                                                                                    
aaCGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGAt                                                                                                                                                >  2:504629/1‑149 (MQ=255)
 aagcacgcacACCACAGCTTTTTTTTCTTCAGCTACGCCTTCGGTTCGTAACAGCGTTTCGGGGCTTCAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGAtt                                                                                                                                               <  1:127864/147‑1 (MQ=255)
 aCGCAGGCACCCCACAGCTTTTGTTTATTCTGGAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGAtt                                                                                                                                               <  1:925350/149‑1 (MQ=255)
 aCGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCCCCCGGCAATGCGCGCGCTGGTCGTCCTTGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGAtt                                                                                                                                               >  2:971030/1‑149 (MQ=255)
    cAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGaa                                                                                                                                            <  1:668725/149‑1 (MQ=255)
    cAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGGGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGaa                                                                                                                                            <  2:1433299/149‑1 (MQ=255)
         aCCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCtt                                                                                                                                       >  2:1461806/1‑149 (MQ=255)
           cccACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTg                                                                                                                                     >  1:819458/1‑149 (MQ=255)
           cccACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTg                                                                                                                                     >  1:604265/1‑149 (MQ=255)
            ccACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGc                                                                                                                                    <  2:991126/149‑1 (MQ=255)
             cacaGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGcc                                                                                                                                   <  1:1247047/149‑1 (MQ=255)
              acaGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGccc                                                                                                                                  >  2:420572/1‑149 (MQ=255)
                 gCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGcccccc                                                                                                                               >  1:558583/1‑149 (MQ=255)
                        tttATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACAt                                                                                                                        >  2:1430091/1‑149 (MQ=255)
                        tttATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACAt                                                                                                                        >  2:56481/1‑149 (MQ=255)
                          tATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATgg                                                                                                                      <  1:359000/149‑1 (MQ=255)
                           aTTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGc                                                                                                                     >  1:1526670/1‑149 (MQ=255)
                                gCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATgg                                                                                                                      <  1:552601/143‑1 (MQ=255)
                                gCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATgg                                                                                                                      >  2:552601/1‑143 (MQ=255)
                                gCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGt                                                                                                               >  2:1562049/1‑149 (MQ=255)
                                gCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGt                                                                                                               >  2:788991/1‑149 (MQ=255)
                                 cAACGCCTGCGGTTCGTAACAGCGTTTCGGGTCGACAAACAACCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGAGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTATGCCCCCCCTTACATGGccagtgt                                                                                                               >  1:1238644/1‑144 (MQ=255)
                                 cAACGCCGGCGGTTCGTAACCGCGCTTCGGGGCGACAATCACCCGGCAATGCGCGCGCGGGTCGGCCTGGGGGAGTGTCAGCGTGGGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGt                                                                                                               <  2:114118/149‑1 (MQ=255)
                                  aaCGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTa                                                                                                              >  2:30787/1‑149 (MQ=255)
                                  aaCGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTa                                                                                                              <  2:474946/149‑1 (MQ=255)
                                   aCGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGCACGCTTTGCCCCCCGTTACATGGCCTTTGTAc                                                                                                             >  1:875533/1‑149 (MQ=255)
                                   aCGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTAc                                                                                                             <  2:953522/149‑1 (MQ=255)
                                    cGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCGTTGCCCCCCGTTACATGGCCTTTGTACt                                                                                                            >  1:856451/1‑149 (MQ=255)
                                       cTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCg                                                                                                         <  1:542628/149‑1 (MQ=255)
                                        tGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGt                                                                                                        >  2:1140331/1‑149 (MQ=255)
                                           ggTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGCGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTcc                                                                                                     <  2:395930/149‑1 (MQ=255)
                                            gTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCct                                                                                                    <  1:504629/149‑1 (MQ=255)
                                               cGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCcttct                                                                                                 <  2:216999/149‑1 (MQ=255)
                                                gTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCcttctt                                                                                                >  1:814303/1‑149 (MQ=255)
                                                 tAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCttcttc                                                                                               >  2:73212/1‑149 (MQ=255)
                                                 tAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCg                                                                                              <  1:788991/149‑1 (MQ=255)
                                                 tAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGGCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCttcttc                                                                                               <  1:810534/149‑1 (MQ=255)
                                                 tAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGACGTCCTGGGTGACTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGCAAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCttattc                                                                                               >  1:939561/1‑149 (MQ=255)
                                                 gaaCAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCg                                                                                              <  1:1562049/148‑1 (MQ=255)
                                                   aCAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCGTCTTCgg                                                                                             >  2:559574/1‑149 (MQ=255)
                                                    cAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCGGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAAACTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGTCCTTTGTACTGCGTTCCTTCTTCggt                                                                                            <  1:62760/149‑1 (MQ=255)
                                                     tgCGTTTCGGGGCGACAATCACCCGTCACTGCGCGCGCTGGTCGTGCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGGCGGTAGATTGTACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtg                                                                                           <  1:304774/148‑1 (MQ=255)
                                                     aGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtg                                                                                           <  2:536379/149‑1 (MQ=255)
                                                         tttCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTca                                                                                       >  1:1321707/1‑149 (MQ=255)
                                                         tttCGGGGCGACAATCACACGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTca                                                                                       >  2:1179761/1‑149 (MQ=255)
                                                             ggggCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCagcagc                                                                                  >  2:1312795/1‑149 (MQ=255)
                                                              gggCGACAATCACCCGGCAATGCGCGCGCGGGTCGTCCTGGGGGGGTGTCAGCGGGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGCACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCagcagc                                                                                  <  1:971030/149‑1 (MQ=255)
                                                                     aaacACCCGGCAATGGGCGCGCTGGTCGTCCAGGGTGAGTGTCAGCGTGTGAAAACCATACGGCAGCTTCGGCGGTAGAATGAACGCTTTGGCCCCCGGTTCATTGACTTTGTACTGCGTTGCTTCATCGGTGGTCAGCAGCCAACtat                                                                           >  2:785319/4‑149 (MQ=255)
                                                                                aaTGCGCGCGCTGGTCGGCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTg                                                                <  2:1209310/149‑1 (MQ=255)
                                                                                       gcgcTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCTCTGCCCTcca                                                         <  1:919778/149‑1 (MQ=255)
                                                                                        cgcTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTccac                                                        <  2:604265/149‑1 (MQ=255)
                                                                                             gtcgtcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTccaccacca                                                   >  2:302668/1‑149 (MQ=255)
                                                                                             gtcgtcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTccaccacca                                                   >  1:489140/1‑149 (MQ=255)
                                                                                                       gTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCACTCCACCACCATCGGCAtttt                                         >  2:981979/1‑149 (MQ=255)
                                                                                                         gAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCAtttttt                                       >  2:813620/1‑149 (MQ=255)
                                                                                                          agggtCAGCGTGTGATTACCTTCCGGCCGGTGCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCAttttttt                                      <  2:91009/146‑1 (MQ=255)
                                                                                                          aGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCAttttttt                                      >  2:429020/1‑149 (MQ=255)
                                                                                                            tgtCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGc                                    >  1:747899/1‑149 (MQ=255)
                                                                                                            tgtCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGGCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGc                                    <  1:1188897/149‑1 (MQ=255)
                                                                                                              tCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCg                                  >  1:1157021/1‑149 (MQ=255)
                                                                                                               cAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGc                                 <  2:419795/149‑1 (MQ=255)
                                                                                                                aGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCt                                >  2:788518/1‑149 (MQ=255)
                                                                                                                  cGTGTGATTAACTTCCGGCAGCTTCGGCGGTTGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTGCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTgg                              <  2:875533/149‑1 (MQ=255)
                                                                                                                   gtgtGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGt                             >  2:150943/1‑149 (MQ=255)
                                                                                                                    ggtgATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGAGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTgtat                            <  1:195692/148‑4 (MQ=255)
                                                                                                                             cTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATg                   >  1:861538/1‑149 (MQ=255)
                                                                                                                              ttCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGa                  >  1:1317336/1‑149 (MQ=255)
                                                                                                                                   gctgCTTGGTCGGTTGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCTGAGAGCTATATTCGCCGCTTCCCTCCACCACCATCTGCATTTTTTTGCCGCTGGTATAAACCATGACATTc             <  2:939561/146‑1 (MQ=255)
                                                                                                                                    cAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCg            <  1:73212/149‑1 (MQ=255)
                                                                                                                                     aGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCgg           <  2:1321707/149‑1 (MQ=255)
                                                                                                                                       tttCGTCGGTAGATTGAACGCTTTGCGCCTCGTTACTTGGGCTTTGTACTGGGTTCCTTCTTCTGTGGTCAGCAGCGAGATATTTTCGCCGCTGGCCTCCATGACCATCGGCATTTTTTTGCGGCTGGTATAATCGATGACATTCggtt         <  2:1238644/148‑3 (MQ=255)
                                                                                                                                        ttCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGAc        >  1:527415/1‑149 (MQ=255)
                                                                                                                                         tcgtcgGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGACGCTGGTATAAACCATGACATTCGGGACt       >  2:312068/1‑149 (MQ=255)
                                                                                                                                          cgtcgGTAGATTGAACTCTTTGCCCCCCGTTACATTGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTg      <  1:56481/149‑1 (MQ=255)
                                                                                                                                          cgtcgGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTg      <  2:819458/149‑1 (MQ=255)
                                                                                                                                          cgtcgGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTg      <  1:948876/149‑1 (MQ=255)
                                                                                                                                          cgtcgGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTATGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCGGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTg      <  1:1430091/149‑1 (MQ=255)
                                                                                                                                            tcgTTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGc    >  1:760543/1‑149 (MQ=255)
                                                                                                                                            tcgGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGc    <  2:665600/149‑1 (MQ=255)
                                                                                                                                            tcgGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGc    <  1:302668/149‑1 (MQ=255)
                                                                                                                                              ggTAGATTGAACGCTTTGCCACCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCAACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGCGt  >  1:1104928/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
AACGCAGGCACCCCACAGCTTTTGTTTATTCAGCAACGCCTGCGGTTCGTAACAGCGTTTCGGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTGCCGCTGGTATAAACCATGACATTCGGGACTGGCGT  >  NC_000913/3549627‑3549917

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