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.

TGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTG  >  NC_000913/3570575‑3570873
                                                                                                                                                      |                                                                                                                                                     
  aTTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGgcagc                                                                                                                                                     >  2:1200014/1‑150 (MQ=255)
  aTTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGgcagc                                                                                                                                                     <  1:952748/150‑1 (MQ=255)
  aTTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGACGGATGCAGCGCAATCTTCGCCACCGGGTTGTAACCCCAGTAGTTTCTTAGCCCCTTGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGgcagc                                                                                                                                                     <  1:228122/150‑1 (MQ=255)
     gCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGt                                                                                                                                                  >  2:354512/1‑150 (MQ=255)
      cATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGC‑GCAGtt                                                                                                                                                 <  2:333426/150‑1 (MQ=255)
       aTCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTc                                                                                                                                                >  1:1560150/1‑150 (MQ=255)
        tCGCGAAACTCATCCAGCGCCTTTTCTGGCGAGCTGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAAGCCAAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTcc                                                                                                                                               <  2:437997/150‑1 (MQ=255)
        tCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGC‑GCAGtttg                                                                                                                                               <  1:1073305/150‑3 (MQ=255)
        tCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTcc                                                                                                                                               <  1:520066/150‑1 (MQ=255)
         cgcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGC‑GCAGtttgc                                                                                                                                              <  1:692464/150‑4 (MQ=255)
         cgcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCAGGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGC‑GCAGtttgc                                                                                                                                              <  1:1548870/150‑4 (MQ=255)
          gcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAg                                                                                                                                             <  1:1496201/150‑1 (MQ=255)
          gcgAAACTCATCCAGCGCCGTTTCT‑GCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAgc                                                                                                                                            >  1:1634681/1‑150 (MQ=255)
              aaCTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGcgcg                                                                                                                                         >  2:541614/1‑150 (MQ=255)
              aaCTCATCAAGCGCCGTTTCT‑GCGAGCAGGCATACGCCGGATGCAGCGGAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCACATGCGTTGCAGACGTGGTTCACAGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCgg                                                                                                                                        >  2:934176/1‑150 (MQ=255)
                    tCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTaaag                                                                                                                                   >  2:162968/1‑148 (MQ=255)
                      cAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGcc                                                                                                                                 >  2:77493/1‑150 (MQ=255)
                      cAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGCTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGcc                                                                                                                                 <  1:1105479/150‑1 (MQ=255)
                         cgcCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCaa                                                                                                                              >  2:1047059/1‑150 (MQ=255)
                            cGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTg                                                                                                                           >  1:183125/1‑150 (MQ=255)
                                tCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTc                                                                                                                       <  1:618236/150‑1 (MQ=255)
                                 cTGGCGAGCAGGCATACGCCGGATGCAGCGCCAACATCGCCACCGGGTTGTAACCCCAGTAGTTACATAGCCCCAAGCGTTGCAGACGTGGTTCACTGGCCAACTGC‑GCCACTGGCAGAAGGGCCCGCGCGGTAATGGCCAATGGTTTCa                                                                                                                      >  1:1197149/1‑150 (MQ=255)
                                 cTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCa                                                                                                                      <  2:646833/150‑1 (MQ=255)
                                   ggCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaaa                                                                                                                    >  1:1677864/1‑150 (MQ=255)
                                   ggCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaaa                                                                                                                    >  1:451340/1‑150 (MQ=255)
                                    gCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGAGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAAt                                                                                                                   <  1:1047059/150‑1 (MQ=255)
                                     cGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATa                                                                                                                  <  2:471437/150‑1 (MQ=255)
                                      gAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAg                                                                                                                 <  2:1692999/150‑1 (MQ=255)
                                       aGCAGGCATACGCCGGATGCAGCGCAGACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGt                                                                                                                >  2:275327/1‑150 (MQ=255)
                                        gCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTGACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGtt                                                                                                               <  2:128405/150‑1 (MQ=255)
                                        gCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGtt                                                                                                               <  2:1214332/150‑1 (MQ=255)
                                          aGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGa                                                                                                             <  2:26406/150‑1 (MQ=255)
                                           ggCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGat                                                                                                            >  2:1477229/1‑150 (MQ=255)
                                            gCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatc                                                                                                           <  1:1314651/150‑1 (MQ=255)
                                             cATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatca                                                                                                          <  1:1168215/150‑1 (MQ=255)
                                              aTACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatcat                                                                                                         >  1:571949/1‑150 (MQ=255)
                                               tACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatcatc                                                                                                        <  1:1015317/150‑1 (MQ=255)
                                                  gCCGGATGCAGCGCAAACATCGACACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCAcc                                                                                                     <  2:102922/150‑1 (MQ=255)
                                                   ccGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCg                                                                                                    <  2:1232309/150‑1 (MQ=255)
                                                     ggATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGa                                                                                                  <  2:828552/150‑1 (MQ=255)
                                                     ggATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGa                                                                                                  >  2:1616813/1‑150 (MQ=255)
                                                      gATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGAt                                                                                                 <  1:1646724/150‑1 (MQ=255)
                                                               cAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAggg                                                                                        <  1:628084/150‑1 (MQ=255)
                                                               cAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAggg                                                                                        <  1:270411/150‑1 (MQ=255)
                                                                     tCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGGGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCGTtat                                                                                  >  1:173815/1‑150 (MQ=255)
                                                                      cGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTtata                                                                                 <  1:984391/150‑1 (MQ=255)
                                                                       gCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATaa                                                                                >  2:756951/1‑150 (MQ=255)
                                                                        ccACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAg                                                                               <  1:227159/150‑1 (MQ=255)
                                                                         cACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGAGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGt                                                                              <  2:93702/150‑1 (MQ=255)
                                                                            cGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGcc                                                                           <  1:1616813/150‑1 (MQ=255)
                                                                                 tGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGa                                                                      >  1:776421/1‑150 (MQ=255)
                                                                                 tGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGa                                                                      <  1:275327/150‑1 (MQ=255)
                                                                                    aaCCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGAtct                                                                   >  1:1231501/1‑150 (MQ=255)
                                                                                      ccccAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCg                                                                 >  2:315370/1‑150 (MQ=255)
                                                                                       cccAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCTGCAGTTCCAGCGCGGTAATGCCAAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGa                                                                <  2:938102/150‑1 (MQ=255)
                                                                                       cccAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGa                                                                <  2:501069/150‑1 (MQ=255)
                                                                                             agtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTACAGCGCTGTAATGCCAAATTGTTTCAAATATTTTATCATCACCGGATGCACGAGGGCTTTATAATTGCCACGTATCTAGACCGCGa                                                          >  1:1061248/1‑150 (MQ=255)
                                                                                                tACTTAGCCCCATGCGTTGCAGACGTGGTTAACTGGCAAACTGC‑TCCACTGGCAGCATTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGAGCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCGCGACCGGGAtcg                                                       >  2:281564/1‑149 (MQ=255)
                                                                                                 aCTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCTCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGAtctc                                                      <  2:1150712/150‑1 (MQ=255)
                                                                                                  cTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATTGTTGATCATCACCGGATGCCCGAGGTCTTTATAAGTGCCACGGATCTCGACCGGGATCTcc                                                     <  1:989716/150‑1 (MQ=255)
                                                                                                  cTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTcc                                                     >  2:1023102/1‑150 (MQ=255)
                                                                                                    tAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCgg                                                   >  2:1564314/1‑150 (MQ=255)
                                                                                                    tAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCAACTGGCAGCAGTTCCATCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCgg                                                   >  1:1131574/1‑150 (MQ=255)
                                                                                                     cgCCCCATGCGATGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCACTTGCTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCggg                                                  <  1:224761/149‑1 (MQ=255)
                                                                                                     aGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCggg                                                  >  1:1715109/1‑150 (MQ=255)
                                                                                                      gCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGt                                                 <  1:756951/150‑1 (MQ=255)
                                                                                                           aTGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAg                                            <  2:1414115/150‑1 (MQ=255)
                                                                                                               gTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACg                                        >  2:1593424/1‑150 (MQ=255)
                                                                                                                  gCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTa                                     <  2:1560150/150‑1 (MQ=255)
                                                                                                                    aGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAAt                                   <  2:183125/150‑1 (MQ=255)
                                                                                                                      aCGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATcc                                 <  2:1392164/150‑1 (MQ=255)
                                                                                                                      aCGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATcc                                 >  1:1688360/1‑150 (MQ=255)
                                                                                                                        gTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATGGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCtt                               <  2:1139579/150‑1 (MQ=255)
                                                                                                                        gTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGGAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGACACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCtt                               >  2:1108578/1‑150 (MQ=255)
                                                                                                                         tGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCttt                              <  1:937324/150‑1 (MQ=255)
                                                                                                                          ggTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGGAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTg                             <  1:411535/150‑1 (MQ=255)
                                                                                                                           gTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGa                            >  1:487216/1‑150 (MQ=255)
                                                                                                                           gTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGa                            <  2:1324906/150‑1 (MQ=255)
                                                                                                                                cTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGa                            >  2:256112/1‑145 (MQ=255)
                                                                                                                                cTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGa                            <  1:256112/145‑1 (MQ=255)
                                                                                                                                 gctggaaCTGCTGCCAGTGGC‑GCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGCCATggg                      >  1:1092759/6‑150 (MQ=255)
                                                                                                                                  ggCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATTCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGC‑TTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCt                    <  1:541614/150‑1 (MQ=255)
                                                                                                                                  ggCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGc                     <  1:763470/150‑1 (MQ=255)
                                                                                                                                  ggCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGc                     >  2:1141501/1‑150 (MQ=255)
                                                                                                                                   gCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCt                    <  1:960594/150‑1 (MQ=255)
                                                                                                                                    cAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCtt                   <  2:822201/150‑1 (MQ=255)
                                                                                                                                     aaaCTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTc                  <  2:1509761/150‑1 (MQ=255)
                                                                                                                                      aaCTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCa                 >  1:1522299/1‑150 (MQ=255)
                                                                                                                                       actgctgcCAGTGGC‑GCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCAt                >  2:97706/7‑150 (MQ=255)
                                                                                                                                       aCTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCAt                <  1:760571/150‑1 (MQ=255)
                                                                                                                                       aCTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCAt                >  1:452561/1‑150 (MQ=255)
                                                                                                                                       aCTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCAt                <  1:1484894/150‑1 (MQ=255)
                                                                                                                                         tGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTTATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCCCGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATaa              >  1:1245381/1‑150 (MQ=255)
                                                                                                                                         tGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATaa              <  2:1029904/150‑1 (MQ=255)
                                                                                                                                          gc‑gcCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATaaa             <  1:929265/150‑1 (MQ=255)
                                                                                                                                           c‑gcCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAat            <  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)
                                                                                                                                                      |                                                                                                                                                     
TGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCTG  >  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: ‑