Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,549,792 (C)6→5 coding (279/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,7870C.97.5% 356.2 / 2.2 83coding (284/2085 nt)malQ4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base . (50/29);  total (52/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTG  >  NC_000913/3549688‑3549882
                                                                                                   |                                                                                               
ggggCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGc                                                                                                 >  1:1539443/1‑100 (MQ=255)
ggggCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGc                                                                                                 >  2:1442594/1‑100 (MQ=255)
 gggCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGcc                                                                                                <  2:2417309/100‑1 (MQ=255)
 gggCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGcc                                                                                                >  2:2115572/1‑100 (MQ=255)
  ggCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGccc                                                                                               <  2:2635568/100‑1 (MQ=255)
   gCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGcccc                                                                                              <  1:1089024/100‑1 (MQ=255)
       cAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGcccccgtt                                                                                          >  1:1198948/1‑97 (MQ=255)
           cACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACAt                                                                                     >  2:1966449/1‑100 (MQ=255)
           cACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACAt                                                                                     >  1:1727165/1‑100 (MQ=255)
            aCCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATg                                                                                    >  1:246121/1‑100 (MQ=255)
            aCCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATg                                                                                    <  1:1425271/100‑1 (MQ=255)
            aCCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATg                                                                                    >  1:1492930/1‑100 (MQ=255)
             cccGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATgg                                                                                   >  2:1088658/1‑100 (MQ=255)
                 gCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCtt                                                                               <  1:567102/100‑1 (MQ=255)
                   taTGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTg                                                                             <  1:175345/99‑1 (MQ=255)
                   aaTGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTg                                                                             <  2:178143/100‑1 (MQ=255)
                   aaTGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTg                                                                             <  2:350762/100‑1 (MQ=255)
                   aaTGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTg                                                                             <  2:467304/100‑1 (MQ=255)
                    aTGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTga                                                                            >  2:175342/1‑99 (MQ=255)
                     tGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGt                                                                            <  2:817324/99‑1 (MQ=255)
                     tGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTa                                                                           >  1:2304729/1‑100 (MQ=255)
                      gcgcgcgcTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTAc                                                                          >  1:279984/1‑100 (MQ=255)
                       cgcgcgcTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACt                                                                         >  2:1367312/1‑100 (MQ=255)
                           cgcTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGt                                                                     <  2:1346936/100‑1 (MQ=255)
                                 tcgtcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCcttc                                                               <  2:213792/100‑1 (MQ=255)
                                   gtcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAGCGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCcttctt                                                             <  1:2023252/100‑1 (MQ=255)
                                   gtcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCcttctt                                                             >  1:1039241/1‑100 (MQ=255)
                                   gtcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCcttctt                                                             <  2:335864/100‑1 (MQ=255)
                                    tcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCttcttc                                                            <  2:1358841/100‑1 (MQ=255)
                                    tcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCttcttc                                                            >  1:2738727/1‑100 (MQ=255)
                                    tcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCttcttc                                                            <  1:1219515/100‑1 (MQ=255)
                                        gggTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtg                                                        >  1:177677/1‑100 (MQ=255)
                                        gggTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtg                                                        >  2:1926174/1‑100 (MQ=255)
                                         ggTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtgg                                                       >  1:1995666/1‑100 (MQ=255)
                                         ggTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtgg                                                       >  1:1995657/1‑100 (MQ=255)
                                          gtaaGTCTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtgg                                                       <  1:403945/97‑1 (MQ=255)
                                          gTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTT‑ACCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtggt                                                      >  2:2045257/1‑100 (MQ=255)
                                          gTGAGTGGCAGCGTGTGATAACCTGCCGGCAGCTTCGACGGTAGATTGAGCGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtgg                                                       <  2:179492/99‑1 (MQ=255)
                                            gAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTca                                                    >  2:225337/1‑100 (MQ=255)
                                             aGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTcag                                                   <  1:1500855/100‑1 (MQ=255)
                                             aGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTcag                                                   >  1:1761990/1‑100 (MQ=255)
                                               tgtCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTcagca                                                 >  2:820085/1‑100 (MQ=255)
                                                  cAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCagcagc                                               >  2:173742/1‑100 (MQ=255)
                                                  cAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGcagcc                                              >  2:1775326/1‑100 (MQ=255)
                                                  cAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGcagcc                                              >  2:2534589/1‑100 (MQ=255)
                                                  cAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGcagcc                                              >  2:2534429/1‑100 (MQ=255)
                                                     cGTGTGATAACCTTCCGGCAGGTGCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCATTTGTACTGCGTTCCTTCTTCGGTGGTCAGcagccagc                                           <  2:980532/100‑1 (MQ=255)
                                                          gATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTatat                                      >  1:1496949/1‑100 (MQ=255)
                                                           aTAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGC‑TTGGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATAtt                                     <  1:1317507/100‑1 (MQ=255)
                                                            tAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCGCCGTTACATGGCCTTTGTAGTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTc                                    >  2:1510946/1‑100 (MQ=255)
                                                             aaCCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcg                                   >  1:1878543/1‑100 (MQ=255)
                                                             aaCCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcg                                   <  1:311022/100‑1 (MQ=255)
                                                             aaCCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcg                                   >  2:517433/1‑100 (MQ=255)
                                                              aCCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcgc                                  <  1:1367321/100‑1 (MQ=255)
                                                               ccTTCCGGCGGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGTTATATTcgcc                                 <  2:1155320/100‑1 (MQ=255)
                                                               ccTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcgcc                                 >  2:779531/1‑100 (MQ=255)
                                                                cTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcgccg                                <  2:1727151/100‑1 (MQ=255)
                                                                cTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcgccg                                >  1:298780/1‑100 (MQ=255)
                                                                 ttCCGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcgccgc                               <  1:2723933/100‑1 (MQ=255)
                                                                   ccGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTg                             >  1:2578576/1‑100 (MQ=255)
                                                                   ccGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTg                             >  1:2183570/1‑100 (MQ=255)
                                                                    cGGCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTGCTGCGTTCCTTCTTCGGTAGTCAGCAGCCAGCTATATTCGCCGCTGc                            <  1:946174/100‑1 (MQ=255)
                                                                     ggCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGcc                           >  2:1693804/1‑100 (MQ=255)
                                                                     ggCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGcc                           >  2:822290/1‑100 (MQ=255)
                                                                     ggCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTCCTTCGGTGGTCAGCAGCCCGCTATATTCGCCGCTGcc                           >  2:479972/1‑100 (MQ=255)
                                                                      gCAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGccc                          >  2:2270813/1‑100 (MQ=255)
                                                                       cAGCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCt                         <  2:2738708/100‑1 (MQ=255)
                                                                         gCTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTcc                       >  1:1186687/1‑100 (MQ=255)
                                                                          cgtcgtggGTAGATTGAACGCTTTG‑CCCCCGTTACGGGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGTTATATTcgcc                                 <  2:1155464/87‑1 (MQ=255)
                                                                          cTTCGTCGGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTcca                      <  1:825324/100‑1 (MQ=255)
                                                                            tcgtcgGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTccacc                    >  1:142138/1‑100 (MQ=255)
                                                                             cgtcgGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTccacca                   <  1:2300621/100‑1 (MQ=255)
                                                                                cgGTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTccaccacca                >  2:2359349/1‑100 (MQ=255)
                                                                                 ngTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCAt               >  2:941069/2‑100 (MQ=255)
                                                                                 ggTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCTGTGGTCAGCATCCAGCTATATTCGCCGCTGCCCTCCACCACCAt               >  2:2586565/1‑100 (MQ=255)
                                                                                 ggTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCAt               <  1:1772130/100‑1 (MQ=255)
                                                                                 ggTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCAt               >  1:500911/1‑100 (MQ=255)
                                                                                 ggTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCAt               >  2:940360/1‑100 (MQ=255)
                                                                                 ggTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCAt               >  2:654055/1‑100 (MQ=255)
                                                                                 ggTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGACAGCAGCCAGCTATATTCGCCGCTGCCCTccaccacca                >  2:2586180/1‑99 (MQ=255)
                                                                                  gTAGATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATc              >  2:215954/1‑100 (MQ=255)
                                                                                     gATTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGc           >  1:1813503/1‑100 (MQ=255)
                                                                                      aTTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCa          >  1:629559/1‑100 (MQ=255)
                                                                                      aTTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGACGCTGCCATCCACCACCATCGGCa          >  1:918104/1‑100 (MQ=255)
                                                                                      aTTGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTCCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCa          >  1:748552/1‑100 (MQ=255)
                                                                                       ttGAACGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCAt         >  2:1424257/1‑100 (MQ=255)
                                                                                          aaCGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCAtttt      >  2:2492188/1‑100 (MQ=255)
                                                                                          aaCGCTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCAtttt      >  1:1884858/1‑100 (MQ=255)
                                                                                              cTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTg  <  2:177674/100‑1 (MQ=255)
                                                                                              cTTTG‑CCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTg  <  1:2115589/100‑1 (MQ=255)
                                                                                                   |                                                                                               
GGGGCGACAATCACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTG  >  NC_000913/3549688‑3549882

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