Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,570,724 Δ1 bp 5.6% coding (596/1974 nt) glgX ← glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,570,7240A. 5.6% 337.2 / 11.8 108coding (596/1974 nt)glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base A (45/56);  new base . (2/4);  total (47/61)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

Only 100 of 113 total aligned reads displayed.

TGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTG  >  NC_000913/3570575‑3570873
                                                                                                                                                      |                                                                                                                                                     
  aTTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgcagc                                                                                                                                                     >  2:1200014/1‑150 (MQ=255)
  aTTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgcagc                                                                                                                                                     <  1:952748/150‑1 (MQ=255)
  aTTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGACGGATGCAGCGCAATCTTCGCCACCGGGTTGTAACCCCAGTAGTTTCTTAGCCCCTTGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgcagc                                                                                                                                                     <  1:228122/150‑1 (MQ=255)
     gCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGt                                                                                                                                                  >  2:354512/1‑150 (MQ=255)
      cATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGCGCAGtt                                                                                                                                                 <  2:333426/150‑1 (MQ=255)
       aTCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTc                                                                                                                                                >  1:1560150/1‑150 (MQ=255)
        tCGCGAAACTCATCCAGCGCCTTTTCTGGCGAGCTGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAAGCCAAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTcc                                                                                                                                               <  2:437997/150‑1 (MQ=255)
        tCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGCGCAGtttg                                                                                                                                               <  1:1073305/150‑3 (MQ=255)
        tCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTcc                                                                                                                                               <  1:520066/150‑1 (MQ=255)
         cgcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGCGCAGtttgc                                                                                                                                              <  1:692464/150‑4 (MQ=255)
         cgcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCAGGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGCGCAGtttgc                                                                                                                                              <  1:1548870/150‑4 (MQ=255)
          gcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAg                                                                                                                                             <  1:1496201/150‑1 (MQ=255)
          gcgAAACTCATCCAGCGCCGTTTCTGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAgc                                                                                                                                            >  1:1634681/1‑150 (MQ=255)
              aaCTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGcgcg                                                                                                                                         >  2:541614/1‑150 (MQ=255)
              aaCTCATCAAGCGCCGTTTCTGCGAGCAGGCATACGCCGGATGCAGCGGAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCACATGCGTTGCAGACGTGGTTCACAGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCgg                                                                                                                                        >  2:934176/1‑150 (MQ=255)
                    tCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTaaag                                                                                                                                   >  2:162968/1‑148 (MQ=255)
                      cAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGcc                                                                                                                                 >  2:77493/1‑150 (MQ=255)
                      cAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGCTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGcc                                                                                                                                 <  1:1105479/150‑1 (MQ=255)
                         cgcCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCaa                                                                                                                              >  2:1047059/1‑150 (MQ=255)
                            cGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTg                                                                                                                           >  1:183125/1‑150 (MQ=255)
                                tCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTc                                                                                                                       <  1:618236/150‑1 (MQ=255)
                                 cTGGCGAGCAGGCATACGCCGGATGCAGCGCCAACATCGCCACCGGGTTGTAACCCCAGTAGTTACATAGCCCCAAGCGTTGCAGACGTGGTTCACTGGCCAACTGCGCCACTGGCAGAAGGGCCCGCGCGGTAATGGCCAATGGTTTCa                                                                                                                      >  1:1197149/1‑150 (MQ=255)
                                 cTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCa                                                                                                                      <  2:646833/150‑1 (MQ=255)
                                   ggCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaaa                                                                                                                    >  1:1677864/1‑150 (MQ=255)
                                   ggCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaaa                                                                                                                    >  1:451340/1‑150 (MQ=255)
                                    gCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGAGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAAt                                                                                                                   <  1:1047059/150‑1 (MQ=255)
                                     cGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATa                                                                                                                  <  2:471437/150‑1 (MQ=255)
                                      gAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAg                                                                                                                 <  2:1692999/150‑1 (MQ=255)
                                       aGCAGGCATACGCCGGATGCAGCGCAGACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGt                                                                                                                >  2:275327/1‑150 (MQ=255)
                                        gCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTGACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGtt                                                                                                               <  2:128405/150‑1 (MQ=255)
                                        gCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGtt                                                                                                               <  2:1214332/150‑1 (MQ=255)
                                          aGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGa                                                                                                             <  2:26406/150‑1 (MQ=255)
                                           ggCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGat                                                                                                            >  2:1477229/1‑150 (MQ=255)
                                            gCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatc                                                                                                           <  1:1314651/150‑1 (MQ=255)
                                             cATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatca                                                                                                          <  1:1168215/150‑1 (MQ=255)
                                              aTACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatcat                                                                                                         >  1:571949/1‑150 (MQ=255)
                                               tACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatcatc                                                                                                        <  1:1015317/150‑1 (MQ=255)
                                                  gCCGGATGCAGCGCAAACATCGACACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCAcc                                                                                                     <  2:102922/150‑1 (MQ=255)
                                                   ccGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCg                                                                                                    <  2:1232309/150‑1 (MQ=255)
                                                     ggATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGa                                                                                                  <  2:828552/150‑1 (MQ=255)
                                                     ggATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGa                                                                                                  >  2:1616813/1‑150 (MQ=255)
                                                      gATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGAt                                                                                                 <  1:1646724/150‑1 (MQ=255)
                                                               cAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAggg                                                                                        <  1:628084/150‑1 (MQ=255)
                                                               cAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAggg                                                                                        <  1:270411/150‑1 (MQ=255)
                                                                     tCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGGGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCGTtat                                                                                  >  1:173815/1‑150 (MQ=255)
                                                                      cGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTtata                                                                                 <  1:984391/150‑1 (MQ=255)
                                                                       gCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATaa                                                                                >  2:756951/1‑150 (MQ=255)
                                                                        ccACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAg                                                                               <  1:227159/150‑1 (MQ=255)
                                                                         cACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGAGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGt                                                                              <  2:93702/150‑1 (MQ=255)
                                                                            cGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGcc                                                                           <  1:1616813/150‑1 (MQ=255)
                                                                                 tGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGa                                                                      >  1:776421/1‑150 (MQ=255)
                                                                                 tGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGa                                                                      <  1:275327/150‑1 (MQ=255)
                                                                                    aaCCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGAtct                                                                   >  1:1231501/1‑150 (MQ=255)
                                                                                      ccccAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCg                                                                 >  2:315370/1‑150 (MQ=255)
                                                                                       cccAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCTGCAGTTCCAGCGCGGTAATGCCAAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGa                                                                <  2:938102/150‑1 (MQ=255)
                                                                                       cccAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGa                                                                <  2:501069/150‑1 (MQ=255)
                                                                                             agtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTACAGCGCTGTAATGCCAAATTGTTTCAAATATTTTATCATCACCGGATGCACGAGGGCTTTATAATTGCCACGTATCTAGACCGCGa                                                          >  1:1061248/1‑150 (MQ=255)
                                                                                                tACTTAGCCCCATGCGTTGCAGACGTGGTTAACTGGCAAACTGCTCCACTGGCAGCATTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGAGCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCGCGACCGGGAtcg                                                       >  2:281564/1‑149 (MQ=255)
                                                                                                 aCTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCTCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGAtctc                                                      <  2:1150712/150‑1 (MQ=255)
                                                                                                  cTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATTGTTGATCATCACCGGATGCCCGAGGTCTTTATAAGTGCCACGGATCTCGACCGGGATCTcc                                                     <  1:989716/150‑1 (MQ=255)
                                                                                                  cTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTcc                                                     >  2:1023102/1‑150 (MQ=255)
                                                                                                    tAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCgg                                                   >  2:1564314/1‑150 (MQ=255)
                                                                                                    tAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCAACTGGCAGCAGTTCCATCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCgg                                                   >  1:1131574/1‑150 (MQ=255)
                                                                                                     cgCCCCATGCGATGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCACTTGCTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCggg                                                  <  1:224761/149‑1 (MQ=255)
                                                                                                     aGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCggg                                                  >  1:1715109/1‑150 (MQ=255)
                                                                                                      gCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGt                                                 <  1:756951/150‑1 (MQ=255)
                                                                                                           aTGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAg                                            <  2:1414115/150‑1 (MQ=255)
                                                                                                               gTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACg                                        >  2:1593424/1‑150 (MQ=255)
                                                                                                                  gCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTa                                     <  2:1560150/150‑1 (MQ=255)
                                                                                                                    aGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAAt                                   <  2:183125/150‑1 (MQ=255)
                                                                                                                      aCGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATcc                                 <  2:1392164/150‑1 (MQ=255)
                                                                                                                      aCGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATcc                                 >  1:1688360/1‑150 (MQ=255)
                                                                                                                        gTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATGGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCtt                               <  2:1139579/150‑1 (MQ=255)
                                                                                                                        gTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGGAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGACACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCtt                               >  2:1108578/1‑150 (MQ=255)
                                                                                                                         tGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCttt                              <  1:937324/150‑1 (MQ=255)
                                                                                                                          ggTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGGAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTg                             <  1:411535/150‑1 (MQ=255)
                                                                                                                           gTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGa                            >  1:487216/1‑150 (MQ=255)
                                                                                                                           gTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGa                            <  2:1324906/150‑1 (MQ=255)
                                                                                                                                cTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGa                            >  2:256112/1‑145 (MQ=255)
                                                                                                                                cTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGa                            <  1:256112/145‑1 (MQ=255)
                                                                                                                                 gctggaaCTGCTGCCAGTGGCGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGCCATggg                      >  1:1092759/6‑150 (MQ=255)
                                                                                                                                  ggCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATTCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCt                    <  1:541614/150‑1 (MQ=255)
                                                                                                                                  ggCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGc                     <  1:763470/150‑1 (MQ=255)
                                                                                                                                  ggCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGc                     >  2:1141501/1‑150 (MQ=255)
                                                                                                                                   gCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCt                    <  1:960594/150‑1 (MQ=255)
                                                                                                                                    cAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCtt                   <  2:822201/150‑1 (MQ=255)
                                                                                                                                     aaaCTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTc                  <  2:1509761/150‑1 (MQ=255)
                                                                                                                                      aaCTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCa                 >  1:1522299/1‑150 (MQ=255)
                                                                                                                                       actgctgcCAGTGGCGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCAt                >  2:97706/7‑150 (MQ=255)
                                                                                                                                       aCTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCAt                <  1:760571/150‑1 (MQ=255)
                                                                                                                                       aCTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCAt                >  1:452561/1‑150 (MQ=255)
                                                                                                                                       aCTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCAt                <  1:1484894/150‑1 (MQ=255)
                                                                                                                                         tGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTTATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCCCGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATaa              >  1:1245381/1‑150 (MQ=255)
                                                                                                                                         tGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATaa              <  2:1029904/150‑1 (MQ=255)
                                                                                                                                          gcgcCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATaaa             <  1:929265/150‑1 (MQ=255)
                                                                                                                                           cgcCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAat            <  1:1156474/150‑1 (MQ=255)
                                                                                                                                              ccACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAatga          <  2:1121246/150‑1 (MQ=255)
                                                                                                                                              ccACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAatga          >  1:1149263/1‑150 (MQ=255)
                                                                                                                                                   ggCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGtg     >  2:603008/1‑150 (MQ=255)
                                                                                                                                                      agcagTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGtgctg  >  2:620757/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTG  >  NC_000913/3570575‑3570873

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: