Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,570,724 Δ1 bp 8.4% coding (596/1974 nt) glgX ← limit dextrin alpha‑1,6‑glucohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,570,7240A. 8.4% 305.7 / 20.2 95coding (596/1974 nt)glgXlimit dextrin alpha‑1,6‑glucohydrolase
Reads supporting (aligned to +/- strand):  ref base A (36/51);  new base . (6/2);  total (42/53)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 101 total aligned reads displayed.

TGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCT  >  NC_000913/3570575‑3570872
                                                                                                                                                      |                                                                                                                                                    
tGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgca                                                                                                                                                      <  2:1143901/150‑1 (MQ=255)
tGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgca                                                                                                                                                      <  2:357539/150‑1 (MQ=255)
 gATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgcag                                                                                                                                                     >  1:1245510/1‑150 (MQ=255)
 gATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgcag                                                                                                                                                     >  1:1240177/1‑150 (MQ=255)
      cATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGtt                                                                                                                                                >  2:1513608/1‑150 (MQ=255)
        tCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTcc                                                                                                                                              >  1:775989/1‑150 (MQ=255)
        tCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTcc                                                                                                                                              <  2:753711/150‑1 (MQ=255)
         cgcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGCGCAGtttgc                                                                                                                                             <  2:1519716/150‑4 (MQ=255)
         cgcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCa                                                                                                                                             >  1:1306252/1‑150 (MQ=255)
         cgcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCa                                                                                                                                             >  2:563114/1‑150 (MQ=255)
          tcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAg                                                                                                                                            >  1:998441/2‑150 (MQ=255)
          gcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCACATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGCGCAGtttgcc                                                                                                                                            <  2:661291/150‑5 (MQ=255)
           cgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAgc                                                                                                                                           >  1:862757/1‑150 (MQ=255)
              aaCTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGcgcg                                                                                                                                        <  1:1261047/150‑1 (MQ=255)
               aCTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGGTGTAACCCCAGTATTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCgg                                                                                                                                       <  1:1173825/150‑1 (MQ=255)
                  cATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTaa                                                                                                                                    >  1:1184454/1‑150 (MQ=255)
                     ccAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGGTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGc                                                                                                                                 >  1:260016/1‑150 (MQ=255)
                       aGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGccc                                                                                                                               <  1:611793/150‑1 (MQ=255)
                         cgcCGTTTCCGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCaa                                                                                                                             <  2:46842/150‑1 (MQ=255)
                            cGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTg                                                                                                                          <  2:502786/150‑1 (MQ=255)
                            cGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTg                                                                                                                          >  2:851238/1‑150 (MQ=255)
                               ttCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGttt                                                                                                                       <  2:220747/150‑1 (MQ=255)
                               ttCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGttt                                                                                                                       <  2:1147864/150‑1 (MQ=255)
                                tCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTc                                                                                                                      <  2:539703/150‑1 (MQ=255)
                                tCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTc                                                                                                                      <  2:751625/150‑1 (MQ=255)
                                 cTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCa                                                                                                                     <  2:174983/150‑1 (MQ=255)
                                 cTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCa                                                                                                                     <  1:94844/150‑1 (MQ=255)
                                 cTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCa                                                                                                                     >  2:94844/1‑150 (MQ=255)
                                  tGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaa                                                                                                                    <  2:965593/150‑1 (MQ=255)
                                   ggCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaaa                                                                                                                   <  1:410204/150‑1 (MQ=255)
                                    gCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAAt                                                                                                                  <  1:435779/150‑1 (MQ=255)
                                    gCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAAt                                                                                                                  <  2:831737/150‑1 (MQ=255)
                                            gCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGGTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatc                                                                                                          >  1:613511/1‑150 (MQ=255)
                                             cATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCAAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatca                                                                                                         >  2:1261998/1‑150 (MQ=255)
                                               tACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatcatc                                                                                                       >  1:16391/1‑150 (MQ=255)
                                               tACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAATAGTTGAtcatca                                                                                                      >  2:62314/1‑150 (MQ=255)
                                                aCGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGAtcatca                                                                                                      >  1:990795/1‑150 (MQ=255)
                                                 cGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCAc                                                                                                     >  1:1027969/1‑150 (MQ=255)
                                                  gCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCCCTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCAcc                                                                                                    >  2:1310968/1‑150 (MQ=255)
                                                    ctgATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCgg                                                                                                  <  2:828277/148‑1 (MQ=255)
                                                     ggATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGa                                                                                                 <  1:670387/150‑1 (MQ=255)
                                                       aTGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATg                                                                                               <  2:421675/150‑1 (MQ=255)
                                                        tGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGc                                                                                              <  1:210537/150‑1 (MQ=255)
                                                        tGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGc                                                                                              >  1:583371/1‑150 (MQ=255)
                                                             cgcAAACATCGCCACGGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAg                                                                                         <  2:74465/150‑1 (MQ=255)
                                                             cgcAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAg                                                                                         <  1:574029/150‑1 (MQ=255)
                                                               cAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAggg                                                                                       <  1:949486/150‑1 (MQ=255)
                                                                 aaCATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAgg                                                                                        >  1:844537/1‑147 (MQ=255)
                                                                 aaCATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAgg                                                                                        <  2:844537/147‑1 (MQ=255)
                                                                  aCATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCtt                                                                                    <  2:16391/150‑1 (MQ=255)
                                                                   cATCGCCACCGGGTTGTAAGCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGGTCACTGGCAACCTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCATTTGTTTCAAATAGTTGATCATGACCGGATGCCCGAGGGCttt                                                                                   <  2:428659/150‑1 (MQ=255)
                                                                      cGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTtata                                                                                <  1:1443632/150‑1 (MQ=255)
                                                                         cACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCCGATGCCCGAGGGCTTTATAAGt                                                                             >  2:1485558/1‑150 (MQ=255)
                                                                          aCCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTg                                                                            <  2:750152/150‑1 (MQ=255)
                                                                          aCCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTg                                                                            <  1:849563/150‑1 (MQ=255)
                                                                                  gTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGAt                                                                    >  1:95241/1‑150 (MQ=255)
                                                                                   tAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGAtc                                                                   <  1:1513608/150‑1 (MQ=255)
                                                                                    aaCCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGAtct                                                                  <  2:1240177/150‑1 (MQ=255)
                                                                                    aaCCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGAtct                                                                  <  2:1245510/150‑1 (MQ=255)
                                                                                       cccAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGa                                                               >  1:791584/1‑150 (MQ=255)
                                                                                            tagtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCggg                                                          <  2:260016/150‑1 (MQ=255)
                                                                                            tagtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCggg                                                          <  1:1528389/150‑1 (MQ=255)
                                                                                             agtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGa                                                         >  1:258429/1‑150 (MQ=255)
                                                                                                tACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGAtct                                                      >  1:801205/1‑150 (MQ=255)
                                                                                                   ttAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGGGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGCTTTATAAGTGCCACGGATCCCGGTCGAGAtccgtgg                                                  <  2:896199/150‑6 (MQ=255)
                                                                                                     aGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCggg                                                 <  2:1431573/150‑1 (MQ=255)
                                                                                                       ccccATGCGTTGCAGACGTGGTTCCCTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTg                                               <  1:1310968/150‑1 (MQ=255)
                                                                                                       ccccATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTg                                               >  1:420987/1‑150 (MQ=255)
                                                                                                       ccccATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTg                                               <  1:851238/150‑1 (MQ=255)
                                                                                                           aTGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAg                                           <  2:862757/150‑1 (MQ=255)
                                                                                                             gCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTa                                         <  1:602194/150‑1 (MQ=255)
                                                                                                               gTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACg                                       >  2:1448516/1‑150 (MQ=255)
                                                                                                               gTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACg                                       <  2:550102/150‑1 (MQ=255)
                                                                                                               gTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCAGGGTGCAAGTACg                                       <  2:1306252/150‑1 (MQ=255)
                                                                                                                 tGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGtt                                     <  1:915481/150‑1 (MQ=255)
                                                                                                                          ggTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTg                            <  1:1111023/150‑1 (MQ=255)
                                                                                                                            ttCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGAc                          >  1:1464205/1‑150 (MQ=255)
                                                                                                                              cACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACAt                        >  2:1249166/1‑150 (MQ=255)
                                                                                                                                cTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGCCCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATgg                      >  1:1006597/1‑150 (MQ=255)
                                                                                                                                 tGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATggg                     <  1:1045720/150‑1 (MQ=255)
                                                                                                                                 gctggaaCTGCTGCCAGTGGCGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATggg                     >  2:488290/6‑150 (MQ=255)
                                                                                                                                 gctggaaCTGCTGCCAGTGGCGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATggg                     >  2:316037/6‑150 (MQ=255)
                                                                                                                                 gctggaaCTGCTGCCAGTGGCGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATggg                     >  1:1297539/6‑150 (MQ=255)
                                                                                                                                  ggCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTg                            <  2:1058097/142‑1 (MQ=255)
                                                                                                                                  ggCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTg                            >  1:1058097/1‑142 (MQ=255)
                                                                                                                                  ggCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGc                    >  1:211063/1‑150 (MQ=255)
                                                                                                                                  ggCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGc                    <  1:1249166/150‑1 (MQ=255)
                                                                                                                                    ggaaCTGCTGCCAGTGGCGCAGTTCCCGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCtt                  >  1:1436367/3‑150 (MQ=255)
                                                                                                                                    ggaaCTGCTGCCAGTGGCGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCtt                  >  2:28677/3‑150 (MQ=255)
                                                                                                                                    cAAACTGCGCCACTGGCAGCAGTTCCAGCGTGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCtt                  >  1:1259447/1‑150 (MQ=255)
                                                                                                                                      aaCTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCa                <  1:1448516/150‑1 (MQ=255)
                                                                                                                                        cTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGCCCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATa              >  2:421559/1‑150 (MQ=255)
                                                                                                                                        cTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATa              >  1:639286/1‑150 (MQ=255)
                                                                                                                                        cTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATa              <  2:1027969/150‑1 (MQ=255)
                                                                                                                                         tgctgcCAGTGGCGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACCTTAATCATTTGACATCGGCTTCATaa             >  2:1197010/5‑150 (MQ=255)
                                                                                                                                              ccACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCGGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATTTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAatga         <  2:1184454/150‑1 (MQ=255)
                                                                                                                                                 cTGTCAGCAGTTCGAGCGCGGTAATCCCCATTTGTGTCAAATAGTTGATCATCAGCGGAGGCCCGAGGGCTTTATAAGTGCCTCGTATCTCGACCGGGTTCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATTAATGATgg      <  2:1528318/150‑1 (MQ=255)
                                                                                                                                                   ggCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGtg    >  2:612114/1‑150 (MQ=255)
                                                                                                                                                     cagcagTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGtgct  >  1:479750/1‑150 (MQ=255)
                                                                                                                                                     cagcagTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGtcct  >  2:306746/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCT  >  NC_000913/3570575‑3570872

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