Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,740,839 C→G 21.2% intergenic (+76/+677) 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,8390CG21.2% 128.5 / 31.9 66intergenic (+76/+677)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (25/27);  new base G (7/7);  total (32/34)
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.56e-01

GGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTG  >  NC_002947/4740695‑4740983
                                                                                                                                                 |                                                                                                                                                  
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                                                                                                                                 <  5:374177/130‑1 (MQ=37)
                               aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAAGAGGc                                                                                                                                 >  6:374177/1‑130 (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                                                                                                                                    <  2:386189/115‑1 (MQ=37)
                                           tCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGC‑‑CGCGAaga                                                                                                                                    >  1:386189/1‑115 (MQ=37)
                                                  ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTTCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCACAGATTTTATccgccc                                                                                           >  8:161000/1‑149 (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)
                                                                                                                  ttttATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATc                             >  1:308012/1‑149 (MQ=255)
                                                                                                                              tgTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTa                 >  5:300724/2‑149 (MQ=255)
                                                                                                                                tGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                <  6:406545/148‑1 (MQ=255)
                                                                                                                                  ccGGCCTCTACGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACTGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTaccac             <  7:192418/149‑1 (MQ=255)
                                                                                                                                   cgcccTCTTCGCGGCCTGCCC‑‑CGCGATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTaccacc            <  3:87251/146‑1 (MQ=255)
                                                                                                                                        tctTCGCGGGCAAGGC‑‑CGCGAAGAGGCCGGCACAGACCCCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGCACTTGTTGGACACCATCACCGGCAAGCTACCACCTGGCt       >  1:168214/1‑149 (MQ=255)
                                                                                                                                        tctTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCCGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt       >  3:311544/1‑149 (MQ=255)
                                                                                                                                        tctTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCt       >  5:156556/1‑149 (MQ=255)
                                                                                                                                         ctTCGCGGGCAAGGC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                >  4:100715/1‑139 (MQ=255)
                                                                                                                                         ctTCGCGGGCAAGGC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTac                <  3:100715/139‑1 (MQ=255)
                                                                                                                                         ctTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTc      >  1:268226/1‑149 (MQ=255)
                                                                                                                                          ttCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTg          <  4:332430/144‑1 (MQ=255)
                                                                                                                                          ttCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTg          >  3:332430/1‑144 (MQ=255)
                                                                                                                                          ttCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCg     >  3:285681/1‑149 (MQ=255)
                                                                                                                                            cgcgGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCg     <  4:300727/147‑1 (MQ=255)
                                                                                                                                             gcgGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAATACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt   <  1:340994/148‑1 (MQ=255)
                                                                                                                                             gcgGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgtg  >  3:395375/1‑149 (MQ=255)
                                                                                                                                             gcgGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGgt   <  5:363325/148‑1 (MQ=255)
                                                                                                                                              cgGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                                                                   >  4:193144/1‑83 (MQ=255)
                                                                                                                                              cgGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                                                                   <  3:193144/83‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCC‑‑CGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTGTTGGCCACCATCGCCGGCACGCTACCACCTGGCTCGGTG  >  NC_002947/4740695‑4740983

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: