Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,570,714:1 +T 6.1% coding (606/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,7141.T 6.1% 320.7 / 17.0 82coding (606/1974 nt)glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base . (39/38);  new base T (3/2);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

GCAGTGCTTTGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAA  >  NC_000913/3570566‑3570862
                                                                                                                                                     |                                                                                                                                                    
gCAGTGCTTTGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTAGTTAGCCCCATGCGTTGCAGACGGGGGTCACTGGCAAACTgc                                                                                                                                                       >  1:172146/1‑149 (MQ=255)
  aGTGCTTTGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTgcgc                                                                                                                                                    <  1:861792/149‑1 (MQ=255)
  aGTGCTTTGATTGCATCGCGAAACTACTACAGCGCCGTCTCTGGCTAGCAGGCAAACGCCGGATCCAGCGCAAACATCGCCACCGGGTTGTAACCACAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTgcgc                                                                                                                                                    <  1:1124125/149‑1 (MQ=255)
          gATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgca                                                                                                                                            >  2:739772/1‑149 (MQ=255)
          gATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgca                                                                                                                                            >  1:152116/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:369425/149‑1 (MQ=255)
                                       tttCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGt                                                                                                               >  1:1135214/1‑149 (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                                                    >  1:287203/1‑149 (MQ=255)
                                                                                                  cAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGAc                                                    <  2:303409/149‑1 (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                                         <  2:208435/149‑1 (MQ=255)
                                                                                                             tAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCg                                         <  1:28626/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)
                                                                                                                                                     |                                                                                                                                                    
GCAGTGCTTTGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAA  >  NC_000913/3570566‑3570862

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: