Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,572,318 T→C 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,3180TC98.0% 346.2 / ‑3.3 102E395E (GAA→GAG) glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (43/57);  total (44/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 102 total aligned reads displayed.

AATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTT  >  NC_000913/3572223‑3572412
                                                                                               |                                                                                              
aaTTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCaa                                                                                            >  2:2103227/1‑100 (MQ=255)
aaTTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCaa                                                                                            <  2:1628116/100‑1 (MQ=255)
   tcgttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATcc                                                                                         >  2:1390101/1‑100 (MQ=255)
    cgttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATcc                                                                                         <  2:1377815/99‑1 (MQ=255)
     gttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAg                                                                                       <  1:2868877/100‑1 (MQ=255)
     gttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAg                                                                                       <  2:494884/100‑1 (MQ=255)
       tcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTa                                                                                     >  1:2650040/1‑100 (MQ=255)
       tcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTa                                                                                     <  1:985240/100‑1 (MQ=255)
        cgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGt                                                                                      <  2:836174/98‑1 (MQ=255)
         gggATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGCCGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTaa                                                                                    >  1:887036/1‑99 (MQ=255)
         gggATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTaaa                                                                                   >  2:1124378/1‑100 (MQ=255)
         gggATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTaaa                                                                                   >  1:2668920/1‑100 (MQ=255)
          ggATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAg                                                                                  >  1:233984/1‑100 (MQ=255)
           gATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAgc                                                                                 <  2:2842641/100‑1 (MQ=255)
              ccACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGcgcg                                                                              >  2:2269658/1‑100 (MQ=255)
               cACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCGGTAAAGCGCGt                                                                             >  1:259607/1‑100 (MQ=255)
               cACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGt                                                                             <  1:1209703/100‑1 (MQ=255)
               cACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGt                                                                             <  1:420724/100‑1 (MQ=255)
                  tCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCACTAAAGCGCGTTAc                                                                          <  1:1720968/100‑1 (MQ=255)
                   cccccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAATGCGCGTTAcc                                                                         >  1:441145/1‑100 (MQ=255)
                   cccccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTAcc                                                                         >  2:2775483/1‑100 (MQ=255)
                   cccccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTAcc                                                                         <  1:505328/100‑1 (MQ=255)
                     cccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACcga                                                                       >  1:2418862/1‑100 (MQ=255)
                     cccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACcga                                                                       >  1:1457659/1‑100 (MQ=255)
                        tctTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACcgacga                                                                    >  1:1624790/1‑100 (MQ=255)
                          tttACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAgg                                                                  <  2:2590224/100‑1 (MQ=255)
                          tttACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAgg                                                                  <  2:2565727/100‑1 (MQ=255)
                              cGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGt                                                              <  2:2548314/100‑1 (MQ=255)
                                 cTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTAc                                                           >  1:2548056/1‑100 (MQ=255)
                                 cTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTAc                                                           >  2:811343/1‑100 (MQ=255)
                                  tGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACt                                                          <  2:2068660/100‑1 (MQ=255)
                                   gTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTg                                                         <  2:2987259/100‑1 (MQ=255)
                                   gTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTg                                                         <  1:2503411/100‑1 (MQ=255)
                                   gTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTg                                                         <  1:1056010/100‑1 (MQ=255)
                                     aGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGAc                                                       <  1:1493793/100‑1 (MQ=255)
                                      gTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACt                                                      >  2:2656403/1‑100 (MQ=255)
                                        cgcgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTc                                                    <  1:2486711/100‑1 (MQ=255)
                                         gcgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCa                                                   >  2:15688/1‑100 (MQ=255)
                                          cgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCAc                                                  >  1:1260773/1‑100 (MQ=255)
                                          cgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCAc                                                  <  1:1345880/100‑1 (MQ=255)
                                             tAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGc                                               <  1:2138514/100‑1 (MQ=255)
                                             tAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGc                                               <  1:681371/100‑1 (MQ=255)
                                               aaTCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGa                                             >  1:2554723/1‑100 (MQ=255)
                                                aTCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGAc                                            <  2:2668896/100‑1 (MQ=255)
                                                aTCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGAc                                            <  2:2994355/100‑1 (MQ=255)
                                                   aTTGACGCCACCGCATCGACGCGCAGCGCATCAATGCCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCa                                          <  2:428814/99‑1 (MQ=255)
                                                   aTTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGATTTCACGGCGACCAt                                         <  2:1610795/100‑1 (MQ=255)
                                                   aTTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCAt                                         >  2:1548637/1‑100 (MQ=255)
                                                    ttGACTCCACCGCATCGACGCGCAGCGCATCAAGACCAAAACGTTCAATTCAGTAATGTGGGTGACCGACGAGGAAGTTACTGACTTCACGGCGACCATa                                        >  2:1811246/1‑100 (MQ=255)
                                                    ttGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATa                                        >  1:2313009/1‑100 (MQ=255)
                                                    ttGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATa                                        >  2:1386543/1‑100 (MQ=255)
                                                      gACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGCCGAGGAAGTCACTGACTTCACGGCCAc                                            >  2:239035/1‑94 (MQ=255)
                                                      gACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGt                                      <  2:1700211/100‑1 (MQ=255)
                                                      gACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGt                                      <  2:1967694/100‑1 (MQ=255)
                                                      gACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGt                                      <  2:1968354/100‑1 (MQ=255)
                                                      gACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGt                                      <  2:400172/100‑1 (MQ=255)
                                                        cGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTg                                    <  1:217782/100‑1 (MQ=255)
                                                        cGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACCCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTg                                    >  2:1870670/1‑100 (MQ=255)
                                                          ccACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTa                                  <  2:2554700/100‑1 (MQ=255)
                                                          ccACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTa                                  <  1:736362/100‑1 (MQ=255)
                                                          ccACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTa                                  <  1:2544523/100‑1 (MQ=255)
                                                            aCCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGa                                <  2:2418840/100‑1 (MQ=255)
                                                              cGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATc                              <  1:1471070/100‑1 (MQ=255)
                                                               gCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCa                             <  1:1178227/100‑1 (MQ=255)
                                                                cATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAg                            >  1:653862/1‑100 (MQ=255)
                                                                cATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAg                            >  1:1240414/1‑100 (MQ=255)
                                                                 atctgcacgcagcgcaTCTATACCAAATCGCTCAATCCAGTAAAGTGAGTTATTGACGTGGAAGTTTTTGTTTTAACGGCGACCATAGTTGTAGATc                              <  2:2940096/90‑1 (MQ=255)
                                                                 aTCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGc                           <  2:638438/100‑1 (MQ=255)
                                                                  tCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCg                          >  2:911191/1‑100 (MQ=255)
                                                                  tCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCg                          <  1:2678452/100‑1 (MQ=255)
                                                                   cGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCgt                         >  1:1271093/1‑100 (MQ=255)
                                                                    gACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCgtg                        <  1:2218571/100‑1 (MQ=255)
                                                                       gcgcagcgcaTCAATACCAAAACGCTCAATCCAGTAACGCGCGTTACCGACGAGGCAGTTACTCACTTCCCGGCGACCATAGTTGTAGATCAGCGTGtt                      >  2:1330860/1‑99 (MQ=255)
                                                                        cgcagcgcaTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTcc                    >  2:1553837/1‑100 (MQ=255)
                                                                        cgcagcgcaTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTcc                    <  1:57502/100‑1 (MQ=255)
                                                                         gcagcgcaTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCa                   >  1:2620736/1‑100 (MQ=255)
                                                                          cagcgcaTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAg                  <  2:233981/100‑1 (MQ=255)
                                                                          cagcgcaTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAg                  <  2:1399717/100‑1 (MQ=255)
                                                                          cagcgcaTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAg                  <  2:1399671/100‑1 (MQ=255)
                                                                            gcgcaTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTc                <  2:2848227/100‑1 (MQ=255)
                                                                             cgcaTCAATACCAAAACGCTCAAACCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTcc               <  2:1055414/100‑1 (MQ=255)
                                                                              gcaTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCt              >  2:487892/1‑100 (MQ=255)
                                                                               caTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGGTCAGCGTGTTCCAGTCCtg             <  2:1563791/100‑1 (MQ=255)
                                                                               caTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtg             <  1:3027025/100‑1 (MQ=255)
                                                                               caTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtg             >  2:603927/1‑100 (MQ=255)
                                                                                aTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtga            <  2:1457647/100‑1 (MQ=255)
                                                                                aTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtga            <  2:441141/100‑1 (MQ=255)
                                                                                aTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtg             <  1:1390113/99‑1 (MQ=255)
                                                                                 tCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgat           >  1:2927912/1‑100 (MQ=255)
                                                                                  cAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgatg          >  1:469023/1‑100 (MQ=255)
                                                                                  cAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgatg          <  1:20728/100‑1 (MQ=255)
                                                                                    aTACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTgatgat        <  1:1435435/100‑1 (MQ=255)
                                                                                     tACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATa       >  1:2412693/1‑100 (MQ=255)
                                                                                     tACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATa       <  1:2656426/100‑1 (MQ=255)
                                                                                       ccAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGc     >  1:2288034/1‑100 (MQ=255)
                                                                                       ccAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGc     >  1:1577391/1‑100 (MQ=255)
                                                                                        cAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGcc    >  2:1913454/1‑100 (MQ=255)
                                                                                         aaaaCGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATGGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCt   >  2:2038788/1‑100 (MQ=255)
                                                                                         aaaaCGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTCGATCAGCGTGTTCCAGTCCTGATGATAGCCt   >  2:2038810/1‑100 (MQ=255)
                                                                                          aaaCGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCtt  <  2:1577378/100‑1 (MQ=255)
                                                                                               |                                                                                              
AATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTT  >  NC_000913/3572223‑3572412

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: ‑