Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,740,819 2 bp→CT 100% intergenic (+56/+696) 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,8190TC99.1% 183.2 / ‑6.5 59intergenic (+56/+697)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (31/27);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*NC_0029474,740,8200GT98.4% 190.9 / ‑5.1 59intergenic (+57/+696)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (31/27);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACC  >  NC_002947/4740676‑4740953
                                                                                                                                                ||                                                                                                                                     
ctcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGCGGGCAGCC‑TTTTTTATTTGTCCTGTGcc                                                                                                                                  >  8:297460/1‑149 (MQ=255)
 tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCg                                                                                                                                 >  5:111406/1‑149 (MQ=255)
       ccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCctct                                                                                                                           >  2:278969/1‑149 (MQ=255)
          tACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCttag                                                                                                                        >  2:111517/1‑147 (MQ=255)
             aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTATTCGCgg                                                                                                                     >  7:2262/1‑149 (MQ=255)
                aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                                  >  3:192894/1‑149 (MQ=255)
                   ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTg                                                                                                                <  8:108859/148‑1 (MQ=255)
                        tCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAG‑CATGTTTTATTTCTCCTGTGCCGGCCTCGTCGCGGACGTGCCcgc                                                                                                           >  5:264976/1‑148 (MQ=255)
                          cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                        <  4:283863/149‑1 (MQ=255)
                          cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                        <  5:267105/149‑1 (MQ=255)
                            cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                      <  8:265641/149‑1 (MQ=38)
                            cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                      <  2:258312/149‑1 (MQ=38)
                            cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                      <  4:124251/149‑1 (MQ=38)
                             gcCAGCTGTACACCGGCGAGCACAAGCGCGACATCGTTGCCCTTAAAGACCCATCAGTAAGCTTATGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                      <  1:269196/149‑1 (MQ=255)
                               cAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGc                                                                                                   <  7:26117/149‑1 (MQ=255)
                               cAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCATAGCGGCCGTGGCCGCGAAGAGGc                                                                                                   >  7:113028/1‑149 (MQ=255)
                                   tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCgg                                                                                                >  8:11802/1‑148 (MQ=255)
                                        acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                      <  3:42094/137‑1 (MQ=255)
                                        acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                      >  4:42094/1‑137 (MQ=255)
                                         ccGGCTAGCAGAACCGCGACATGGTTCCCCTGAAAGA‑CCGCAAGTAAGCGCCAGAACCCGAACTCAAGAAGGCGGCCCTGGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGCCCTCTGCGCGGCCTGGCCCGCTAAGAGGCCGGCACAGac                                                                                         <  6:281923/149‑1 (MQ=255)
                                                gCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCACTTAGCGGCCATGCCCCCGAACAGGCCGGGCCCGCCACcaca                                                                                   >  1:315837/1‑148 (MQ=255)
                                                 cAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACaa                                                                                  >  8:141892/1‑148 (MQ=255)
                                                        gcgACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGACTTGGCCTCGAATAGGCCGGAAACGACACCACAAAACCGCAg                                                                          >  3:199077/1‑149 (MQ=255)
                                                         cgACATCGTTCCCCTGAAACA‑CCCCACGTACGCTCCAGTACCCGAACGCAACACGGCTGCCCGCGGGCAGCC‑TTTTTTATTTGTCCTGTCCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGt                                                                         <  4:56299/149‑1 (MQ=255)
                                                                     ccTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCCCAGACACCACAAAACCACAGGATTTCTccgccc                                                             >  6:62349/1‑149 (MQ=255)
                                                                                cGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAg                                                   >  6:159989/1‑149 (MQ=255)
                                                                                  cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCAAAAAACCTCAGTTTTTGTCCGCCCCAGCCCGCAGcc                                                 >  4:167973/1‑149 (MQ=255)
                                                                                         gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                                          >  7:220821/1‑149 (MQ=255)
                                                                                                 acccGACCGCAAGAAGTCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                                   <  8:242220/146‑1 (MQ=255)
                                                                                                 acccGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                                   <  7:141892/146‑1 (MQ=255)
                                                                                                    cGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg                               >  1:234980/1‑149 (MQ=255)
                                                                                                    cGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACg                               <  1:111517/149‑1 (MQ=255)
                                                                                                       aCGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCGGGCCTCTTCGCGGCCTTGCCCGCTATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                             <  8:2262/148‑1 (MQ=255)
                                                                                                        ataaaaaagGCTGCCCGAGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgcg                            >  3:86712/7‑148 (MQ=255)
                                                                                                           aagaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCATCGCCGCCATGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCACAGGTTTTATACGCCCGCGGACGCAGCCCGTTCAGGGTATTGAACGGCGCAta                        >  2:265485/1‑148 (MQ=255)
                                                                                                            agaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc                        >  5:119538/1‑148 (MQ=255)
                                                                                                            agaagGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCAATTAGCGGCCGTGCCCGCGAAGAGGGCGGCACAGCCAACACAAAACCGCAGGTTTTAAGCGCCCCCGGACGCAGCCCGTTCAGGGTATTTGACCgcgccacc                       >  8:131703/1‑145 (MQ=255)
                                                                                                              aagGCTGCCCGAGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                                   >  5:269286/1‑135 (MQ=39)
                                                                                                              aagGCTGCCCGAGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                                   <  6:269286/135‑1 (MQ=255)
                                                                                                                ggCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca                   <  8:220821/149‑1 (MQ=255)
                                                                                                                ggCTGCCCGAGGGCAGCC‑TTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca                   >  3:257062/1‑149 (MQ=255)
                                                                                                                  cTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCcaccac                  >  8:300504/1‑148 (MQ=255)
                                                                                                                  cTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACg                 <  2:234980/149‑1 (MQ=255)
                                                                                                                     cccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGaa               <  4:192894/148‑1 (MQ=255)
                                                                                                                     cccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc              <  6:152704/149‑1 (MQ=255)
                                                                                                                     cccTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTGCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCACGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAAc              <  7:11802/149‑1 (MQ=255)
                                                                                                                       cTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACtt            >  6:157787/1‑149 (MQ=255)
                                                                                                                       cTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTCTCCGCCCGCTCACGCAGCACTTTCAGGGTATTTAACGGCGCATCCACACCCAACtt            >  8:286642/1‑149 (MQ=255)
                                                                                                                         ggggCAGCC‑TTTTTTATTTGACCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgt          <  6:111406/148‑1 (MQ=255)
                                                                                                                          gggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgtt         <  7:300504/149‑1 (MQ=255)
                                                                                                                           ggCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg        <  7:272641/149‑1 (MQ=255)
                                                                                                                           ggCAGCC‑TTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg        <  4:257062/149‑1 (MQ=255)
                                                                                                                            gCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTAATAGCGGCCATTCCCGCGAAGAGGCCGGCACAGACACCACAAAACCAAAGTTTTTATACGCCCGCGCACGCAGCCCGTTAAGGGTATTTAAAGGCGCAACCACCACGAACatgtttg       >  4:148403/1‑147 (MQ=255)
                                                                                                                             ctgCCTTTTTTTTTTGGCCTGG‑CCCGGCCCCTGCGGGGCCTGGCCCCCTAAGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGc      <  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)
                                                                                                                                                ||                                                                                                                                     
CTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACC  >  NC_002947/4740676‑4740953

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: ‑