Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,548,507 A→G 6.2% F522L (TTC→CTC)  malQ ← 4‑alpha‑glucanotransferase (amylomaltase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,548,5070AG 6.2% 182.0 / 2.5 64F522L (TTC→CTC) malQ4‑alpha‑glucanotransferase (amylomaltase)
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base G (2/2);  total (31/33)
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

CAGCCCCAGGGTTTTGCCCAGCGTTAGATCCCCGCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGAT  >  NC_000913/3548364‑3548653
                                                                                                                                               |                                                                                                                                                  
cAGCCCCAGGGTTTTGCCCAGCGTTAGATCCCCGCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATaga                                                                                                                                               >  2:1280472/1‑149 (MQ=255)
cAGCCCCAGGGTTTTGCCCAGCGTTAGATCCCCGCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGGCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGGCTTCTCGTGGGCGTTTTCGAAATaga                                                                                                                                               <  2:196861/149‑1 (MQ=255)
   ccccAGGGTTTTGCCCAGCGTTAGATCCCCGCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCTTCGCCTGCTGCGGATACGCTGTCGGTGCCCGGACCGTCTTCTCGTGGGCGTTTTCGAAATAGAGcc                                                                                                                                            >  1:136054/1‑148 (MQ=255)
   ccccAGGGTTTTGCCCAGCGTTAGATCCCCGCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAGATAGAGCa                                                                                                                                            >  2:1040720/1‑149 (MQ=255)
   ccccAGGGTTTTGCCCAGCGTTAGATCCCCGCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCACGTGGTCGTTTTCGAGATAGAGCa                                                                                                                                            >  2:1036606/1‑149 (MQ=255)
         ggTTTTGCCCAGCGTTAGATCCCCGCACTCCCAGTACCCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGt                                                                                                                                      <  1:1524614/149‑1 (MQ=255)
              tGCCCAGCGTTAGATCCCCGCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAagag                                                                                                                                 <  1:1208531/149‑1 (MQ=255)
                cccAGCGTTAGATCCCCGCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTa                                                                                                                               <  1:659942/149‑1 (MQ=255)
                  cAGCGTTAGATCCCCGCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCTTGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCGTTCGGTGCCCGGACCGTCTTCTCGAGGTCGTTGTCGAAATAGAGCACTTTGTAAGAGTaca                                                                                                                             >  1:476641/1‑149 (MQ=255)
                       ttAGATCCCCGCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCCCGGAACGTCTTCTCTTGGTCGTTTTCGAAATAGACCACTTTGTAAGAGTACACACCg                                                                                                                        >  2:1375664/1‑149 (MQ=255)
                             ccccGCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGctgctg                                                                                                                  >  2:564015/1‑149 (MQ=255)
                              cccGCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTTTAAGAGTACACACCGCttctgc                                                                                                                 >  1:1086785/1‑149 (MQ=255)
                               cagCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCTGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTgcg                                                                                                                <  1:682533/147‑1 (MQ=255)
                                 gCACTCCCAGTAACCGCGCAGCGTTGGCAGTTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCa                                                                                                              <  1:574696/149‑1 (MQ=255)
                                 gCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGCCTGCTCCGGATACGCTTTCGGTGCACGGAACGTCGTCTCGTGGTCGTTTTCGAAATAGCGCACTTTGTAAGAGTACACACCGCTGCTGCGCa                                                                                                              >  1:342726/1‑149 (MQ=255)
                                  cACTCCCAGTATCCGCGCATCGGTGGCAGGGCTTGTGTCGCCGCAACCGCCATCGAGTGCTGCGTCTACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAATTAGAGCTCTTTGTAAGTGTACACACCGCTGCTGCGCAg                                                                                                             <  2:1321982/149‑1 (MQ=255)
                                  cACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCATCCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAg                                                                                                             <  2:175290/149‑1 (MQ=255)
                                   aCTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGGGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGc                                                                                                            <  2:735461/149‑1 (MQ=255)
                                     tCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCtt                                                                                                          <  2:407474/149‑1 (MQ=255)
                                     tCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCtt                                                                                                          >  2:403297/1‑149 (MQ=255)
                                       ccAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTAc                                                                                                        <  1:822368/149‑1 (MQ=255)
                                        cAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTAcc                                                                                                       <  2:1567194/149‑1 (MQ=255)
                                           tAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCgcc                                                                                                    >  2:1480075/1‑147 (MQ=255)
                                              ccGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAAt                                                                                                 <  2:102127/149‑1 (MQ=255)
                                              ccGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAAt                                                                                                 <  2:103579/149‑1 (MQ=255)
                                                   cAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTa                                                                                            <  2:1524384/149‑1 (MQ=255)
                                                    aGCGTTGGCAGTTCATGTGTCGCCGCAACAGCCATCGACTGCTCCGGATACGCTTACGGTTCACGGAACGTCTTCTCGTCGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTAc                                                                                           <  1:1347036/149‑1 (MQ=255)
                                                    aGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTACGGATACGCTTTCGGTGCACGGACCGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTCCACACCGCTGCTGCGCAGCTTACCGACAATCTCTAc                                                                                           >  2:702161/1‑149 (MQ=255)
                                                    aGCGTTGGCAGGGCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTAc                                                                                           <  2:1246932/149‑1 (MQ=255)
                                                        ttGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACcggt                                                                                       <  2:1466653/149‑1 (MQ=255)
                                                           gCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTGTCGGTGCACGGAACGTCTTCTGGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACcggtacg                                                                                    <  2:1331445/149‑1 (MQ=255)
                                                                  aTGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACcg                                                                                         <  1:7425/137‑1 (MQ=255)
                                                                  aTGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACcg                                                                                         >  2:7425/1‑137 (MQ=255)
                                                                   tgtgTCGCCTCAACCGCCATCGACGGCTCCGGATACGCGTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTTCGCAGCTTACCGACAATCTCTACCGGTACGGTACCgag                                                                            <  2:725061/149‑1 (MQ=255)
                                                                      gtCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATc                                                                         <  1:967552/149‑1 (MQ=255)
                                                                         gccgcAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCgagctgtta                                                                      >  2:1525141/1‑143 (MQ=255)
                                                                         gccgcAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACAGACAATCTCTACCGGTACGGTACCGAGATCTTc                                                                      >  2:350166/1‑149 (MQ=255)
                                                                           cgcAACCGCCATCGACTGCTCCGTATACGCTTTCGGTTCACGGAACGTCTTCTCGTGGTCGTTTTCGAGATACAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCAc                                                                    <  1:1036606/149‑1 (MQ=255)
                                                                           cgcAACCGCAATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAGATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCAc                                                                    <  1:1040720/149‑1 (MQ=255)
                                                                                 cGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCa                                                              >  1:164624/1‑149 (MQ=255)
                                                                                         aCTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAg                                                      <  1:1560069/149‑1 (MQ=255)
                                                                                              tCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATg                                                 >  2:187791/1‑149 (MQ=255)
                                                                                               ccGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGa                                                >  1:1554086/1‑149 (MQ=255)
                                                                                                   aTACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGtt                                            <  1:986773/149‑1 (MQ=255)
                                                                                                   aTACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGtt                                            <  1:187791/149‑1 (MQ=255)
                                                                                                    tACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTAGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGttt                                           <  1:1280472/149‑1 (MQ=255)
                                                                                                     aCGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTa                                          >  1:1241717/1‑149 (MQ=255)
                                                                                                       gCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGATTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCCCCAAACACCATACAGAGATGCCGTTTACt                                        >  2:1164322/1‑149 (MQ=255)
                                                                                                         gttCGGTGCACGGAACGTCTTCTCGTGGTCGTATTCGAAATAGCGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACttt                                      <  1:589562/148‑1 (MQ=255)
                                                                                                              gTGCAGGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTGGTAAGAGTACACACCGCTGCTGCGCAGCTTACCTACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGt                                 <  1:1489448/149‑1 (MQ=255)
                                                                                                              gTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGt                                 >  2:1147818/1‑149 (MQ=255)
                                                                                                                    ggAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCCAGTGCCAGa                           >  1:656451/1‑149 (MQ=255)
                                                                                                                      aaCGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAAt                         >  1:1535169/1‑149 (MQ=255)
                                                                                                                       aCGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCAACCGGTACGGTACCGAGATATTCACCAAAAACCATACACCGATGACCTTTACTTACCAGTGACAGAATc                        >  2:280254/1‑149 (MQ=255)
                                                                                                                       aCGTCTTCTCGTGGTCGTGTTAGAACTCGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTTCAGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATc                        <  1:1373396/149‑1 (MQ=255)
                                                                                                                          tcttctCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGCTGACGTTTACTTTCGAGTGCCAGAATCgag                     >  2:1244440/1‑149 (MQ=255)
                                                                                                                               tcGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACCGCGATGACGTGTACTTTCGAGTGCCAGAATCGAGAGCAg                >  2:217817/1‑149 (MQ=255)
                                                                                                                                 gTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGat              >  2:64082/1‑149 (MQ=255)
                                                                                                                                   ggTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGatca            <  1:675522/149‑1 (MQ=255)
                                                                                                                                    gTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAgatcagatcaa           >  1:374549/1‑148 (MQ=255)
                                                                                                                                    gTCGTTTTAGAAATAGAGCACTTTCTAAGAGTACACACCGCTGCTGCGCAGCTTACAGACAAACTCTACCGGTAAGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAAACGAGAGCAGatcat           >  2:815263/1‑149 (MQ=255)
                                                                                                                                       gTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTCCCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCa        <  1:1270173/149‑1 (MQ=255)
                                                                                                                                        ttttCGAAATAGAGCACTTTGTTAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCAc       >  1:539911/1‑149 (MQ=255)
                                                                                                                                             gAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCCCCAATCACCATACAGCGATGCCGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGAt  >  1:384026/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CAGCCCCAGGGTTTTGCCCAGCGTTAGATCCCCGCACTCCCAGTAACCGCGCAGCGTTGGCAGGTCATGTGTCGCCGCAACCGCCATCGACTGCTCCGGATACGCTTTCGGTGCACGGAACGTCTTCTCGTGGTCGTTTTCGAAATAGAGCACTTTGTAAGAGTACACACCGCTGCTGCGCAGCTTACCGACAATCTCTACCGGTACGGTACCGAGATCTTCACCAATCACCATACAGCGATGACGTTTACTTTCGAGTGCCAGAATCGAGAGCAGATCATCCACCGGAT  >  NC_000913/3548364‑3548653

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: