Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 394,897 G→A 6.8% P93P (CCG→CCA fdhA → formaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947394,8970GA 6.8% 173.0 / 2.6 59P93P (CCG→CCAfdhAformaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (30/25);  new base A (2/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CTGCGGTTCTGACCAGCACATGGTCCGTGGCCGCACTACTGCCCAGGTCGGCCTGGTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAGGCCGAATACGTGCTGGTGC  >  NC_002947/394756‑395045
                                                                                                                                              |                                                                                                                                                    
cTGCGGTTCTGACCAGCACATGGTCCGTGGCCGCACTACTGCCCAGGTCGGCCTGGTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACg                                                                                                                                               >  1:197889/1‑149 (MQ=255)
  gCGGTTCTGACCAGCACATGGTCCGTGGCCGCACTACTGCCCAGGTCGGCCTGGTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTc                                                                                                                                             >  8:157366/1‑149 (MQ=255)
   cGGTTCTGACCAGCACATGGTCCGTGGCCGCACTACTGCCCAGGTCGGCCTGGTCCTGGGCCACGAAATCACCTGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTGAACGTc                                                                                                                                             <  2:1386/148‑1 (MQ=255)
           aCCAGCACATGGTCCGTGGCCGCACTACTGCCCAGGTCGGCCTGGTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCg                                                                                                                                      >  1:190512/1‑147 (MQ=255)
           aCCAGCACATGGTCCGTGGCCGCACTACTGCCCAGGTCGGCCTGGTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCg                                                                                                                                      <  2:190512/147‑1 (MQ=255)
                 acaTGGTCCGTGGCCGCACTACTGCCCAGGTCGGCCTGGTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGgccgctgc                                                                                                                              <  1:24557/149‑1 (MQ=255)
                     ggTCCGTGGCCGCACTACTGCCCAGGTCGGCCTGGTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCg                                                                                                                                      >  3:18439/1‑137 (MQ=255)
                     ggTCCGTGGCCGCACTACTGCCCAGGTCGGCCTGGTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCg                                                                                                                                      <  4:18439/137‑1 (MQ=255)
                             gCCGCACTACTGCCCTGGTCGGCCTGGTCCTGGGCCACGAAATCACCGGTGAATTCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCaaa                                                                                                                  <  3:4985/149‑1 (MQ=255)
                             gCCGCACTACTGCCCAGGTCGGCCTGGTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCa                                                                                                                    >  1:244479/1‑147 (MQ=255)
                                    tactGCCCAGGTCGGCCTGGTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATg                                                                                                            >  5:60613/1‑148 (MQ=255)
                                                 ggCCTGGTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCAGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTACAAAGAGATGCACACCGGTGTCt                                                                                               <  8:2882/148‑1 (MQ=255)
                                                       gTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTGGCCTGCGGCCGCTGCCGCTCCTGCAAagag                                                                                                               <  1:24804/126‑1 (MQ=255)
                                                       gTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAagag                                                                                                               >  2:24804/1‑126 (MQ=255)
                                                        tCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCAcc                                                                                       >  4:49495/1‑149 (MQ=255)
                                                         ccTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCg                                                                                      >  7:85287/1‑149 (MQ=255)
                                                         ccTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCg                                                                                      >  6:90319/1‑149 (MQ=255)
                                                                cACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGca                                                                                                          <  5:30200/122‑1 (MQ=255)
                                                                cACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGca                                                                                                          >  6:30200/1‑122 (MQ=255)
                                                                cACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTAAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAAccc                                                                               >  5:218768/1‑149 (MQ=255)
                                                                   gAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGc                                                                            >  4:184554/1‑149 (MQ=255)
                                                                     aaTCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGccc                                                                          <  3:49495/149‑1 (MQ=255)
                                                                         aCCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGGGTCTGCCTCACCGTCAACCCGGCCCgcgc                                                                      >  4:91003/1‑149 (MQ=255)
                                                                           cGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCcg                                                                    <  2:197889/149‑1 (MQ=255)
                                                                            ggTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCcgg                                                                   >  1:19972/1‑149 (MQ=255)
                                                                              ttaaaTCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCTTCAACCCGGCCCGCGCcggcg                                                                 <  4:101195/147‑1 (MQ=255)
                                                                              tGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCcggcg                                                                 >  2:127287/1‑149 (MQ=255)
                                                                               gAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGca                                                                                                          >  4:224088/1‑107 (MQ=255)
                                                                               gAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGca                                                                                                          <  3:224088/107‑1 (MQ=255)
                                                                                aaaTCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGt                                                               <  8:85287/149‑1 (MQ=255)
                                                                                 aaTCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTg                                                              >  8:186249/1‑149 (MQ=255)
                                                                                  aTCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGc                                                             >  4:246613/1‑149 (MQ=255)
                                                                                   tcgtcgAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGcc                                                            <  7:186249/149‑1 (MQ=255)
                                                                                      tcgAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTAc                                                         <  1:127287/149‑1 (MQ=255)
                                                                                      tcgAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTAc                                                         >  6:288597/1‑149 (MQ=255)
                                                                                              ggCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCcggcgg                                                                >  2:255358/1‑134 (MQ=255)
                                                                                              ggCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCcggcgg                                                                <  1:255358/134‑1 (MQ=255)
                                                                                                ccGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGAt                                                                                                             <  4:5703/87‑1 (MQ=255)
                                                                                                ccGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGAt                                                                                                             >  3:5703/1‑87 (MQ=255)
                                                                                                  gcgACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGAc                                             <  7:157366/149‑1 (MQ=255)
                                                                                                        tCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGc                                       <  6:218768/149‑1 (MQ=255)
                                                                                                        tCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGc                                       >  6:141244/1‑149 (MQ=255)
                                                                                                          gAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCATTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGa                                     >  2:62441/1‑149 (MQ=255)
                                                                                                          gAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCATTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGa                                     >  2:58000/1‑149 (MQ=255)
                                                                                                             cgcATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTg                                  <  3:184554/149‑1 (MQ=255)
                                                                                                               cATGCAGATCGGCGACCTGGTATCGGTCCCGGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTg                                  >  1:189150/1‑148 (MQ=255)
                                                                                                               cATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGa                                <  5:90319/149‑1 (MQ=255)
                                                                                                                     gATCGGCGACCTGGTATCGGTCACGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTg                          >  7:101764/1‑149 (MQ=255)
                                                                                                                         ggCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCa                      <  3:246613/149‑1 (MQ=255)
                                                                                                                           cGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAgg                    <  5:288597/149‑1 (MQ=255)
                                                                                                                           cGACCTGGTATCGGTCCCATTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAgg                    <  1:62441/149‑1 (MQ=255)
                                                                                                                           cGACCTGGTATCGGTCCCATTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAgg                    <  1:58000/149‑1 (MQ=255)
                                                                                                                            gACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAGGc                   <  6:60613/149‑1 (MQ=255)
                                                                                                                              ccTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTATCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAGGcc                  >  3:245846/1‑148 (MQ=255)
                                                                                                                                     tCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAGGCCGAATACGt          >  4:98493/1‑149 (MQ=255)
                                                                                                                                       ggTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTCGACCGGTGGCCAGGCCGAATACGTGc        >  4:204374/1‑149 (MQ=255)
                                                                                                                                         tCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAGGCCGAATACGTGCTg      >  5:95497/1‑149 (MQ=255)
                                                                                                                                         tCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAGGCCGAATACGTGCTg      >  5:210402/1‑149 (MQ=255)
                                                                                                                                           ccGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAGGCCGAATACGTGCTgg     <  2:189150/148‑1 (MQ=255)
                                                                                                                                              gTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAGGCCGAATAc            >  7:258675/1‑139 (MQ=255)
                                                                                                                                              gTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAGGCCGAATAc            <  8:258675/139‑1 (MQ=255)
                                                                                                                                              gTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAGGCCGAATACGTGCTGGTGc  >  8:198427/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CTGCGGTTCTGACCAGCACATGGTCCGTGGCCGCACTACTGCCCAGGTCGGCCTGGTCCTGGGCCACGAAATCACCGGTGAAATCGTCGAGAAGGGCCGCGACGTCGAGCGCATGCAGATCGGCGACCTGGTATCGGTCCCGTTCAACGTCGCCTGCGGCCGCTGCCGCTCCTGCAAAGAGATGCACACCGGTGTCTGCCTCACCGTCAACCCGGCCCGCGCCGGCGGTGCCTACGGCTACGTCGACATGGGCGACTGGACCGGTGGCCAGGCCGAATACGTGCTGGTGC  >  NC_002947/394756‑395045

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: