Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,740,841 T→A 13.0% intergenic (+78/+675) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,740,8410TA13.0% 130.3 / 7.5 63intergenic (+78/+675)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (21/30);  new base A (4/5);  total (28/35)
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.29e-01

ATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGT  >  NC_002947/4740692‑4740984
                                                                                                                                                      |                                                                                                                                                 
aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                                                                   >  3:192894/1‑149 (MQ=255)
   ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTg                                                                                                                                                <  8:108859/148‑1 (MQ=255)
        tCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCATGTTTTATTTCTCCTGTGCCGGCCTCGTCGCGGACGTGCCcgc                                                                                                                                          >  5:264976/1‑148 (MQ=255)
          cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                       <  5:267105/149‑1 (MQ=255)
          cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                                                       <  4:283863/149‑1 (MQ=255)
            cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                     <  4:124251/149‑1 (MQ=38)
            cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                     <  2:258312/149‑1 (MQ=38)
            cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                                                     <  8:265641/149‑1 (MQ=38)
             gcCAGCTGTACACCGGCGAGCACAAGCGCGACATCGTTGCCCTTAAAGACCCATCAGTAAGCTTATGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                     <  1:269196/149‑1 (MQ=255)
               cAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGc                                                                                                                                  <  7:26117/149‑1 (MQ=255)
               cAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATAGCGGCCGTGGCCGCGAAGAGGc                                                                                                                                  >  7:113028/1‑149 (MQ=255)
                   tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCgg                                                                                                                               >  8:11802/1‑148 (MQ=255)
                        acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                     >  4:42094/1‑137 (MQ=255)
                        acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                                                     <  3:42094/137‑1 (MQ=255)
                         ccGGCTAGCAGAACCGCGACATGGTTCCCCTGAAAGACCGCAAGTAAGCGCCAGAACCCGAACTCAAGAAGGCGGCCCTGGGGCAGCCTTTTTTATTTGTCCTGTGCCGCCCTCTGCGCGGCCTGGCCCGCTAAGAGGCCGGCACAGac                                                                                                                        <  6:281923/149‑1 (MQ=255)
                                gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCATGCCCCCGAACAGGCCGGGCCCGCCACcaca                                                                                                                  >  1:315837/1‑148 (MQ=255)
                                 cAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACaa                                                                                                                 >  8:141892/1‑148 (MQ=255)
                                        gcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGACTTGGCCTCGAATAGGCCGGAAACGACACCACAAAACCGCAg                                                                                                         >  3:199077/1‑149 (MQ=255)
                                         cgACATCGTTCCCCTGAAACACCCCACGTACGCTCCAGTACCCGAACGCAACACGGCTGCCCGCGGGCAGCCTTTTTTATTTGTCCTGTCCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGt                                                                                                        <  4:56299/149‑1 (MQ=255)
                                                     ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCCCAGACACCACAAAACCACAGGATTTCTccgccc                                                                                            >  6:62349/1‑149 (MQ=255)
                                                                cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAg                                                                                  >  6:159989/1‑149 (MQ=255)
                                                                  cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCAAAAAACCTCAGTTTTTGTCCGCCCCAGCCCGCAGcc                                                                                >  4:167973/1‑149 (MQ=255)
                                                                         gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                                                         >  7:220821/1‑149 (MQ=255)
                                                                                 acccGACCGCAAGAAGTCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                                                                  <  8:242220/146‑1 (MQ=255)
                                                                                 acccGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                                                                  <  7:141892/146‑1 (MQ=255)
                                                                                    cGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg                                                              >  1:234980/1‑149 (MQ=255)
                                                                                    cGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg                                                              <  1:111517/149‑1 (MQ=255)
                                                                                       aCGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCGGGCCTCTTCGCGGCCTTGCCCGCTATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                                                            <  8:2262/148‑1 (MQ=255)
                                                                                        ataaaaaagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgcg                                                           >  3:86712/7‑148 (MQ=255)
                                                                                           aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATCGCCGCCATGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCACAGGTTTTATACGCCCGCGGACGCAGCCCGTTCAGGGTATTGAACGGCGCAta                                                       >  2:265485/1‑148 (MQ=255)
                                                                                            agaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc                                                       >  5:119538/1‑148 (MQ=255)
                                                                                            agaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCAATTAGCGGCCGTGCCCGCGAAGAGGGCGGCACAGCCAACACAAAACCGCAGGTTTTAAGCGCCCCCGGACGCAGCCCGTTCAGGGTATTTGACCgcgccacc                                                      >  8:131703/1‑145 (MQ=255)
                                                                                              aagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                                                                  >  5:269286/1‑135 (MQ=39)
                                                                                              aagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                                                                  <  6:269286/135‑1 (MQ=255)
                                                                                                ggCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca                                                  <  8:220821/149‑1 (MQ=255)
                                                                                                ggCTGCCCGAGGGCAGCCTTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca                                                  >  3:257062/1‑149 (MQ=255)
                                                                                                  cTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCcaccac                                                 >  8:300504/1‑148 (MQ=255)
                                                                                                  cTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACg                                                <  2:234980/149‑1 (MQ=255)
                                                                                                     cccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGaa                                              <  4:192894/148‑1 (MQ=255)
                                                                                                     cccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc                                             <  6:152704/149‑1 (MQ=255)
                                                                                                     cccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTGCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCACGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc                                             <  7:11802/149‑1 (MQ=255)
                                                                                                       cTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt                                           >  6:157787/1‑149 (MQ=255)
                                                                                                       cTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTCTCCGCCCGCTCACGCAGCACTTTCAGGGTATTTAACGGCGCATCCACACCCAACtt                                           >  8:286642/1‑149 (MQ=255)
                                                                                                         ggggCAGCCTTTTTTATTTGACCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgt                                         <  6:111406/148‑1 (MQ=255)
                                                                                                          gggCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt                                        <  7:300504/149‑1 (MQ=255)
                                                                                                           ggCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg                                       <  7:272641/149‑1 (MQ=255)
                                                                                                           ggCAGCCTTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg                                       <  4:257062/149‑1 (MQ=255)
                                                                                                            gCAGCCTTTTTTATTTGTCCTGTGCCGGCCTAATAGCGGCCATTCCCGCGAAGAGGCCGGCACAGACACCACAAAACCAAAGTTTTTATACGCCCGCGCACGCAGCCCGTTAAGGGTATTTAAAGGCGCAACCACCACGAACatgtttg                                      >  4:148403/1‑147 (MQ=255)
                                                                                                             ctgCCTTTTTTTTTTGGCCTGGCCCGGCCCCTGCGGGGCCTGGCCCCCTAAGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGc                                     <  1:278969/147‑1 (MQ=255)
                                                                                                               agccttttttATTTGTCCTGTGCCGGCCTCTTCGCCGCCTTGCCCGCGAAGAGGCCCCCACAGACACCACAAAACCTCACTTTTTCTACGCCCGCTCACGCACCCCTTTCAGGCTAATGAACGGCGCATACACCACCAACTTGTTGGcc                                    >  2:226719/5‑149 (MQ=255)
                                                                                                                   ttttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccacc                                 >  2:231430/1‑148 (MQ=255)
                                                                                                                     ttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGccacc                                 <  1:231430/146‑1 (MQ=255)
                                                                                                                       ttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGCCACCACAAAACCTCAGTTTTTGTCCGCCCGCACACGCAGCCCGTTAAGGGTATTTAACGGCGCATCCACCACGCACttgttg                                       >  4:179812/1‑138 (MQ=38)
                                                                                                                       ttATTTGTCCTGTGCCGCCCTCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg                                       <  3:179812/138‑1 (MQ=255)
                                                                                                                                ctgTGCCGGCCTAATAGCGGCCATGCCCGCCAAGAGGCCGGCACAGACACCACAAAAACTCAGTTTTTATCCGCCCCATCCCCCAGCCCTTTAAGGGTATTGAAACGCGCATCCACCACGAACATGTTGGCCCCCAAAGCCGGCACGc                    >  8:291297/3‑148 (MQ=255)
                                                                                                                                     ccGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCgg                         >  5:104116/1‑138 (MQ=255)
                                                                                                                                     ccGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCgg                         <  6:104116/138‑1 (MQ=255)
                                                                                                                                     ccGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTacca               <  6:119538/148‑1 (MQ=255)
                                                                                                                                         ccTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTgg          >  7:220601/1‑149 (MQ=255)
                                                                                                                                               cgcgGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt    <  4:86712/149‑1 (MQ=255)
                                                                                                                                                gcgGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCATTGTATTGAACGGCGCATCCACCACGAACTTGTTGACCACCATCGCCGGCACGCTACCACCTGGCTCGgt    <  1:265485/148‑1 (MQ=255)
                                                                                                                                                cgcgCCTTGCCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCgg     <  3:167973/145‑1 (MQ=255)
                                                                                                                                                 gggCGTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtgt  <  5:62349/148‑1 (MQ=255)
                                                                                                                                                 gcgCCTTGGCCCGGTAATAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt    <  2:315837/146‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
ATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTGT  >  NC_002947/4740692‑4740984

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: