Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,740,805 Δ1 bp 100% intergenic (+42/+711) 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,8050T.96.3% 228.8 / ‑1.1 54intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (27/25);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-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.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGC  >  NC_002947/4740658‑4740949
                                                                                                                                                    |                                                                                                                                                
gCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtt                                                                                                                                                 >  3:283863/1‑149 (MQ=255)
      gaaagaaaTGCTCTCCATCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTTCACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                           <  7:118111/149‑3 (MQ=255)
      gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                           >  5:152704/1‑147 (MQ=255)
      gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                           >  2:269196/1‑147 (MQ=255)
       aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                                          >  1:258312/1‑146 (MQ=255)
        aagaaaTGCTCTCCAGCCCGTACAAGATCCGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAATCGCGACATCGTTGCCCTGAAAGA‑CCCCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                                         <  3:233680/149‑5 (MQ=255)
         agaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatttg                                                                                                                                        <  6:90274/149‑6 (MQ=255)
           aaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCACTCGGGCAGCCttttttatttgt                                                                                                                                       >  7:260958/1‑142 (MQ=255)
                  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                                                                                                    <  5:267105/149‑1 (MQ=255)
                                            cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGaa                                                                                                    <  4:283863/149‑1 (MQ=255)
                                              cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                  <  2:258312/149‑1 (MQ=38)
                                              cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                                  <  8:265641/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                                                                                                  >  4:42094/1‑137 (MQ=255)
                                                          acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                                  <  3:42094/137‑1 (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)
                                                                    aGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGAGGGCAGCCttc                                                                                                                                                <  5:167335/82‑2 (MQ=255)
                                                                    aGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGAGGGCAGCCttc                                                                                                                                                >  6:167335/1‑81 (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                               <  6:269286/135‑1 (MQ=255)
                                                                                                                                aagGCTGCCCGAGGGCAGCC‑TTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTg                               >  5:269286/1‑135 (MQ=39)
                                                                                                                                  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)
                                                                                                                                                    |                                                                                                                                                
GCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGA‑CCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGC  >  NC_002947/4740658‑4740949

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