Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,571,851 C→T R551H (CGC→CAC)  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,571,8510CT99.2% 497.7 / ‑4.9 130R551H (CGC→CAC) glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (65/64);  total (66/64)
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.90e-01

Only 100 of 131 total aligned reads displayed.

TGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATG  >  NC_000913/3571704‑3571996
                                                                                                                                                   |                                                                                                                                                 
  gACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAgg                                                                                                                                              <  2:752993/151‑1 (MQ=255)
  gACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTCTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAgg                                                                                                                                              <  2:374359/151‑1 (MQ=255)
        gtggtgCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGATTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTcgcg                                                                                                                                        <  1:1533152/151‑1 (MQ=255)
         tggtgCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTcgcg                                                                                                                                        >  2:35030/1‑150 (MQ=255)
         tggtgCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTcgcg                                                                                                                                        >  3:149797/1‑150 (MQ=255)
           gtgCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAAt                                                                                                                                     >  2:848687/1‑151 (MQ=255)
              ccAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTc                                                                                                                                  <  1:1204235/151‑1 (MQ=255)
                    aTCGCGGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCAt                                                                                                                            >  1:1161698/1‑151 (MQ=255)
                    aTCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCAt                                                                                                                            >  3:316127/1‑151 (MQ=255)
                      cgccgcCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATg                                                                                                                           <  2:990731/150‑1 (MQ=255)
                         cgccTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTc                                                                                                                       >  4:597555/1‑151 (MQ=255)
                          gccTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCg                                                                                                                      <  2:655719/151‑1 (MQ=255)
                          gccTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCg                                                                                                                      >  2:164091/1‑151 (MQ=255)
                            cTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGAAGGCGTGCAGGTTCGCGCATTTCTGCCATGCGTc                                                                                                                       >  1:1152357/1‑148 (MQ=255)
                              tCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGccc                                                                                                                   <  3:199177/150‑1 (MQ=255)
                               ccAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTCCCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGccc                                                                                                                   <  3:290490/149‑1 (MQ=255)
                               ccAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCgg                                                                                                                 >  1:1577002/1‑151 (MQ=255)
                                    aGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGc                                                                                                            <  1:1323494/151‑1 (MQ=255)
                                      aTGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCgg                                                                                                          <  2:1692317/151‑1 (MQ=255)
                                        gCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCGCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTc                                                                                                        <  1:418295/151‑1 (MQ=255)
                                        gCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGt                                                                                                         <  2:592016/150‑1 (MQ=255)
                                        gCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTc                                                                                                        <  1:777279/151‑1 (MQ=255)
                                              cGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAAt                                                                                                  >  1:976918/1‑151 (MQ=255)
                                               gAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATc                                                                                                 >  2:1734226/1‑151 (MQ=255)
                                               gAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATc                                                                                                 <  2:117987/151‑1 (MQ=255)
                                                     ggCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGAttt                                                                                            <  2:1208193/150‑1 (MQ=255)
                                                       cGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGAtttttt                                                                                         <  2:1502348/151‑1 (MQ=255)
                                                          cATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTAc                                                                                       <  1:566207/150‑1 (MQ=255)
                                                              gTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTgg                                                                                  <  1:946321/151‑1 (MQ=255)
                                                              gTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTgg                                                                                  <  2:1484233/151‑1 (MQ=255)
                                                               ttCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTCCCGTgg                                                                                  <  1:341394/150‑1 (MQ=255)
                                                                  cACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGAcc                                                                               >  2:1428762/1‑150 (MQ=255)
                                                                  cACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCa                                                                              <  2:1610799/151‑1 (MQ=255)
                                                                   aCTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCAc                                                                             >  2:1353341/1‑151 (MQ=255)
                                                                         ggCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTtcatc                                                                       <  2:1524408/151‑1 (MQ=255)
                                                                            ccTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTcatcat                                                                     <  2:1403308/150‑1 (MQ=255)
                                                                            ccTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATg                                                                    <  1:133516/151‑1 (MQ=255)
                                                                             cTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTc                                                                                                                       >  2:1183306/1‑99 (MQ=255)
                                                                              tGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCg                                                                  >  1:961266/1‑151 (MQ=255)
                                                                              tGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCg                                                                  >  2:607896/1‑151 (MQ=255)
                                                                              tGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCg                                                                  >  2:1147210/1‑151 (MQ=255)
                                                                                  cAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACa                                                               >  2:939953/1‑150 (MQ=255)
                                                                                   aaaTTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCg                                                                                                       >  3:44689/1‑109 (MQ=255)
                                                                                   aaaTTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCg                                                                                                       <  4:44688/109‑1 (MQ=255)
                                                                                   aaaTTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAAc                                                             <  1:1144316/151‑1 (MQ=255)
                                                                                    aaTTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACg                                                            >  1:1250682/1‑151 (MQ=255)
                                                                                     aTTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACgg                                                           <  1:208563/151‑1 (MQ=255)
                                                                                     aTTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACgg                                                           >  2:1189508/1‑151 (MQ=255)
                                                                                     aTTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACgg                                                           >  2:328892/1‑151 (MQ=255)
                                                                                      ttCGCTACCCATGTACAGTAGTTTCTTGCTGGGCAATGCCGACATTAAGCCATAGTAAGCGCGCAGGTTCGTTGATTTCTGCCATGCGACGCTCGGCGTGAGGTAGCGAATGGAATTTTTACCGTGGACCACTTCATTATTCGACAACAGc                                                          >  1:1205614/1‑151 (MQ=255)
                                                                                       tCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGc                                                          >  2:1738191/1‑150 (MQ=255)
                                                                                        cGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCa                                                         <  2:1380171/150‑1 (MQ=255)
                                                                                          ttACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACAGCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAgg                                                       <  1:1668597/150‑1 (MQ=255)
                                                                                            aCCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGAc                                                     <  2:605425/150‑1 (MQ=255)
                                                                                            aCCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACg                                                    >  1:1196519/1‑151 (MQ=255)
                                                                                             cccATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACg                                                    >  1:171966/1‑150 (MQ=255)
                                                                                              ccATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGaa                                                  <  2:89590/151‑1 (MQ=255)
                                                                                               cATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGa                                                   >  2:680007/1‑149 (MQ=255)
                                                                                                aTGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGt                                                >  2:1559544/1‑151 (MQ=255)
                                                                                                aTGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGt                                                <  2:1278437/151‑1 (MQ=255)
                                                                                                  gAACAGTAGTTTCTTGCCCGGGAATGCCCGCATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGtt                                               <  1:1705893/150‑1 (MQ=255)
                                                                                                  gAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGtt                                               <  3:289217/150‑1 (MQ=255)
                                                                                                   aaCAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGtttt                                             >  2:677825/1‑151 (MQ=255)
                                                                                                    aCAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTc                                            >  2:276670/1‑151 (MQ=255)
                                                                                                        tagtTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAgtg                                        >  1:609599/1‑151 (MQ=255)
                                                                                                        tagtTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAgt                                         >  4:499501/1‑150 (MQ=255)
                                                                                                         agtTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAgtgt                                       >  4:580980/1‑151 (MQ=255)
                                                                                                         agtTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAgtgt                                       >  2:383352/1‑151 (MQ=255)
                                                                                                            ttCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTCCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGt                                    >  2:1449091/1‑151 (MQ=255)
                                                                                                                 gCCCGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTag                               <  2:770665/151‑1 (MQ=255)
                                                                                                                  cccGGGAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGt                                 >  2:813296/1‑148 (MQ=255)
                                                                                                                     gggAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAgaga                            >  1:778644/1‑150 (MQ=255)
                                                                                                                     gggAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGaa                           <  1:979511/151‑1 (MQ=255)
                                                                                                                      ggAATGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGaa                           >  4:297392/1‑150 (MQ=255)
                                                                                                                        aaTGGCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCTATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAgaga                            <  1:39411/147‑1 (MQ=255)
                                                                                                                        aaTGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATcc                        <  2:139580/151‑1 (MQ=255)
                                                                                                                        aaTGCCCACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATcc                        >  1:386598/1‑151 (MQ=255)
                                                                                                                            cccACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGaa                    <  4:299449/151‑1 (MQ=255)
                                                                                                                            cccACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGa                     >  1:330310/1‑150 (MQ=255)
                                                                                                                             ccACATCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAg                   <  1:472258/151‑1 (MQ=255)
                                                                                                                              catatccatccCTAGTCGGCGTGCAGGTTCGCGAATTTCTGACATGCGTCGCCCGGCATGCGGTCGAGCATCGATTTTTTACCGTTGACCACTTCATCATGCCAAAACTGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATTTCGaa                    <  1:1047168/146‑1 (MQ=255)
                                                                                                                              cacaTCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAgg                  <  1:558398/151‑1 (MQ=255)
                                                                                                                              cacaTCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAgg                  <  3:534084/151‑1 (MQ=255)
                                                                                                                              cacaTCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAgg                  >  1:405744/1‑151 (MQ=255)
                                                                                                                               acaTCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGt                 >  1:1649524/1‑151 (MQ=255)
                                                                                                                                caTCCAGCCATAGTAGGCGTGCAGGTTCGCTAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTc                <  1:1003328/151‑1 (MQ=255)
                                                                                                                                caTCCAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTc                <  1:421382/151‑1 (MQ=255)
                                                                                                                                   ccAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGTAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAg              <  1:816154/150‑1 (MQ=255)
                                                                                                                                   ccAGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAg              <  1:785788/150‑1 (MQ=255)
                                                                                                                                     aGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGttt           <  1:1615682/151‑1 (MQ=255)
                                                                                                                                     aGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGttt           <  1:1086923/151‑1 (MQ=255)
                                                                                                                                     aGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGttt           <  1:173943/151‑1 (MQ=255)
                                                                                                                                     aGCCATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGttt           >  2:749794/1‑151 (MQ=255)
                                                                                                                                       ccATAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTAt         <  2:1594101/151‑1 (MQ=255)
                                                                                                                                         aTAGTAGGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCg       >  1:1012377/1‑151 (MQ=255)
                                                                                                                                            gtagGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGtga    <  2:340645/151‑1 (MQ=255)
                                                                                                                                            gtagGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGtg     >  1:938380/1‑150 (MQ=255)
                                                                                                                                             tagGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGtgat   >  2:1666684/1‑151 (MQ=255)
                                                                                                                                             tagGCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGtgat   >  2:559940/1‑151 (MQ=255)
                                                                                                                                               ggCGTGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGtgatg  >  3:410147/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGGACACCGTGGTGCCAGTTATCGCCGCCTTCCAACAGATGCCAGTCGAGGCTGGCGTCATGGTTCCACTCGCGGCCCTGGGCAAATTCGTTACCCATGAACAGTAGTTTCTTGCCCGGGAATGCCCACATCCAGCCATAGTAGGCGCGCAGGTTCGCGAATTTCTGCCATGCGTCGCCCGGCATGCGGTCGAGAATCGATTTTTTACCGTGGACCACTTCATCATGCGACAACGGCAGGACGAAGTTTTCAGTGTAGTTGTAGAGAATCCCGAAGGTCAGTTTATCGTGATG  >  NC_000913/3571704‑3571996

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