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)
           cACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACAt                                                                                     >  2:1966449/1‑100 (MQ=255)
           cACCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACAt                                                                                     >  1:1727165/1‑100 (MQ=255)
            aCCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATg                                                                                    >  1:246121/1‑100 (MQ=255)
            aCCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATg                                                                                    <  1:1425271/100‑1 (MQ=255)
            aCCCGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATg                                                                                    >  1:1492930/1‑100 (MQ=255)
             cccGGCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATgg                                                                                   >  2:1088658/1‑100 (MQ=255)
                 gCAATGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCtt                                                                               <  1:567102/100‑1 (MQ=255)
                   taTGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTg                                                                             <  1:175345/99‑1 (MQ=255)
                   aaTGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTg                                                                             <  2:178143/100‑1 (MQ=255)
                   aaTGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTg                                                                             <  2:350762/100‑1 (MQ=255)
                   aaTGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTg                                                                             <  2:467304/100‑1 (MQ=255)
                    aTGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTga                                                                            >  2:175342/1‑99 (MQ=255)
                     tGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGt                                                                            <  2:817324/99‑1 (MQ=255)
                     tGCGCGCGCTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTa                                                                           >  1:2304729/1‑100 (MQ=255)
                      gcgcgcgcTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTAc                                                                          >  1:279984/1‑100 (MQ=255)
                       cgcgcgcTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACt                                                                         >  2:1367312/1‑100 (MQ=255)
                           cgcTGGTCGTCCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGt                                                                     <  2:1346936/100‑1 (MQ=255)
                                 tcgtcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCcttc                                                               <  2:213792/100‑1 (MQ=255)
                                   gtcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAGCGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCcttctt                                                             <  1:2023252/100‑1 (MQ=255)
                                   gtcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCcttctt                                                             >  1:1039241/1‑100 (MQ=255)
                                   gtcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCcttctt                                                             <  2:335864/100‑1 (MQ=255)
                                    tcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCttcttc                                                            <  2:1358841/100‑1 (MQ=255)
                                    tcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCttcttc                                                            >  1:2738727/1‑100 (MQ=255)
                                    tcCTGGGTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCttcttc                                                            <  1:1219515/100‑1 (MQ=255)
                                        gggTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtg                                                        >  1:177677/1‑100 (MQ=255)
                                        gggTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtg                                                        >  2:1926174/1‑100 (MQ=255)
                                         ggTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtgg                                                       >  1:1995666/1‑100 (MQ=255)
                                         ggTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtgg                                                       >  1:1995657/1‑100 (MQ=255)
                                          gtaaGTCTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtgg                                                       <  1:403945/97‑1 (MQ=255)
                                          gTGAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTACCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtggt                                                      >  2:2045257/1‑100 (MQ=255)
                                          gTGAGTGGCAGCGTGTGATAACCTGCCGGCAGCTTCGACGGTAGATTGAGCGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCggtgg                                                       <  2:179492/99‑1 (MQ=255)
                                            gAGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTca                                                    >  2:225337/1‑100 (MQ=255)
                                             aGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTcag                                                   <  1:1500855/100‑1 (MQ=255)
                                             aGTGTCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTcag                                                   >  1:1761990/1‑100 (MQ=255)
                                               tgtCAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTcagca                                                 >  2:820085/1‑100 (MQ=255)
                                                  cAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCagcagc                                               >  2:173742/1‑100 (MQ=255)
                                                  cAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGcagcc                                              >  2:1775326/1‑100 (MQ=255)
                                                  cAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGcagcc                                              >  2:2534589/1‑100 (MQ=255)
                                                  cAGCGTGTGATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGcagcc                                              >  2:2534429/1‑100 (MQ=255)
                                                     cGTGTGATAACCTTCCGGCAGGTGCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCATTTGTACTGCGTTCCTTCTTCGGTGGTCAGcagccagc                                           <  2:980532/100‑1 (MQ=255)
                                                          gATAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTatat                                      >  1:1496949/1‑100 (MQ=255)
                                                           aTAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTGGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATAtt                                     <  1:1317507/100‑1 (MQ=255)
                                                            tAACCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCGCCGTTACATGGCCTTTGTAGTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTc                                    >  2:1510946/1‑100 (MQ=255)
                                                             aaCCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcg                                   >  1:1878543/1‑100 (MQ=255)
                                                             aaCCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcg                                   <  1:311022/100‑1 (MQ=255)
                                                             aaCCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcg                                   >  2:517433/1‑100 (MQ=255)
                                                              aCCTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcgc                                  <  1:1367321/100‑1 (MQ=255)
                                                               ccTTCCGGCGGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGTTATATTcgcc                                 <  2:1155320/100‑1 (MQ=255)
                                                               ccTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcgcc                                 >  2:779531/1‑100 (MQ=255)
                                                                cTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcgccg                                <  2:1727151/100‑1 (MQ=255)
                                                                cTTCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcgccg                                >  1:298780/1‑100 (MQ=255)
                                                                 ttCCGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTcgccgc                               <  1:2723933/100‑1 (MQ=255)
                                                                   ccGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTg                             >  1:2578576/1‑100 (MQ=255)
                                                                   ccGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTg                             >  1:2183570/1‑100 (MQ=255)
                                                                    cGGCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTGCTGCGTTCCTTCTTCGGTAGTCAGCAGCCAGCTATATTCGCCGCTGc                            <  1:946174/100‑1 (MQ=255)
                                                                     ggCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGcc                           >  2:1693804/1‑100 (MQ=255)
                                                                     ggCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGcc                           >  2:822290/1‑100 (MQ=255)
                                                                     ggCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTCCTTCGGTGGTCAGCAGCCCGCTATATTCGCCGCTGcc                           >  2:479972/1‑100 (MQ=255)
                                                                      gCAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGccc                          >  2:2270813/1‑100 (MQ=255)
                                                                       cAGCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCt                         <  2:2738708/100‑1 (MQ=255)
                                                                         gCTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTcc                       >  1:1186687/1‑100 (MQ=255)
                                                                          cgtcgtggGTAGATTGAACGCTTTGCCCCCGTTACGGGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGTTATATTcgcc                                 <  2:1155464/87‑1 (MQ=255)
                                                                          cTTCGTCGGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTcca                      <  1:825324/100‑1 (MQ=255)
                                                                            tcgtcgGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTccacc                    >  1:142138/1‑100 (MQ=255)
                                                                             cgtcgGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTccacca                   <  1:2300621/100‑1 (MQ=255)
                                                                                cgGTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTccaccacca                >  2:2359349/1‑100 (MQ=255)
                                                                                 ngTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCAt               >  2:941069/2‑100 (MQ=255)
                                                                                 ggTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCTGTGGTCAGCATCCAGCTATATTCGCCGCTGCCCTCCACCACCAt               >  2:2586565/1‑100 (MQ=255)
                                                                                 ggTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCAt               <  1:1772130/100‑1 (MQ=255)
                                                                                 ggTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCAt               >  1:500911/1‑100 (MQ=255)
                                                                                 ggTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCAt               >  2:940360/1‑100 (MQ=255)
                                                                                 ggTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCAt               >  2:654055/1‑100 (MQ=255)
                                                                                 ggTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGACAGCAGCCAGCTATATTCGCCGCTGCCCTccaccacca                >  2:2586180/1‑99 (MQ=255)
                                                                                  gTAGATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATc              >  2:215954/1‑100 (MQ=255)
                                                                                     gATTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGc           >  1:1813503/1‑100 (MQ=255)
                                                                                      aTTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCa          >  1:629559/1‑100 (MQ=255)
                                                                                      aTTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGACGCTGCCATCCACCACCATCGGCa          >  1:918104/1‑100 (MQ=255)
                                                                                      aTTGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTCCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCa          >  1:748552/1‑100 (MQ=255)
                                                                                       ttGAACGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCAt         >  2:1424257/1‑100 (MQ=255)
                                                                                          aaCGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCAtttt      >  2:2492188/1‑100 (MQ=255)
                                                                                          aaCGCTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCAtttt      >  1:1884858/1‑100 (MQ=255)
                                                                                              cTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTg  <  2:177674/100‑1 (MQ=255)
                                                                                              cTTTGCCCCCGTTACATGGCCTTTGTACTGCGTTCCTTCTTCGGTGGTCAGCAGCCAGCTATATTCGCCGCTGCCCTCCACCACCATCGGCATTTTTTTg  <  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: