Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 983,272 G→A 6.9% D95N (GAC→AAC)  fdx → ferredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947983,2720GA 6.9% 215.4 / 4.6 72D95N (GAC→AAC) fdxferredoxin
Reads supporting (aligned to +/- strand):  ref base G (29/38);  new base A (3/2);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAAAGCAAGCCTGAGGTCGATC  >  NC_002947/983126‑983411
                                                                                                                                                  |                                                                                                                                           
aGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGAc                                                                                                                                           <  6:192254/149‑1 (MQ=255)
aGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGAc                                                                                                                                           <  3:239508/149‑1 (MQ=255)
 ggCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTCTTCGTCAGCGACGAAAAcc                                                                                                                                          <  5:171500/149‑1 (MQ=255)
  gCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTCTTCGTCAGCGACGAAAACCt                                                                                                                                         <  4:226234/149‑1 (MQ=255)
    cTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGa                                                                                                                                       >  6:255819/1‑149 (MQ=255)
       tACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCa                                                                                                                                    <  5:175540/149‑1 (MQ=255)
          cACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCg                                                                                                                                 <  4:98418/149‑1 (MQ=255)
                ccACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATcc                                                                                                                           <  5:255819/149‑1 (MQ=255)
                ccACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATcc                                                                                                                           <  2:164461/149‑1 (MQ=255)
                 cacaGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATccc                                                                                                                          >  2:259072/1‑149 (MQ=255)
                      cATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGt                                                                                                                     >  4:162395/1‑149 (MQ=255)
                        tCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTAt                                                                                                                   >  6:172412/1‑149 (MQ=255)
                          gTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTc                                                                                                                 >  2:104347/1‑149 (MQ=255)
                            gcgcAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGc                                                                                                               >  8:210237/1‑149 (MQ=255)
                             cgcAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCACGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCt                                                                                                              >  3:188919/1‑149 (MQ=255)
                                  gggCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAAcc                                                                                                         <  5:172412/149‑1 (MQ=255)
                                      ttCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCAtgc                                                                                                     >  4:155592/1‑149 (MQ=255)
                                      ttCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCAtgc                                                                                                     >  1:213701/1‑149 (MQ=255)
                                        cGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCAtgc                                                                                                     <  2:33101/147‑1 (MQ=255)
                                         gATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCAtgctgc                                                                                                  >  2:60483/1‑149 (MQ=255)
                                                tCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCa                                                                                           >  2:25592/1‑149 (MQ=255)
                                                tCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCa                                                                                           >  5:27762/1‑149 (MQ=255)
                                                 cGAAGAGGCGGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCAc                                                                                          <  6:132050/149‑1 (MQ=255)
                                                   aaGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCg                                                                                        <  3:162395/149‑1 (MQ=255)
                                                         cTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGa                                                                                  <  6:91717/149‑1 (MQ=255)
                                                          tGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCAtgct                                                                                                    <  6:58937/130‑1 (MQ=255)
                                                          tGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCAtgct                                                                                                    >  5:58937/1‑130 (MQ=255)
                                                          tGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAg                                                                                 >  2:88865/1‑149 (MQ=255)
                                                            aCGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCGAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAgg                                                                                <  2:13730/148‑1 (MQ=255)
                                                                  tGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCCACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTa                                                                         >  8:155146/1‑149 (MQ=255)
                                                                    gAAGAGGACATGCTGGAGAAAGCGGGGGGCCTGGAAGCGCAATGGCGCGTGGGTTGCCAGGTTTTCGTCGCTGACGATGACGTGACCATCGAGATCCCCTAGTATTCGCTCAACGATGCTGGCGAAGCTCCGCTGTGAGGGCCTGTAc                                                                        <  1:255005/148‑1 (MQ=255)
                                                                    gAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTAc                                                                        <  1:259072/148‑1 (MQ=255)
                                                                     aaGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACAt                                                                      >  3:70406/1‑149 (MQ=255)
                                                                     aaGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACAt                                                                      <  1:88865/149‑1 (MQ=255)
                                                                         ggACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGt                                                                  <  6:222320/149‑1 (MQ=255)
                                                                              tGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGa                                                              >  5:88643/1‑148 (MQ=255)
                                                                              tGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGa                                                              >  7:39734/1‑148 (MQ=255)
                                                                               gCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGt                                                                          <  3:96285/135‑1 (MQ=255)
                                                                               gCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGt                                                                          >  4:96285/1‑135 (MQ=255)
                                                                                  tgaaaacgCCTGGGGCCTGTAAGCCCAATCGCGCGTGGGTTGTACCCCTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATgg                                                         <  1:224716/142‑1 (MQ=255)
                                                                                   gACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATggg                                                        <  1:60483/149‑1 (MQ=255)
                                                                                   gACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATgg                                                         <  6:27762/148‑1 (MQ=255)
                                                                                    aCAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGt                                                       >  3:201910/1‑149 (MQ=255)
                                                                                      aaaGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGtt                                                      >  3:147405/1‑148 (MQ=255)
                                                                                      aaaGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGtt                                                      <  4:147405/148‑1 (MQ=255)
                                                                                       aaGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGa                                                    <  3:155592/149‑1 (MQ=255)
                                                                                           cTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTa                                                <  4:188919/149‑1 (MQ=255)
                                                                                                  cTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACtt                                             >  8:167562/1‑145 (MQ=255)
                                                                                                   tGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAg                                                                                 <  3:80478/108‑1 (MQ=255)
                                                                                                   tGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAg                                                                                 >  4:80478/1‑108 (MQ=255)
                                                                                                    ggAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACtt                                             >  7:191375/1‑143 (MQ=255)
                                                                                                     gAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACtt                                             <  3:63725/142‑1 (MQ=255)
                                                                                                     gAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACtt                                             >  4:63725/1‑142 (MQ=255)
                                                                                                     gAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATc                                       <  7:210237/148‑1 (MQ=255)
                                                                                                          gcAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATc                                 <  2:213701/149‑1 (MQ=255)
                                                                                                             aTCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGa                                                                    >  2:230469/1‑111 (MQ=255)
                                                                                                             aTCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGa                                                                    <  1:230469/111‑1 (MQ=255)
                                                                                                              tCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGc                             >  8:124474/1‑149 (MQ=255)
                                                                                                                gcgcCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCt                            <  8:39734/148‑1 (MQ=255)
                                                                                                                   ccTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCa                        <  6:88643/149‑1 (MQ=255)
                                                                                                                       ggTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGaaa                    <  4:201910/149‑1 (MQ=255)
                                                                                                                         tggtcAGGTCTTCGTCAGCGACGAAAACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAaag                   >  2:230138/5‑148 (MQ=255)
                                                                                                                          ggtcAGGTCTTCGTCAGCGACGACAACCTGACCATCGAGATCCCCCAGTAGTCGCTCACCCATGCTGCCGACGCACCGCACTGAGGGCCTGTACAGGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAaagc                  >  5:164000/4‑148 (MQ=255)
                                                                                                                           gCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAaagca                 <  4:70406/148‑1 (MQ=255)
                                                                                                                            tcAGGTCTTCGTCAGCGACGAAAACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAaagcaag               >  2:150053/2‑149 (MQ=255)
                                                                                                                             cAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATgg                                                         <  7:70055/106‑1 (MQ=255)
                                                                                                                             cAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATgg                                                         >  8:70055/1‑106 (MQ=255)
                                                                                                                             cAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAaagcaagc              <  3:8790/149‑1 (MQ=255)
                                                                                                                                   ttCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAAAGCAAGCCTGAgg        >  7:39200/1‑149 (MQ=255)
                                                                                                                                   ttCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAAAGCAAGCCTGAg         <  1:25592/148‑1 (MQ=255)
                                                                                                                                     cgtcgCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAAAGCAAGCCTGAGGt       <  7:167562/148‑1 (MQ=255)
                                                                                                                                         gCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAAAGCAAGCCTGAGGTCGATc  <  8:39200/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
AGGCCTGTACCACTTGCCACTGCATCGTGCGCAAGGGCTTCGATTCGCTCGAAGAGGCTGACGAGCTGGAAGAGGACATGCTGGACAAAGCCTGGGGCCTGGAAGCGCAATCGCGCCTGGGTTGCCAGGTTTTCGTCGCTGACGAAGACCTGACCATCGAGATCCCCAAGTATTCGCTCAACCATGCTGCCGAAGCACCGCACTGAGGGCCTGTACATGAGTCTGAAATGGGTTGATGTACTTGAGATCGCAATCCAGCTTGCAGAAAGCAAGCCTGAGGTCGATC  >  NC_002947/983126‑983411

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑