Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,570,724 Δ1 bp 5.8% 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.8% 236.7 / 9.6 86coding (596/1974 nt)glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base A (40/41);  new base . (3/2);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

GATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATG  >  NC_000913/3570576‑3570867
                                                                                                                                                     |                                                                                                                                               
gATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGgca                                                                                                                                                 >  1:152116/1‑149 (MQ=255)
gATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGgca                                                                                                                                                 >  2:739772/1‑149 (MQ=255)
 aTTGCATCGCGAAACTCATCCCGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGGTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTTGCAAACTGC‑GCCACTGgcag                                                                                                                                                <  2:466704/149‑1 (MQ=255)
  ttGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCAACCGGGTTGTAACCCCAGTAGTTACTTAGACCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGgcagc                                                                                                                                               >  1:572915/1‑149 (MQ=255)
     cATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGt                                                                                                                                            <  1:4634/149‑1 (MQ=255)
      aTCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACAACAGTAGTTTCTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGtt                                                                                                                                           >  2:17763/1‑149 (MQ=255)
      aTCGCGAAACTCATCCAGCGCCGGTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGtt                                                                                                                                           <  1:273686/149‑1 (MQ=255)
       tCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGC‑GCAGttt                                                                                                                                          <  2:287208/149‑2 (MQ=255)
       tCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTc                                                                                                                                          >  2:1171762/1‑149 (MQ=255)
          cgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGGGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAg                                                                                                                                       >  1:805179/1‑149 (MQ=255)
           gAAACTCATACAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGATGCAGACGTGGTTCACTGGCAATCTGC‑GCCACTGGCAGCAGTTCCAgc                                                                                                                                      >  2:19512/1‑149 (MQ=255)
             aaCTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGC‑GCAGtttgccagt                                                                                                                                    <  2:1119194/149‑8 (MQ=255)
              aCTCATCCAGCGCCGTTTCTGGCGAGCAGGTATCCGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGAGGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGcgcg                                                                                                                                   >  1:910050/1‑149 (MQ=255)
              aCTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGcgcg                                                                                                                                   <  2:1061158/149‑1 (MQ=255)
              aCTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGcgcg                                                                                                                                   >  1:788779/1‑149 (MQ=255)
                     cAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGc                                                                                                                            >  2:310354/1‑149 (MQ=255)
                     cAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGc                                                                                                                            >  1:310904/1‑149 (MQ=255)
                         gcCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCTCCGGGTTGTAACCCCAGTAGTTACTTAGCCCCTGGCGTTGCAGACGTGGTGCACTGGCAAACTGA‑GCCACTGGCAGCAGGTCCAGCGCGGTACTGCCCaa                                                                                                                        >  2:515362/1‑149 (MQ=255)
                          ccGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGCGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAAt                                                                                                                       >  2:565733/1‑149 (MQ=255)
                           cGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAAtt                                                                                                                      <  2:751652/149‑1 (MQ=255)
                            gTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGTTTGTAACCACAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTTTCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTg                                                                                                                     <  1:166212/149‑1 (MQ=255)
                             tttCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGt                                                                                                                    >  1:1135214/1‑149 (MQ=255)
                             tttCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGt                                                                                                                    <  1:369425/149‑1 (MQ=255)
                              ttCTGCCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCTCCGGGTTGTAACCCCAGTAGTTCCTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGtt                                                                                                                   <  2:1093027/149‑1 (MQ=255)
                               tCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGttt                                                                                                                  <  1:725426/149‑1 (MQ=255)
                                cTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTc                                                                                                                 <  1:111632/149‑1 (MQ=255)
                                   gCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCGAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaaa                                                                                                              <  1:167272/149‑1 (MQ=255)
                                   gCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaaa                                                                                                              <  1:468460/149‑1 (MQ=255)
                                   gCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaaa                                                                                                              <  2:458064/149‑1 (MQ=255)
                                        cAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGtt                                                                                                         <  2:503581/149‑1 (MQ=255)
                                          ggCATTCGCCGGATGCAGCGCTATCATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGa                                                                                                       <  1:863425/149‑1 (MQ=255)
                                          ggCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGa                                                                                                       <  1:268016/149‑1 (MQ=255)
                                           gCATACGCCGGATGCAGCGCAACCATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGat                                                                                                      >  2:553772/1‑149 (MQ=255)
                                                 gCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCAc                                                                                                >  1:571441/1‑149 (MQ=255)
                                                   cGGATGCAGCGCACACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTATTGGCCAATTGTTTCAAATAGTTGATCATCACCg                                                                                              >  1:640180/1‑149 (MQ=255)
                                                    ggATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGAGCGGTAATGACCAATTGTTTCAAATAGTTGATCATCACCgg                                                                                             >  2:1122277/1‑149 (MQ=255)
                                                     gATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGa                                                                                            <  1:352793/149‑1 (MQ=255)
                                                          aGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGGTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGccc                                                                                       <  2:310904/149‑1 (MQ=255)
                                                          aGCGCAAACATCGCCACCGGGGTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGccc                                                                                       <  1:264431/149‑1 (MQ=255)
                                                             gcAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAg                                                                                    >  2:29628/1‑149 (MQ=255)
                                                              cAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAgg                                                                                   <  1:860225/149‑1 (MQ=255)
                                                              cAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAgg                                                                                   >  1:619549/1‑149 (MQ=255)
                                                                       ccACCGGGTTGTAACCCCAGTAGTTAATTAGCCCCATGCTTTGCAGACCGGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCTCGGTAATGCCCAATTGTTTCAAATAGTTGATAATCACCGCATGCCCGAGGGCTTTATaa                                                                          >  2:895677/1‑149 (MQ=255)
                                                                            gggTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGcc                                                                     <  1:615374/149‑1 (MQ=255)
                                                                             ggTTGTACCCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGccc                                                                    >  1:465184/1‑148 (MQ=255)
                                                                              gTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCAc                                                                   <  2:3442/149‑1 (MQ=255)
                                                                                tGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCTCTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACgg                                                                 >  2:196875/1‑149 (MQ=255)
                                                                                     cgcctgtagtTACTTTGCCGCTTGCGTTGGTGAAGTGGTTCAGTGGGAAACTGC‑GCCACTGGCAGCAGGTCCAGCGCGGTAATGCCCAATTGTGTCAAATCGTGGATCATCAGCGGATGCCCGAGGGCTTTATAAGTGCCACGGAtctc                                                            <  2:784169/144‑1 (MQ=255)
                                                                                      cccAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCg                                                           >  1:466927/1‑149 (MQ=255)
                                                                                        cAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGAc                                                         <  2:303409/149‑1 (MQ=255)
                                                                                        cAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGAc                                                         >  1:287203/1‑149 (MQ=255)
                                                                                           tagtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTGCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGTCCGAGGGCTTTATAAGTGCCACGGATCTCGCCCgg                                                      >  1:788021/1‑149 (MQ=255)
                                                                                            agtTTCTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAAGTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCAGCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCggg                                                     <  1:528115/149‑1 (MQ=255)
                                                                                             gtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCCATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGa                                                    >  2:851587/1‑149 (MQ=255)
                                                                                             gtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATAATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGa                                                    >  2:305028/1‑149 (MQ=255)
                                                                                              ttACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCTAACTTC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGAt                                                   <  2:958167/149‑1 (MQ=255)
                                                                                                aCTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTGCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGAGCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGCCCGGGAtct                                                 >  2:1068881/1‑149 (MQ=255)
                                                                                                   tAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCg                                              <  1:28626/149‑1 (MQ=255)
                                                                                                   tAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCg                                              <  2:208435/149‑1 (MQ=255)
                                                                                                     gcaccATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCTGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCggg                                            <  2:361570/146‑1 (MQ=255)
                                                                                                          aggCGTGGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGGTCCAGCGCGGTAATGCCCAATTGATTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCaa                                       <  2:711077/147‑1 (MQ=255)
                                                                                                          aTGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCaa                                       <  1:739772/149‑1 (MQ=255)
                                                                                                              gTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTAc                                   <  1:279540/149‑1 (MQ=255)
                                                                                                              gTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTAc                                   >  2:737213/1‑149 (MQ=255)
                                                                                                               ttGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACg                                  <  2:329904/149‑1 (MQ=255)
                                                                                                                tGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGt                                 <  1:718452/149‑1 (MQ=255)
                                                                                                                 tcAGACGTGGATCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAACTAGTTGATCATCACCGGCTGCCCGAGGGCTTTATAAGTGCCACGGAGCTCGACCGGTATCTCCGGGTGCAAGTACGtt                                <  1:396271/148‑1 (MQ=255)
                                                                                                                   aGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCCCCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTaa                              >  1:417064/1‑149 (MQ=255)
                                                                                                                     aCGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATc                            <  2:1154204/149‑1 (MQ=255)
                                                                                                                     aCGTGGTTCACGGGCAAACCGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATc                            <  1:1082029/149‑1 (MQ=255)
                                                                                                                       gTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAGTGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATTAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTTCGTTAATCCt                          >  2:873638/1‑149 (MQ=255)
                                                                                                                        tGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCtt                         >  2:463837/1‑149 (MQ=255)
                                                                                                                           ttCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGa                      >  2:823557/1‑149 (MQ=255)
                                                                                                                                 ggCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCAGGGTGCAAGTACGTTAATCCTTTGACATggg                >  2:465470/1‑149 (MQ=255)
                                                                                                                                 ctggaaCTGCTGCCAGTGGC‑GCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCTTCACCGGATGCCCGTGGGCTTTATAAGTGCCACGGATCTCGCCCGGGATCTCGGGGTGGCAGGAGGTTAATCCTTTGACTTggg                >  1:70423/5‑149 (MQ=255)
                                                                                                                                 ctggaaCTGCTGCCAGTGGC‑GCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATAACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATggg                >  2:600796/5‑149 (MQ=255)
                                                                                                                                  tggaaCTGCTGCCAGTGGC‑GCAGTTCCAGCGCGGTAATGCCCATTTGTTGCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCTTGTACGTTAATCCTTTGCCAGGGGc               >  2:1001870/4‑149 (MQ=255)
                                                                                                                                  gCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGc               <  1:751955/149‑1 (MQ=255)
                                                                                                                                   cAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCATATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGAGCTCCGGGTGCAAGTACGTTAATGCTTTGACATGGGCt              >  1:108200/1‑149 (MQ=255)
                                                                                                                                    aaaCTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCtt             >  1:609423/1‑149 (MQ=255)
                                                                                                                                    aaaCTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACAGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCtt             >  1:353716/1‑149 (MQ=255)
                                                                                                                                      aCTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAATTAGTTGATCATCACCGGATGCCCGCGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCa           >  2:1127206/1‑149 (MQ=255)
                                                                                                                                          c‑gcCACTGGCAGCAGTTCCAGCGCGGTAAGGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATaaa       <  1:1086260/149‑1 (MQ=255)
                                                                                                                                            gcCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAat      >  1:936995/1‑149 (MQ=255)
                                                                                                                                             ccACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGCCCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAatg     >  1:836183/1‑149 (MQ=255)
                                                                                                                                             ccACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAatg     >  1:40898/1‑149 (MQ=255)
                                                                                                                                              cACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAatga    >  1:507120/1‑149 (MQ=255)
                                                                                                                                               aCTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATATTTGATCATCACCTGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAatgat   <  1:17763/149‑1 (MQ=255)
                                                                                                                                               aCTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAatgat   <  2:1135214/149‑1 (MQ=255)
                                                                                                                                                cTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGGTAATCCTTTGACATGGGCTTCATAAatgatg  <  1:310354/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATG  >  NC_000913/3570576‑3570867

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