Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,450,129 C→G 9.7% A68G (GCC→GGC)  nicA → nicotinate dehydrogenase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,450,1290CG 9.7% 186.2 / 10.1 62A68G (GCC→GGC) nicAnicotinate dehydrogenase subunit A
Reads supporting (aligned to +/- strand):  ref base C (26/30);  new base G (3/3);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

ACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATC  >  NC_002947/4449990‑4450275
                                                                                                                                           |                                                                                                                                                   
acCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccggccg                                                                                                                                            <  2:116314/149‑1 (MQ=255)
 ccccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGc                                                                                                                                           <  5:302253/149‑1 (MQ=255)
  cccGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGc                                                                                                                                           >  4:88361/1‑148 (MQ=255)
      ctgctgctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTAGCCACCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAata                                                                                                                                       >  1:119428/1‑148 (MQ=255)
           gctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaagggg                                                                                                                                       <  8:45435/143‑7 (MQ=255)
           gctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaagggg                                                                                                                                       >  7:45435/1‑137 (MQ=255)
           gctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccggcc                                                                                                                                             <  8:62390/137‑1 (MQ=255)
           gctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCgccggcc                                                                                                                                             >  7:62390/1‑137 (MQ=255)
            ctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATcaccac                                                                                                                                <  3:85328/149‑1 (MQ=255)
            ctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATcacca                                                                                                                                 <  3:310465/148‑1 (MQ=255)
            ctgATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATcacc                                                                                                                                  <  3:39584/147‑1 (MQ=255)
                 ccTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTc                                                                                                                            >  2:157956/1‑148 (MQ=255)
                  cTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGATGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCg                                                                                                                           >  8:58204/1‑148 (MQ=255)
                  cTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCg                                                                                                                           >  6:233287/1‑148 (MQ=255)
                   tGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaag                                                                                                                                          <  7:315568/132‑4 (MQ=255)
                   tGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaag                                                                                                                                          >  8:315568/1‑129 (MQ=255)
                     cgcAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGaa                                                                                                                         <  1:139131/147‑1 (MQ=255)
                         aCGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTg                                                                                                                   <  6:40293/149‑1 (MQ=255)
                                   ccTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAgg                                                                                                         <  2:119428/149‑1 (MQ=255)
                                   ccTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAgg                                                                                                         >  3:282710/1‑149 (MQ=255)
                                    cTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGc                                                                                                        <  7:325926/149‑1 (MQ=255)
                                       aaCGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccg                                                                                                      >  7:17389/1‑148 (MQ=255)
                                       aaCGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccg                                                                                                      >  7:22716/1‑148 (MQ=255)
                                          ggCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGccgcacc                                                                                                  >  8:274994/1‑149 (MQ=255)
                                            ccccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaag                                                                                                                                          >  8:320819/1‑104 (MQ=255)
                                            ccccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGGCgccggccaag                                                                                                                                          <  7:320819/107‑4 (MQ=255)
                                              ccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATc                                                                                                                                     <  7:293220/110‑1 (MQ=255)
                                              ccAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATc                                                                                                                                     >  8:293220/1‑110 (MQ=255)
                                               cAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTgg                                                                                                                  <  6:227620/128‑1 (MQ=255)
                                               cAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTgg                                                                                                                  >  5:227620/1‑128 (MQ=255)
                                                   tACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGt                                                                                         <  6:218163/149‑1 (MQ=255)
                                                     cGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGc                                                                                       <  4:282710/149‑1 (MQ=255)
                                                            ggCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGc                                                                                <  3:88361/149‑1 (MQ=255)
                                                                tGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTg                                                                                        >  1:57567/1‑137 (MQ=255)
                                                                tGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTg                                                                                        <  2:57567/137‑1 (MQ=255)
                                                                gtggCGAATGCGGGGCGTGTACGGGGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTc                                                                            <  7:50432/147‑1 (MQ=255)
                                                                 gggCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCa                                                                           >  2:14353/1‑149 (MQ=255)
                                                                   gCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATc                                                                         <  1:14353/149‑1 (MQ=255)
                                                                     gAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATcga                                                                       >  4:191745/1‑149 (MQ=255)
                                                                      aaTGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCAGCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGGTCATcgac                                                                      <  7:274994/149‑1 (MQ=255)
                                                                                      ggTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCCGGCGGCCCAgt                                                      >  4:146489/1‑149 (MQ=255)
                                                                                      ggTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAgt                                                      >  2:239958/1‑149 (MQ=255)
                                                                                        tGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCTACGAGCAGGCGGCCCAgtgt                                                    <  8:22716/149‑1 (MQ=255)
                                                                                        tGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAgtgt                                                    <  8:17389/149‑1 (MQ=255)
                                                                                          atcatcGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGtgtg                                                   <  3:191745/148‑1 (MQ=255)
                                                                                            catcGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCAACACCCTCGAAGGGCTTGGAAGCAAGGCCGCCCCTCACCCCGTGCAAAAGGCGTTAATCGACGAGCAGGCGGCCCAGTGTGGt                                                 >  6:20532/1‑148 (MQ=255)
                                                                                              tcGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAg                                                                >  7:223112/1‑131 (MQ=255)
                                                                                              tcGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAg                                                                <  8:223112/131‑1 (MQ=255)
                                                                                                 aCGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTg                                            >  5:178958/1‑148 (MQ=255)
                                                                                                 aCGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGc                                           <  2:41803/149‑1 (MQ=255)
                                                                                                       tGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAAc                                    <  7:58204/149‑1 (MQ=255)
                                                                                                        ggCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCATTGTGGTTACTGCATGAACg                                   <  5:20532/149‑1 (MQ=255)
                                                                                                        ggCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACg                                   <  3:25501/149‑1 (MQ=255)
                                                                                                          ccgcccgcTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGc                                                                  >  2:172445/1‑117 (MQ=255)
                                                                                                          ccgcccgcTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGc                                                                  <  1:172445/117‑1 (MQ=255)
                                                                                                           cgcccgcTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCa                                >  2:103250/1‑149 (MQ=255)
                                                                                                            gcccgcTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGTCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCAt                               >  3:141082/1‑149 (MQ=255)
                                                                                                            gcccgcTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCAt                               >  1:331396/1‑149 (MQ=255)
                                                                                                                    cTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGa                       <  1:157956/149‑1 (MQ=255)
                                                                                                                     tGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCCCCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATAGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGAc                      >  5:45721/1‑149 (MQ=255)
                                                                                                                            tCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCa                                <  1:229755/132‑1 (MQ=255)
                                                                                                                            tCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCa                                >  2:229755/1‑132 (MQ=255)
                                                                                                                                          gccgccGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGACACGGCATGATCATGACCGCCAAGGCCCTGCTCGATc  <  1:239958/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
ACCCCGCTGCTGCTGATCCTGCGCAACGACCTGTGCCTGAACGGCCCCAAGTACGGCTGCGGCCTGGGCGAATGCGGGGCGTGTACGGTGATCATCGACGGCGTGGCCGCCCGCTCCTGCGTCATCCCCTTGGCCGGCGCCGCCGGCCGCAATATCACCACCCTCGAAGGGCTGGGCAGCAAGGCCGCACCGCACCCCGTGCAACAGGCGTTCATCGACGAGCAGGCGGCCCAGTGTGGTTACTGCATGAACGGCATGATCATGACCGCCAAGGCCCTGCTCGATC  >  NC_002947/4449990‑4450275

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: