Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,572,318 T→C 30.6% E395E (GAA→GAG)  glgB ← 1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,572,3180TC30.6% 216.5 / 136.2 163E395E (GAA→GAG) glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base T (53/60);  new base C (28/22);  total (81/82)
Fisher's exact test for biased strand distribution p-value = 3.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

Only 100 of 163 total aligned reads displayed.

CAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACG  >  NC_000913/3572221‑3572416
                                                                                                 |                                                                                                  
 naaTTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCa                                                                                                 >  2:1889287/2‑100 (MQ=255)
 aaaTTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCa                                                                                                 >  2:1731170/1‑100 (MQ=255)
 aaaTTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCa                                                                                                 <  1:625295/100‑1 (MQ=255)
  aaTTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCaa                                                                                                >  2:2579939/1‑100 (MQ=255)
   aTTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAAt                                                                                               <  1:1132179/100‑1 (MQ=255)
   aTTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAAt                                                                                               <  2:3459520/100‑1 (MQ=255)
       gttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAg                                                                                           >  1:135789/1‑100 (MQ=255)
       gttcgGGATCCACTCCCCCTCTTTACGGCTGGAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTCCAATCCAg                                                                                           <  1:3288022/100‑1 (MQ=255)
        ttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGt                                                                                          <  1:3344526/100‑1 (MQ=255)
           gggATCCACTCCCCCTCTTTACGGCTGTAGCCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTaa                                                                                        >  1:1378833/1‑99 (MQ=255)
             gATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAgc                                                                                     >  1:1514349/1‑100 (MQ=255)
              aTCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAgcg                                                                                    >  1:3414427/1‑100 (MQ=255)
              aTCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAgcg                                                                                    >  1:2947499/1‑100 (MQ=255)
              aTCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAACACCAAAACGTTCAATCCAGTAAAgcg                                                                                    >  1:2419845/1‑100 (MQ=255)
                ccACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGcgcg                                                                                  <  2:1389572/100‑1 (MQ=255)
                 cACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGt                                                                                 >  2:662903/1‑100 (MQ=255)
                 cACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGt                                                                                 >  2:742774/1‑100 (MQ=255)
                 cACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGt                                                                                 >  2:2313147/1‑100 (MQ=255)
                 cACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGt                                                                                 >  2:640524/1‑100 (MQ=255)
                  aCTCCCCCTCTTTACGGCTGTATTCGCGATATATCATGGACGCCACCTCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGtt                                                                                <  1:2331983/100‑1 (MQ=255)
                  aCTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGtt                                                                                <  1:432008/100‑1 (MQ=255)
                   cTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTa                                                                               >  2:2813037/1‑100 (MQ=255)
                    tCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTAc                                                                              <  2:3048697/100‑1 (MQ=255)
                     cccccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTAcc                                                                             >  1:279684/1‑100 (MQ=255)
                     cccccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTAcc                                                                             <  2:2406923/100‑1 (MQ=255)
                       cccTCTTTACGGCTGTAGTCGCGATAAAGCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACcg                                                                            >  1:2979459/1‑99 (MQ=255)
                        ccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACcgac                                                                          >  2:652213/1‑100 (MQ=255)
                         ctctTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACcgacg                                                                         >  1:461242/1‑100 (MQ=255)
                          tctTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACcgacga                                                                        <  2:766861/100‑1 (MQ=255)
                           ctTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAg                                                                       <  2:851528/100‑1 (MQ=255)
                           ctTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAg                                                                       <  2:3432328/100‑1 (MQ=255)
                           ctTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAg                                                                       <  2:1362066/100‑1 (MQ=255)
                           ctTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAg                                                                       <  1:3202717/100‑1 (MQ=255)
                            tttACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAgg                                                                      <  1:339518/100‑1 (MQ=255)
                               aCGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAg                                                                   <  2:3414453/100‑1 (MQ=255)
                               aCGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAg                                                                   <  1:21051/100‑1 (MQ=255)
                                cGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGt                                                                  <  2:2118671/100‑1 (MQ=255)
                                 ggCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGa                                                                     <  1:798431/96‑1 (MQ=255)
                                 ggCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGtt                                                                 <  2:1960794/100‑1 (MQ=255)
                                 ggCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGtt                                                                 <  2:33461/100‑1 (MQ=255)
                                  gCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTa                                                                >  1:3079639/1‑100 (MQ=255)
                                   cTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTAc                                                               <  1:1795335/100‑1 (MQ=255)
                                      tAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAATCGCTGTACCGACGAGGATGTTACtta                                                            >  2:508167/1‑98 (MQ=255)
                                       aGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCCGTAAAGCGCGTTACCGACGAGGAAGTTACTGAc                                                           >  2:2926156/1‑100 (MQ=255)
                                       aGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGa                                                            <  2:3167628/99‑1 (MQ=255)
                                       aGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGAc                                                           <  1:1700101/100‑1 (MQ=255)
                                       aGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGAc                                                           >  1:546842/1‑100 (MQ=255)
                                        gTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACt                                                          <  2:1646899/100‑1 (MQ=255)
                                        gTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACt                                                          >  1:2379027/1‑100 (MQ=255)
                                         tCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACtt                                                         >  1:1834024/1‑100 (MQ=255)
                                           gcgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCa                                                       >  2:232566/1‑100 (MQ=255)
                                           gcgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCAAGTAAAGCGAGTTACCGACGAGGAAGTTACTGACTTCa                                                       >  2:3050667/1‑100 (MQ=255)
                                           gcgATAAATCATTGACGCCACCACATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCa                                                       >  2:1463128/1‑100 (MQ=255)
                                            cgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCAc                                                      <  2:461246/100‑1 (MQ=255)
                                              aTAAATCATTGACGCCACCGCATCGACGAGCAGCGAATCAATACCAAAACGCTCAATCCAGTAAAGCGAGTTACCGACGAGGAAGTTACTGACTTCACgg                                                    >  1:3205460/1‑100 (MQ=255)
                                               tAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGc                                                   <  2:2232521/100‑1 (MQ=255)
                                               tAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGc                                                   >  1:2979530/1‑100 (MQ=255)
                                               tAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAGAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGc                                                   <  2:2486403/100‑1 (MQ=255)
                                                aaaTCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCg                                                  <  2:2602198/100‑1 (MQ=255)
                                                  aTCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCTCGGCGAc                                                <  2:1818191/100‑1 (MQ=255)
                                                  aTCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGAc                                                >  2:13388/1‑100 (MQ=255)
                                                  aTCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGAc                                                >  1:2230340/1‑100 (MQ=255)
                                                  aTCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGAc                                                <  2:3352668/100‑1 (MQ=255)
                                                      ttGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATa                                            >  2:1216346/1‑100 (MQ=255)
                                                      ttGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATa                                            <  1:2492236/100‑1 (MQ=255)
                                                         aCGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGtt                                         <  1:744196/100‑1 (MQ=255)
                                                          cGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTg                                        >  2:2456095/1‑100 (MQ=255)
                                                            ccACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTgaa                                      >  2:2068971/1‑98 (MQ=255)
                                                            ccACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTa                                      >  2:1447752/1‑100 (MQ=255)
                                                             cACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAg                                     >  2:3062009/1‑100 (MQ=255)
                                                               ccGCGTCGACGCGCAGTGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGAt                                   <  1:2025754/100‑1 (MQ=255)
                                                               ccGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGAt                                   <  1:640520/100‑1 (MQ=255)
                                                                cGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATc                                  >  2:3289126/1‑100 (MQ=255)
                                                                 gCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCa                                 >  1:1943323/1‑100 (MQ=255)
                                                                  cATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAg                                >  2:2851846/1‑100 (MQ=255)
                                                                   aTCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGc                               >  1:2595875/1‑100 (MQ=255)
                                                                    tCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGcn                              <  1:1889271/100‑2 (MQ=255)
                                                                     cGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCgt                             >  1:1014179/1‑100 (MQ=255)
                                                                       aCGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCgtgt                           <  1:1199507/100‑1 (MQ=255)
                                                                        cgcgCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGtt                          <  2:102171/100‑1 (MQ=255)
                                                                        cgcgCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGtt                          >  2:1231483/1‑100 (MQ=255)
                                                                             agcgcaTCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAg                      <  1:2926703/99‑1 (MQ=255)
                                                                               cgcaTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTcc                   <  1:184229/100‑1 (MQ=255)
                                                                               cacaTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTcc                   <  1:1632411/98‑1 (MQ=255)
                                                                                gcaTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCt                  >  1:3212597/1‑100 (MQ=255)
                                                                                   tCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgat               >  2:759874/1‑100 (MQ=255)
                                                                                   tCAATACCAAAACGCTCAATCTAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgat               >  1:1087642/1‑100 (MQ=255)
                                                                                   tCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgat               >  1:531065/1‑100 (MQ=255)
                                                                                    cAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgatg              <  1:456179/100‑1 (MQ=255)
                                                                                    cAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgatg              >  1:2622642/1‑100 (MQ=255)
                                                                                      aTACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTgatgat            <  1:819735/100‑1 (MQ=255)
                                                                                       tACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATa           <  1:2159897/100‑1 (MQ=255)
                                                                                        aCCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAg          >  1:2653209/1‑100 (MQ=255)
                                                                                         ccAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGc         >  1:1305838/1‑100 (MQ=255)
                                                                                          cAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGcc        <  2:2155842/100‑1 (MQ=255)
                                                                                           aaaaCGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCt       <  1:2313129/100‑1 (MQ=255)
                                                                                           aaaaCGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCt       >  1:1302989/1‑100 (MQ=255)
                                                                                            aaaCGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTATAGATCAGCGTGTTCCAGTCCTGATGATAAcc        >  1:1034313/1‑98 (MQ=255)
                                                                                              aCGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCa    <  1:652209/100‑1 (MQ=255)
                                                                                                gTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCAcg  >  2:2633575/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
CAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACG  >  NC_000913/3572221‑3572416

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