Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,572,318 T→C 20.1% 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,3180TC20.1% 212.7 / 47.8 104E395E (GAA→GAGglgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base T (37/46);  new base C (10/11);  total (47/57)
Fisher's exact test for biased strand distribution p-value = 8.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 108 total aligned reads displayed.

GCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACG  >  NC_000913/3572219‑3572416
                                                                                                   |                                                                                                  
gccAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGt                                                                                                    >  1:968961/1‑100 (MQ=255)
gccAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGt                                                                                                    >  2:1912728/1‑100 (MQ=255)
 ccAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGtt                                                                                                   <  2:797718/100‑1 (MQ=255)
  cAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGtt                                                                                                   >  2:1711309/1‑99 (MQ=255)
  cAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTc                                                                                                  <  1:2083786/100‑1 (MQ=255)
  cAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTc                                                                                                  <  2:1914134/100‑1 (MQ=255)
  cAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTc                                                                                                  <  1:85408/100‑1 (MQ=255)
  cAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTc                                                                                                  <  2:1050942/100‑1 (MQ=255)
  cAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTc                                                                                                  <  2:728490/100‑1 (MQ=255)
      ntcgttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATc                                                                                              >  2:2585323/2‑100 (MQ=255)
       tcgttcgGGATCCTCTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATcc                                                                                             <  1:1910072/100‑1 (MQ=255)
       tcgttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATcc                                                                                             <  2:1182162/100‑1 (MQ=255)
         gttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAg                                                                                           >  1:1557253/1‑100 (MQ=255)
         gttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAg                                                                                           >  1:64725/1‑100 (MQ=255)
         gttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAg                                                                                           <  2:2005107/100‑1 (MQ=255)
         gttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAg                                                                                           >  2:2010924/1‑100 (MQ=255)
          ttcgGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGAAGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGt                                                                                          >  2:2055848/1‑100 (MQ=255)
            cgGGATCCACTCCCCCTCTTTACGGCTGTATTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCACATCAATACCAAAACGTTCAATCCAGTaa                                                                                        >  1:1559371/1‑100 (MQ=255)
              ggATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAg                                                                                      <  2:634988/100‑1 (MQ=255)
              ggATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAg                                                                                      <  2:191904/100‑1 (MQ=255)
              ggATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAg                                                                                      >  2:148126/1‑100 (MQ=255)
              ggATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAg                                                                                      >  1:179747/1‑100 (MQ=255)
                  ccACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGcgcg                                                                                  <  2:171555/100‑1 (MQ=255)
                   cACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGt                                                                                 >  2:791013/1‑100 (MQ=255)
                   cACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGt                                                                                 >  1:1179034/1‑100 (MQ=255)
                   cACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGt                                                                                 >  1:759985/1‑100 (MQ=255)
                     cTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTa                                                                               <  1:1770257/100‑1 (MQ=255)
                     cTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCACAGAACATAAATACCAAAACGTTCAATCCAGTAAAGCGCGTTa                                                                               >  1:1042111/1‑100 (MQ=255)
                      tCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTAc                                                                              <  2:936654/100‑1 (MQ=255)
                        ccccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACcg                                                                            >  2:551387/1‑100 (MQ=255)
                        ccccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACcg                                                                            >  1:1866725/1‑100 (MQ=255)
                        ccccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACcg                                                                            >  1:1894385/1‑100 (MQ=255)
                         cccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACcga                                                                           >  1:496684/1‑100 (MQ=255)
                          ccTCTTTACGGCTGTAGTCGCGATATATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACcgac                                                                          <  1:419208/100‑1 (MQ=255)
                          ccTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACcgac                                                                          <  2:850006/100‑1 (MQ=255)
                           ctctTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACcgacg                                                                         <  1:1965400/100‑1 (MQ=255)
                              tttACGGCTGTAGTCGGTTTGTAGTGTTGACGCCAACTCTTTGATGCGTATCGCTTCATTACCATTTCGTTCAATCCAGTAAAg                                                                                      <  2:1669010/84‑1 (MQ=255)
                                 aCGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAg                                                                   <  1:932225/100‑1 (MQ=255)
                                  cGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGt                                                                  <  2:1204615/100‑1 (MQ=255)
                                    gCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGATACCGACGAGGAAGTTa                                                                >  2:1732719/1‑100 (MQ=255)
                                     cTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGttgc                                                               <  1:1801637/100‑3 (MQ=255)
                                     cTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTAc                                                               <  1:1804061/100‑1 (MQ=255)
                                      tGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACt                                                              <  1:1895470/100‑1 (MQ=255)
                                         aGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGAc                                                           <  1:184024/100‑1 (MQ=255)
                                          gTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACt                                                          >  2:727211/1‑100 (MQ=255)
                                            cgcgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTc                                                        >  1:601229/1‑100 (MQ=255)
                                            cgcgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTc                                                        <  2:1248749/100‑1 (MQ=255)
                                            cgcgATAAATCATTGACGCCACCCCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAATAAAAAGCGTTACCGACGAGGAAGTTACTGACTTc                                                        >  1:1417320/1‑100 (MQ=255)
                                             gcgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCa                                                       <  2:2187033/100‑1 (MQ=255)
                                              cgATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCAc                                                      <  2:31072/100‑1 (MQ=255)
                                               gATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACg                                                     >  1:2561803/1‑100 (MQ=255)
                                                 tAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGc                                                   <  1:1284173/100‑1 (MQ=255)
                                                    aTCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGAc                                                <  1:1083082/100‑1 (MQ=255)
                                                    aTCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGAc                                                <  1:2474879/100‑1 (MQ=255)
                                                     tccatGACGCAACCCCATCAACGCGCAGCGCATCAATACCAACAAACTCAATCCAGTCACGCGCGTTACCGTCCAGGAAGTTACTGACTTCACGGCGAcc                                               <  2:486486/96‑1 (MQ=255)
                                                       aTTGACGCCACCGTATCGACGCGCAGCGGATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCAt                                             <  1:2644915/100‑1 (MQ=255)
                                                       aTTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCAt                                             >  1:235121/1‑100 (MQ=255)
                                                       aTTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCAt                                             >  1:2400785/1‑100 (MQ=255)
                                                       aTTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCAt                                             >  1:2432138/1‑100 (MQ=255)
                                                        ttGACGCCACCGCATCGTCGCGCAGTGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATa                                            <  1:1301546/100‑1 (MQ=255)
                                                        ttGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATa                                            <  2:1722084/100‑1 (MQ=255)
                                                        ttGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATa                                            <  1:1153343/100‑1 (MQ=255)
                                                          gACGCCACCGCATCGACGCGCAGCGCATCAATATCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGt                                          <  1:148127/100‑1 (MQ=255)
                                                          gACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGt                                          <  2:1179023/100‑1 (MQ=255)
                                                            cGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGCAGTTACTGACTTCACGGCGACCATAGTTg                                        >  2:2409615/1‑100 (MQ=255)
                                                            cGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTg                                        <  2:2561773/100‑1 (MQ=255)
                                                              ccACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTa                                      >  1:301210/1‑100 (MQ=255)
                                                              ccACCGCATCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTa                                      <  2:1256628/100‑1 (MQ=255)
                                                               cACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAg                                     <  1:2585354/100‑1 (MQ=255)
                                                               cACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAg                                     <  2:1152473/100‑1 (MQ=255)
                                                               aatcttaTCGACTCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCTCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAg                                     <  1:1887568/94‑1 (MQ=255)
                                                                aCCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGa                                    >  2:1153005/1‑100 (MQ=255)
                                                                 ccGCATCGTCGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAGGCGTGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGt                                       <  2:959056/96‑1 (MQ=255)
                                                                    naTCGACGCGCAGCGCATCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAg                                >  2:2262479/2‑100 (MQ=255)
                                                                      tCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCg                              <  1:791021/100‑1 (MQ=255)
                                                                          cgcgCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGtt                          >  2:307328/1‑100 (MQ=255)
                                                                              cagcgcaTCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAg                      <  2:1583473/100‑1 (MQ=255)
                                                                              cagcgcaTAAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAg                      <  1:1183803/100‑1 (MQ=255)
                                                                               agcgcaTCAATACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGt                     >  1:2277819/1‑100 (MQ=255)
                                                                                gcgcaTCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTc                    >  1:1388867/1‑100 (MQ=255)
                                                                                 cgcaTCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTGCCAGTcc                   <  2:500050/100‑1 (MQ=255)
                                                                                  gcaTCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCt                  <  2:601225/100‑1 (MQ=255)
                                                                                    aTCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtga                <  2:327566/100‑1 (MQ=255)
                                                                                    aTCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtga                <  2:235119/100‑1 (MQ=255)
                                                                                    aTCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtga                <  2:2400755/100‑1 (MQ=255)
                                                                                     tCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgat               >  1:1126365/1‑100 (MQ=255)
                                                                                      cAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCGTAGTTGTAGATCAGCGTGTTCCAGTCCtgatg              <  1:886601/100‑1 (MQ=255)
                                                                                      cAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgatg              >  2:266893/1‑100 (MQ=255)
                                                                                       aaTACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgatga             >  1:2269700/1‑100 (MQ=255)
                                                                                       aaTACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgatga             >  1:753454/1‑100 (MQ=255)
                                                                                       aaTACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCtgatga             >  1:671685/1‑100 (MQ=255)
                                                                                        aTACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTgatgat            >  1:2067805/1‑100 (MQ=255)
                                                                                         tACCAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATa           <  1:2275154/100‑1 (MQ=255)
                                                                                          aCCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAg          <  1:677349/100‑1 (MQ=255)
                                                                                           ccAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGc         >  1:970475/1‑100 (MQ=255)
                                                                                           ccAAAACGCTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCAAAGATGTAGATCAGCGAGTTCCAGTCCTGATGATAGc         >  1:871930/1‑100 (MQ=255)
                                                                                            cAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAg          >  2:1290226/1‑98 (MQ=255)
                                                                                            cAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGACGACCATAGTAGTAGATCAGCGTGTTCCAGTCCTGATGATAGcc        >  1:2135461/1‑100 (MQ=255)
                                                                                               aaCGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTc     <  2:968951/100‑1 (MQ=255)
                                                                                                  gTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCAcg  >  1:960269/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GCCAAATTCGTTCGGGATCCACTCCCCCTCTTTACGGCTGTAGTCGCGATAAATCATTGACGCCACCGCATCGACGCGCAGCGCATCAATACCAAAACGTTCAATCCAGTAAAGCGCGTTACCGACGAGGAAGTTACTGACTTCACGGCGACCATAGTTGTAGATCAGCGTGTTCCAGTCCTGATGATAGCCTTCACG  >  NC_000913/3572219‑3572416

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