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

GATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATG  >  NC_000913/3570576‑3570867
                                                                                                                                                     |                                                                                                                                               
gATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgca                                                                                                                                                 >  1:152116/1‑149 (MQ=255)
gATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgca                                                                                                                                                 >  2:739772/1‑149 (MQ=255)
 aTTGCATCGCGAAACTCATCCCGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGGTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTTGCAAACTGCGCCACTGgcag                                                                                                                                                <  2:466704/149‑1 (MQ=255)
  ttGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCAACCGGGTTGTAACCCCAGTAGTTACTTAGACCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgcagc                                                                                                                                               >  1:572915/1‑149 (MQ=255)
     cATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGt                                                                                                                                            <  1:4634/149‑1 (MQ=255)
      aTCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACAACAGTAGTTTCTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGtt                                                                                                                                           >  2:17763/1‑149 (MQ=255)
      aTCGCGAAACTCATCCAGCGCCGGTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGtt                                                                                                                                           <  1:273686/149‑1 (MQ=255)
       tCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGCGCAGttt                                                                                                                                          <  2:287208/149‑2 (MQ=255)
       tCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTc                                                                                                                                          >  2:1171762/1‑149 (MQ=255)
          cgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGGGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAg                                                                                                                                       >  1:805179/1‑149 (MQ=255)
           gAAACTCATACAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGATGCAGACGTGGTTCACTGGCAATCTGCGCCACTGGCAGCAGTTCCAgc                                                                                                                                      >  2:19512/1‑149 (MQ=255)
             aaCTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGCGCAGtttgccagt                                                                                                                                    <  2:1119194/149‑8 (MQ=255)
              aCTCATCCAGCGCCGTTTCTGGCGAGCAGGTATCCGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGAGGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGcgcg                                                                                                                                   >  1:910050/1‑149 (MQ=255)
              aCTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGcgcg                                                                                                                                   <  2:1061158/149‑1 (MQ=255)
              aCTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGcgcg                                                                                                                                   >  1:788779/1‑149 (MQ=255)
                     cAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGc                                                                                                                            >  2:310354/1‑149 (MQ=255)
                     cAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGc                                                                                                                            >  1:310904/1‑149 (MQ=255)
                         gcCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCTCCGGGTTGTAACCCCAGTAGTTACTTAGCCCCTGGCGTTGCAGACGTGGTGCACTGGCAAACTGAGCCACTGGCAGCAGGTCCAGCGCGGTACTGCCCaa                                                                                                                        >  2:515362/1‑149 (MQ=255)
                          ccGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGCGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAAt                                                                                                                       >  2:565733/1‑149 (MQ=255)
                           cGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAAtt                                                                                                                      <  2:751652/149‑1 (MQ=255)
                            gTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGTTTGTAACCACAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTTTCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTg                                                                                                                     <  1:166212/149‑1 (MQ=255)
                             tttCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGt                                                                                                                    >  1:1135214/1‑149 (MQ=255)
                             tttCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGt                                                                                                                    <  1:369425/149‑1 (MQ=255)
                              ttCTGCCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCTCCGGGTTGTAACCCCAGTAGTTCCTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGtt                                                                                                                   <  2:1093027/149‑1 (MQ=255)
                               tCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGttt                                                                                                                  <  1:725426/149‑1 (MQ=255)
                                cTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTc                                                                                                                 <  1:111632/149‑1 (MQ=255)
                                   gCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCGAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaaa                                                                                                              <  1:167272/149‑1 (MQ=255)
                                   gCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaaa                                                                                                              <  1:468460/149‑1 (MQ=255)
                                   gCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaaa                                                                                                              <  2:458064/149‑1 (MQ=255)
                                        cAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGtt                                                                                                         <  2:503581/149‑1 (MQ=255)
                                          ggCATTCGCCGGATGCAGCGCTATCATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGa                                                                                                       <  1:863425/149‑1 (MQ=255)
                                          ggCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGa                                                                                                       <  1:268016/149‑1 (MQ=255)
                                           gCATACGCCGGATGCAGCGCAACCATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGat                                                                                                      >  2:553772/1‑149 (MQ=255)
                                                 gCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCAc                                                                                                >  1:571441/1‑149 (MQ=255)
                                                   cGGATGCAGCGCACACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTATTGGCCAATTGTTTCAAATAGTTGATCATCACCg                                                                                              >  1:640180/1‑149 (MQ=255)
                                                    ggATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGAGCGGTAATGACCAATTGTTTCAAATAGTTGATCATCACCgg                                                                                             >  2:1122277/1‑149 (MQ=255)
                                                     gATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGa                                                                                            <  1:352793/149‑1 (MQ=255)
                                                          aGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGGTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGccc                                                                                       <  2:310904/149‑1 (MQ=255)
                                                          aGCGCAAACATCGCCACCGGGGTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGccc                                                                                       <  1:264431/149‑1 (MQ=255)
                                                             gcAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAg                                                                                    >  2:29628/1‑149 (MQ=255)
                                                              cAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAgg                                                                                   <  1:860225/149‑1 (MQ=255)
                                                              cAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAgg                                                                                   >  1:619549/1‑149 (MQ=255)
                                                                       ccACCGGGTTGTAACCCCAGTAGTTAATTAGCCCCATGCTTTGCAGACCGGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCTCGGTAATGCCCAATTGTTTCAAATAGTTGATAATCACCGCATGCCCGAGGGCTTTATaa                                                                          >  2:895677/1‑149 (MQ=255)
                                                                            gggTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGcc                                                                     <  1:615374/149‑1 (MQ=255)
                                                                             ggTTGTACCCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGccc                                                                    >  1:465184/1‑148 (MQ=255)
                                                                              gTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCAc                                                                   <  2:3442/149‑1 (MQ=255)
                                                                                tGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCTCTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACgg                                                                 >  2:196875/1‑149 (MQ=255)
                                                                                     cgcctgtagtTACTTTGCCGCTTGCGTTGGTGAAGTGGTTCAGTGGGAAACTGCGCCACTGGCAGCAGGTCCAGCGCGGTAATGCCCAATTGTGTCAAATCGTGGATCATCAGCGGATGCCCGAGGGCTTTATAAGTGCCACGGAtctc                                                            <  2:784169/144‑1 (MQ=255)
                                                                                      cccAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCg                                                           >  1:466927/1‑149 (MQ=255)
                                                                                        cAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGAc                                                         <  2:303409/149‑1 (MQ=255)
                                                                                        cAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGAc                                                         >  1:287203/1‑149 (MQ=255)
                                                                                           tagtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTGCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGTCCGAGGGCTTTATAAGTGCCACGGATCTCGCCCgg                                                      >  1:788021/1‑149 (MQ=255)
                                                                                            agtTTCTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAAGTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCAGCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCggg                                                     <  1:528115/149‑1 (MQ=255)
                                                                                             gtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCCATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGa                                                    >  2:851587/1‑149 (MQ=255)
                                                                                             gtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATAATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGa                                                    >  2:305028/1‑149 (MQ=255)
                                                                                              ttACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCTAACTTCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGAt                                                   <  2:958167/149‑1 (MQ=255)
                                                                                                aCTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTGCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGAGCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGCCCGGGAtct                                                 >  2:1068881/1‑149 (MQ=255)
                                                                                                   tAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCg                                              <  1:28626/149‑1 (MQ=255)
                                                                                                   tAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCg                                              <  2:208435/149‑1 (MQ=255)
                                                                                                     gcaccATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCTGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCggg                                            <  2:361570/146‑1 (MQ=255)
                                                                                                          aggCGTGGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGGTCCAGCGCGGTAATGCCCAATTGATTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCaa                                       <  2:711077/147‑1 (MQ=255)
                                                                                                          aTGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCaa                                       <  1:739772/149‑1 (MQ=255)
                                                                                                              gTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTAc                                   <  1:279540/149‑1 (MQ=255)
                                                                                                              gTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTAc                                   >  2:737213/1‑149 (MQ=255)
                                                                                                               ttGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACg                                  <  2:329904/149‑1 (MQ=255)
                                                                                                                tGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGt                                 <  1:718452/149‑1 (MQ=255)
                                                                                                                 tcAGACGTGGATCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAACTAGTTGATCATCACCGGCTGCCCGAGGGCTTTATAAGTGCCACGGAGCTCGACCGGTATCTCCGGGTGCAAGTACGtt                                <  1:396271/148‑1 (MQ=255)
                                                                                                                   aGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCCCCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTaa                              >  1:417064/1‑149 (MQ=255)
                                                                                                                     aCGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATc                            <  2:1154204/149‑1 (MQ=255)
                                                                                                                     aCGTGGTTCACGGGCAAACCGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATc                            <  1:1082029/149‑1 (MQ=255)
                                                                                                                       gTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAGTGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATTAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTTCGTTAATCCt                          >  2:873638/1‑149 (MQ=255)
                                                                                                                        tGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCtt                         >  2:463837/1‑149 (MQ=255)
                                                                                                                           ttCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGa                      >  2:823557/1‑149 (MQ=255)
                                                                                                                                 ggCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCAGGGTGCAAGTACGTTAATCCTTTGACATggg                >  2:465470/1‑149 (MQ=255)
                                                                                                                                 ctggaaCTGCTGCCAGTGGCGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCTTCACCGGATGCCCGTGGGCTTTATAAGTGCCACGGATCTCGCCCGGGATCTCGGGGTGGCAGGAGGTTAATCCTTTGACTTggg                >  1:70423/5‑149 (MQ=255)
                                                                                                                                 ctggaaCTGCTGCCAGTGGCGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATAACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATggg                >  2:600796/5‑149 (MQ=255)
                                                                                                                                  tggaaCTGCTGCCAGTGGCGCAGTTCCAGCGCGGTAATGCCCATTTGTTGCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCTTGTACGTTAATCCTTTGCCAGGGGc               >  2:1001870/4‑149 (MQ=255)
                                                                                                                                  gCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGc               <  1:751955/149‑1 (MQ=255)
                                                                                                                                   cAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCATATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGAGCTCCGGGTGCAAGTACGTTAATGCTTTGACATGGGCt              >  1:108200/1‑149 (MQ=255)
                                                                                                                                    aaaCTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCtt             >  1:609423/1‑149 (MQ=255)
                                                                                                                                    aaaCTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACAGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCtt             >  1:353716/1‑149 (MQ=255)
                                                                                                                                      aCTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAATTAGTTGATCATCACCGGATGCCCGCGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCa           >  2:1127206/1‑149 (MQ=255)
                                                                                                                                          cgcCACTGGCAGCAGTTCCAGCGCGGTAAGGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATaaa       <  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)
                                                                                                                                                     |                                                                                                                                               
GATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATG  >  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: