Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,740,811 (T)7→6 intergenic (+48/+705) 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% 195.1 / 14.9 55intergenic (+42/+711)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (0/4);  new base . (23/26);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCA  >  NC_002947/4740657‑4740936
                                                                                                                                                    |                                                                                                                                   
cGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGCCATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCt                                                                                                                                     >  1:82068/1‑149 (MQ=255)
 gCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtt                                                                                                                                    >  6:164195/1‑149 (MQ=255)
     tGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGTCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttt                                                                                                                                <  6:33994/149‑1 (MQ=255)
     tGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttt                                                                                                                                <  6:258257/149‑1 (MQ=255)
     tGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                 >  7:362471/1‑148 (MQ=255)
       gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttat                                                                                                                              <  8:21303/149‑3 (MQ=255)
       gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttta                                                                                                                               <  3:353536/148‑2 (MQ=255)
       gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttt                                                                                                                                <  1:26697/147‑1 (MQ=255)
        aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttatt                                                                                                                             >  1:250176/1‑146 (MQ=255)
         aagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttttattt                                                                                                                            <  5:179474/149‑5 (MQ=255)
              aTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGCCttc                                                                                                                                   >  4:169334/1‑136 (MQ=255)
              aTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGCCttc                                                                                                                                   <  3:169334/137‑2 (MQ=255)
                       aGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGcc                                                                                                             >  8:4670/1‑149 (MQ=255)
                          ccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCctct                                                                                                          >  6:255593/1‑149 (MQ=255)
                            gTACAAGATCAGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTc                                                                                                        <  7:73457/149‑1 (MQ=255)
                               cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgc                                                                                                      <  7:4670/148‑1 (MQ=255)
                                  gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGcc                                                                                                  >  5:74195/1‑149 (MQ=255)
                                    tCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                <  5:255593/149‑1 (MQ=255)
                                        ccGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGccc                                                                                            <  5:164195/149‑1 (MQ=255)
                                           tCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgcg                                                                                         <  6:176675/149‑1 (MQ=255)
                                             cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                       <  2:47749/149‑1 (MQ=255)
                                             cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                       <  6:78780/149‑1 (MQ=255)
                                             cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                        <  2:377950/148‑1 (MQ=255)
                                              gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGACCGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAag                                                                                      >  5:284500/1‑149 (MQ=255)
                                              gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAag                                                                                      >  7:278144/1‑149 (MQ=255)
                                              gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCAAAGCCGCCATGGCCGCGaac                                                                                      >  1:365729/1‑148 (MQ=255)
                                                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                    <  2:218665/149‑1 (MQ=38)
                                                gcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                    <  6:74195/149‑1 (MQ=255)
                                                  cAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGc                                                                                  <  5:7080/149‑1 (MQ=255)
                                                    gCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                <  3:136730/149‑1 (MQ=255)
                                                     cTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAAGAGGCCgg                                                                               <  6:284500/149‑1 (MQ=38)
                                                        tACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGc                                                                              <  4:334728/147‑1 (MQ=255)
                                                          cacCGGCGAGCAGAAGTGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGTAGGCTGCCCTCGGGCAGCATTTTTTATTTGTCATGTGCCTGCCTATTCGCGGCCGTGCCCGCGAACGGGCCgggcacg                                                                          >  3:238922/1‑144 (MQ=255)
                                                           acCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGcacata                                                                         >  6:350134/1‑147 (MQ=255)
                                                           acCGGCGAGCACAAGCGCGACAGCGTTGCCCTCAAAGACCGCAAGTACGCTGCAGTAGCCGAACGCAAGATGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGGGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGa                                                                         <  2:250176/149‑1 (MQ=255)
                                                                   gCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACcaca                                                                  >  3:307719/1‑148 (MQ=255)
                                                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGcc                                                                                             <  6:331977/119‑1 (MQ=255)
                                                                     aGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGcc                                                                                             >  5:331977/1‑119 (MQ=255)
                                                                       aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                        >  1:69061/1‑122 (MQ=255)
                                                                       aaGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGa                                                                                        <  2:69061/122‑1 (MQ=255)
                                                                          cgcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCa                                                          <  6:394093/149‑1 (MQ=255)
                                                                           gcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAg                                                         >  1:237287/1‑149 (MQ=255)
                                                                            cgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                       <  4:286047/118‑1 (MQ=255)
                                                                            cgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                       >  3:286047/1‑118 (MQ=255)
                                                                             gACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGtt                                                       >  6:97020/1‑149 (MQ=255)
                                                                                aTCGTTGCCCTCAAAGACCGAAAGTACGCTGCAGTAGCCGAACGCAAGACGGCTGCCCTCGGGCAGCCTTTTTAATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGttttt                                                    <  8:278144/149‑1 (MQ=255)
                                                                                  cGTTGCCCTCAAACCCCGCAAGTAAGCTGCAGTACCCGACCGCAAGAAGGCGGCCCCCGGGCAGCCTTTTTTTTTTGTCCTGTGCCGCCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAAACTGAGgttttgt                                                  <  8:362471/149‑3 (MQ=255)
                                                                                   gtgccccTCAAAGACGCCAAGTAAGCTCCATTAGCCGAACGCAAGACGGCTGCCCTGGGGCCGCGTTTTTTTTTTGCCCTGGGCCGGCCGCTTGGCGGCGTTGCCCGCGATGAGGCCGGCACAGACACCACAAAACCTCAGTTTTtct                                                  <  2:82068/144‑1 (MQ=255)
                                                                                        ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCCCAGACACCCCAAAACCTCAGTTTTTGTccgccc                                            >  1:220644/1‑149 (MQ=255)
                                                                                        ccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgccc                                            >  4:244249/1‑149 (MQ=255)
                                                                                         cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGacac                                                                      <  4:85941/122‑1 (MQ=39)
                                                                                         cTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGacac                                                                      >  3:85941/1‑122 (MQ=39)
                                                                                          tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgcccgc                                          >  6:407454/1‑149 (MQ=255)
                                                                                               gACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGATCAGGCCGGCACAGACACCACAAAACCCCAGATTTTCTCCGCCCCCTAACg                                     >  6:134431/1‑149 (MQ=255)
                                                                                                aCCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGc                                    <  7:57572/149‑1 (MQ=255)
                                                                                                   gCAAGTAAGCTGCAGTCGCCGAACCCAACACGGCTGCCCTCGGGCAGCCTTTTTTTTTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGc                                 <  2:237287/149‑1 (MQ=255)
                                                                                                    cAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGAAGAGGGCGGGCCAGaccc                                                                      >  3:101384/1‑109 (MQ=37)
                                                                                                           gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAg                         <  4:307719/149‑1 (MQ=255)
                                                                                                            cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                        <  5:97020/149‑1 (MQ=255)
                                                                                                            cTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCATTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                        <  7:179594/149‑1 (MQ=255)
                                                                                                                           gCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTa                     <  1:305192/137‑1 (MQ=255)
                                                                                                                           gCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATAGCGGCCATGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCAAAGTTTTTATACGCCCGCTCACGCACCCCATTAAGGGTTTAGAAAggaa           >  2:305192/1‑145 (MQ=255)
                                                                                                                             aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc       >  4:244432/1‑149 (MQ=255)
                                                                                                                             aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc       >  8:65721/1‑149 (MQ=255)
                                                                                                                                acggcggcCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTGGCCCGCTAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccac    <  8:8280/147‑1 (MQ=255)
                                                                                                                                 agGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacc   >  5:331557/1‑149 (MQ=255)
                                                                                                                                  ggCTGCCCGCGGGCTGCTTTTTTTTTTTGTCGTGGGCGGCCCTCTTCGCGGCCTGGCCCGCTAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATccacca  <  2:220644/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                   
CGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTTATTTTTCTGGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCA  >  NC_002947/4740657‑4740936

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: