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.92.5% 200.8 / 9.1 53intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (1/3);  new base . (25/24);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 6.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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.

CGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTG  >  NC_002947/4740657‑4740947
                                                                                                                                                     |                                                                                                                                                
cGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCt                                                                                                                                                  >  1:69553/1‑149 (MQ=255)
   aCTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                               <  5:215451/149‑1 (MQ=255)
     tGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                              >  6:314787/1‑148 (MQ=255)
       gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                           >  4:83793/1‑147 (MQ=255)
       gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                               >  8:367859/1‑145 (MQ=255)
       gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                               <  7:367859/145‑1 (MQ=255)
       gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                               >  4:307203/1‑145 (MQ=255)
       gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                               <  3:307203/145‑1 (MQ=255)
               tGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                                          <  5:394622/142‑4 (MQ=255)
               tGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                                          >  6:394622/1‑139 (MQ=255)
                   ctcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGcc                                                                                                                              >  1:351734/1‑149 (MQ=255)
                     ccAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGTGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCgg                                                                                                                            >  1:292235/1‑149 (MQ=255)
                        gCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                              >  6:106237/1‑129 (MQ=255)
                        gCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                              <  5:106237/129‑1 (MQ=255)
                          ccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCctct                                                                                                                       <  2:375813/149‑1 (MQ=255)
                              aCAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                                   >  1:63652/1‑149 (MQ=255)
                                   aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                           >  6:318793/1‑119 (MQ=255)
                                   aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                                           <  5:318793/121‑3 (MQ=255)
                                      ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGc                                                                                                           >  5:344767/1‑149 (MQ=255)
                                      ggCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGc                                                                                                           >  5:4907/1‑149 (MQ=255)
                                       gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGcc                                                                                                          >  8:348315/1‑149 (MQ=255)
                                             cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGaa                                                                                                  <  5:314787/149‑1 (MQ=255)
                                              gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAACGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGaa                                                                                                  <  1:328402/149‑1 (MQ=255)
                                              gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAag                                                                                                 <  2:292235/149‑1 (MQ=255)
                                               cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAaga                                                                                                <  3:83793/149‑1 (MQ=255)
                                               cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCACTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAaga                                                                                                <  5:13022/149‑1 (MQ=255)
                                                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAagag                                                                                               <  5:58332/149‑1 (MQ=38)
                                                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAagag                                                                                               <  7:108256/149‑1 (MQ=255)
                                                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAagag                                                                                               <  2:69553/149‑1 (MQ=255)
                                                   aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGcc                                                                                            >  4:87251/1‑149 (MQ=255)
                                                          cacCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGcaca                                                                                      <  6:344767/148‑1 (MQ=255)
                                                           acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGaa                                                                                                  <  8:145575/135‑1 (MQ=38)
                                                           acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGaa                                                                                                  >  7:145575/1‑135 (MQ=38)
                                                            ccGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAagag                                                                                               <  6:356868/137‑1 (MQ=38)
                                                            ccGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAagag                                                                                               >  5:356868/1‑137 (MQ=38)
                                                                     aGACGCGCGAAATCGTGGCCGTTAAACCCCCCACTAACGCGGCAGTCCCCCAACCCAACACGGGTGCCCGCGGGCAGCTTTTTTTATTTGGCCTGGGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACaaaa                                                                          <  2:351734/149‑1 (MQ=255)
                                                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                               >  2:225692/1‑114 (MQ=255)
                                                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAAGAGGc                                                                                             >  6:374177/1‑130 (MQ=37)
                                                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAAGAGGc                                                                                             <  5:374177/130‑1 (MQ=37)
                                                                      gAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTTTTGTCCTGTGCCGGCCTCTTCGCGGCCATGCC‑‑CGCGAAGCGGCCGGCCCCGACACCCCAAAAc                                                                         >  3:74225/1‑149 (MQ=255)
                                                                          cgcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCt                                                                       >  5:419855/1‑147 (MQ=255)
                                                                             gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGACAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtt                                                                  >  1:328995/1‑149 (MQ=255)
                                                                                 tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAaga                                                                                                >  1:386189/1‑115 (MQ=37)
                                                                                 tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAaga                                                                                                <  2:386189/115‑1 (MQ=37)
                                                                                        ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCACAGATTTTATccgccc                                                       >  8:161000/1‑149 (MQ=255)
                                                                                         cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGAGGGCAGCCttc                                                                                                                                                >  2:362049/1‑61 (MQ=38)
                                                                                         cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGAGGGCAGCCttc                                                                                                                                                <  1:362049/62‑2 (MQ=38)
                                                                                           gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcg                                                                                                                    <  4:128222/87‑1 (MQ=255)
                                                                                           gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcg                                                                                                                    >  3:128222/1‑87 (MQ=255)
                                                                                           gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAAc                                                                         <  5:195303/128‑1 (MQ=39)
                                                                                           gAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAAc                                                                         >  6:195303/1‑128 (MQ=39)
                                                                                              aGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCAc                                                 <  7:161000/149‑1 (MQ=255)
                                                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGCCCCGCACAGACCCCACAAAACCTCac                                                                    >  5:420508/1‑116 (MQ=255)
                                                                                                             tGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGGCCAGCCACCACAAACCCTCAGTTTTTCTCCGCCCCCTCACGCAGCCCTGTAAggg                                  >  3:300727/1‑149 (MQ=255)
                                                                                                              gCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCCGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGCCACCCCAAAACCTCAGTTTTTATCCCCCCGCTCACGCAGCCCGTTCAggg                                  >  6:157603/1‑148 (MQ=255)
                                                                                                                 gTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                              <  6:419855/149‑1 (MQ=255)
                                                                                                                   aGCCGAACGCAAGAAGGCTGCCCTCGGGCATCCTTTTTTATTTGTCCTGTGCCGGCATCTTCGCGGCCTTGCCCTCGCTATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAtt                              <  1:126550/149‑1 (MQ=255)
                                                                                                                       gAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCGCACGCAGCCCGTTCAGGGAATTGAACgg                        >  5:374007/1‑149 (MQ=255)
                                                                                                                       gAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACATAGCGGGCATGGC‑‑CGCGAATAGGGCGGGCCCGAAAACCCAAAACCACACTTTTTATACCCCCCCGCACCCCCCCCGTTAAAGGTATTTAAAgg                        >  2:340994/1‑149 (MQ=255)
                                                                                                                        aaCGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTGGCGGCCTTGCC‑‑CGCGATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                       <  7:348315/149‑1 (MQ=255)
                                                                                                                            taaaaaagGCTGCCCGAGGGCAGCCTTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCa                    <  2:95878/144‑1 (MQ=255)
                                                                                                                             aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTATACGCCCGCTCACGCAGCCCTTTCAGGGTATTTAACGGCGCATc                  >  6:363325/1‑149 (MQ=255)
                                                                                                                                  agGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                     >  5:169664/1‑126 (MQ=39)
                                                                                                                                  agGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                     <  6:169664/126‑1 (MQ=39)
                                                                                                                                    gcggCCCTCGGGAAGCTTTTTTTATTTGGCCTGGGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCcaccac            <  6:374007/146‑1 (MQ=255)
                                                                                                                                        cccTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAATGGCGCATCCACCACGAAc        >  6:179483/1‑149 (MQ=255)
                                                                                                                                          cTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAATGGCGCATCCACCACGAACtt      <  5:179483/149‑1 (MQ=255)
                                                                                                                                            ctatctgCCTTTTTTTTTTGTCCTGTGCGGGCCTCTTCGCGGCCTTGCC‑‑CGGGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGTTATTGAACGGCTCATCCACCACGAACttg     <  1:263214/142‑1 (MQ=255)
                                                                                                                                              ggCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACttgttg  <  2:63652/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTG  >  NC_002947/4740657‑4740947

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: