Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,566,218 C→G 5.8% G122R (GGT→CGT)  glgP ← glycogen phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,566,2180CG 5.8% 241.8 / 4.5 69G122R (GGT→CGT) glgPglycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (30/35);  new base G (2/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGGTAG  >  NC_000913/3566074‑3566363
                                                                                                                                                |                                                                                                                                                 
cTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAgg                                                                                                                                               <  2:425918/149‑1 (MQ=255)
    cGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGTCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTg                                                                                                                                           <  4:82977/149‑1 (MQ=255)
     gACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTgg                                                                                                                                          <  1:207523/149‑1 (MQ=255)
      aCTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTggg                                                                                                                                         >  1:391548/1‑149 (MQ=255)
          tttCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTggga                                                                                                                                        >  3:330394/1‑145 (MQ=255)
          tttCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCTCGGTAACGGTGGCCTggga                                                                                                                                        <  4:330394/146‑2 (MQ=255)
           ttCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCAtt                                                                                                                                    >  2:507777/1‑149 (MQ=255)
             cTGGCTACCGTTAACGATGTTCTGCTTTAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCAtttt                                                                                                                                  <  3:402680/149‑1 (MQ=255)
              tGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttc                                                                                                                                 <  4:142092/149‑1 (MQ=255)
              tGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttc                                                                                                                                 >  2:483109/1‑149 (MQ=255)
                 cTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttcttc                                                                                                                              <  4:118499/149‑1 (MQ=255)
                   aCCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtcat                                                                                                                            <  3:476299/149‑1 (MQ=255)
                   aCCGTTAACGATGTTCTGCTTGAACATACCGTAGGCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtcat                                                                                                                            <  2:10116/149‑1 (MQ=255)
                    ccGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTtcatc                                                                                                                           >  1:233090/1‑149 (MQ=255)
                             aTGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGctc                                                                                                                  <  3:393635/149‑1 (MQ=255)
                                    gCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCgaga                                                                                                           <  2:391548/149‑1 (MQ=255)
                                    gCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCgaga                                                                                                           <  1:4000/149‑1 (MQ=255)
                                      ttGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGAtt                                                                                                         <  3:514053/149‑1 (MQ=255)
                                      ttGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGAtt                                                                                                         >  4:560652/1‑149 (MQ=255)
                                      ttGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGAtt                                                                                                         <  1:43606/149‑1 (MQ=255)
                                      ttGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGAtt                                                                                                         >  4:144033/1‑149 (MQ=255)
                                                 gTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCg                                                                                              >  4:257508/1‑149 (MQ=255)
                                                     tCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCatcaat                                                                                                                        <  1:547265/119‑1 (MQ=255)
                                                     tCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCatcaat                                                                                                                        >  2:547265/1‑119 (MQ=255)
                                                        tAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAg                                                                                       >  1:406211/1‑149 (MQ=255)
                                                                 ccGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTg                                                                              <  3:144033/149‑1 (MQ=255)
                                                                        cgcgCCCCGGCAACCCTAACGTCGCCCGAGAATCAAGGCAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCt                                                                       >  3:491364/1‑149 (MQ=255)
                                                                        cgcgCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCt                                                                       >  4:123445/1‑149 (MQ=255)
                                                                        cgcgCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCt                                                                       >  3:335522/1‑149 (MQ=255)
                                                                         gcgcCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCtt                                                                      <  2:233090/149‑1 (MQ=255)
                                                                             cccGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTa                                                                  >  4:200647/1‑149 (MQ=255)
                                                                                   aaCCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTcc                                                            <  3:235942/149‑1 (MQ=255)
                                                                                    aCCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCt                                                           <  3:4406/149‑1 (MQ=255)
                                                                                      ccTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCg                                                                                                              <  4:302000/96‑1 (MQ=255)
                                                                                      ccTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCg                                                                                                              >  3:302000/1‑96 (MQ=255)
                                                                                       cTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGc                                                        <  2:524821/149‑1 (MQ=255)
                                                                                        tAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCg                                                       >  4:6691/1‑149 (MQ=255)
                                                                                            gTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGacaa                                                   <  1:325910/149‑1 (MQ=255)
                                                                                             tCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGacaac                                                  >  3:263415/1‑149 (MQ=255)
                                                                                                 cAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATgg                                              >  3:306750/1‑149 (MQ=255)
                                                                                                  agagAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttc                                                                                                                                 <  4:511116/65‑1 (MQ=255)
                                                                                                  agagAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTttc                                                                                                                                 >  3:511116/1‑65 (MQ=255)
                                                                                                  agagAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGc                                             >  2:126033/1‑149 (MQ=255)
                                                                                                     gaATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGtt                                          <  3:560652/149‑1 (MQ=255)
                                                                                                      aaTCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTg                                         <  4:119287/149‑1 (MQ=255)
                                                                                                       aTCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTgg                                        <  1:507777/149‑1 (MQ=255)
                                                                                                        tCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGga                                       <  4:89982/149‑1 (MQ=255)
                                                                                                        tCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGga                                       >  3:208707/1‑149 (MQ=255)
                                                                                                         cAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCga                                                                                                             >  3:253733/1‑78 (MQ=255)
                                                                                                         cAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCga                                                                                                             <  4:253733/78‑1 (MQ=255)
                                                                                                          aaGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGgaga                                     >  1:476434/1‑149 (MQ=255)
                                                                                                          aaGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGgaga                                     >  3:508785/1‑149 (MQ=255)
                                                                                                            ggAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGc                                   <  4:317446/149‑1 (MQ=255)
                                                                                                               agcaagcCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTa                                >  2:266360/1‑149 (MQ=255)
                                                                                                                  aagcCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACgg                             >  3:371583/1‑149 (MQ=255)
                                                                                                                  aagcCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACgg                             >  2:521268/1‑149 (MQ=255)
                                                                                                                   agcCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGc                            >  3:325301/1‑149 (MQ=255)
                                                                                                                     ccgccAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGc                                             >  3:503481/1‑130 (MQ=255)
                                                                                                                     ccgccAGACCTCCCAGGCCACCGTTACCGAGGCCTGGGTTATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGCCAACATGGc                                             <  4:503481/130‑1 (MQ=255)
                                                                                                                      cgccAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCaa                         <  1:3268/149‑1 (MQ=255)
                                                                                                                             cGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCaaaaa                  <  4:306750/149‑1 (MQ=255)
                                                                                                                              gacccAGGCCTCGGTAACGGTGGCCTGGGACATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTg                                         <  1:315849/123‑1 (MQ=38)
                                                                                                                                 ccAGGCCTCGGTAACGGTGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCACCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTcc              >  1:568132/1‑149 (MQ=255)
                                                                                                                                  cAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCa             <  1:521268/149‑1 (MQ=255)
                                                                                                                                       cACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGAc        >  1:505200/1‑149 (MQ=255)
                                                                                                                                       cACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGAc        >  3:4973/1‑149 (MQ=255)
                                                                                                                                        aCCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACa       <  4:371583/149‑1 (MQ=255)
                                                                                                                                         ccGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAg      <  4:39490/149‑1 (MQ=255)
                                                                                                                                             tACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGgtag  <  1:405439/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CTGGCGACTCTTTCTGGCTACCGTTAACGATGTTCTGCTTGAACATACCGTAGTCATAGCGGATGCCGTAACCGCGCCCCGGCAACCCTAACGTCGCCAGAGAATCAAGGAAGCAAGCCGCCAGACGTCCCAGGCCACCGTTACCGAGGCCTGGGTCATTTTCTTCATCAATCAGCTCTTCGAGATTTAACCCCATCGCTTCCAGTGCGCCCTGTACATCTTCGTAAATTCCTAGCGACAACATGGCGTTGGAGAGCGTACGGCCAATCAAAAACTCCATCGACAGGTAG  >  NC_000913/3566074‑3566363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: