Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,451,822 G→T 6.7% A476S (GCC→TCC)  nicB → nicotinate dehydrogenase subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,451,8220GT 6.7% 240.3 / 5.5 75A476S (GCC→TCC) nicBnicotinate dehydrogenase subunit B
Reads supporting (aligned to +/- strand):  ref base G (35/35);  new base T (2/3);  total (37/38)
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.60e-01

AGGAACACGTGTGGGAGCCCAAGGGCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGCCTCATGGATGCGGGGTGTGTCGGCCATGCCCAAC  >  NC_002947/4451676‑4451968
                                                                                                                                                  |                                                                                                                                                  
aGGAACACGTGTGGGAGCCCAAGGGCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGcc                                                                                                                                                  >  2:58319/1‑149 (MQ=255)
aGGAACACGTGTGGGAGCCCAAGGGCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGcc                                                                                                                                                  >  6:143286/1‑149 (MQ=255)
  gAACACGTGTGGGAGCCCAAGGGCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGGCCAATGGCGCGCCCCCCCTGGCCCTCATGCTGACCGGTGCCGt                                                                                                                                                >  7:267869/1‑149 (MQ=255)
    acacGTGTGGGAGCCCAAGGGCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTgg                                                                                                                                              <  2:199044/149‑1 (MQ=255)
      acGTGTGGGAGCCCAAGGGCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGaa                                                                                                                                            >  4:196699/1‑149 (MQ=255)
           tgGGAGCCCAAGGGCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGt                                                                                                                                       <  8:274919/149‑1 (MQ=255)
                gCCCAAGGGCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCgg                                                                                                                                  >  6:75956/1‑149 (MQ=255)
                   cAAGTGCAACGCGTAACTGATTGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGAGCAGCTACCCGTCCAATAGCGGGCCCACCCTGGCCCTGCTGCTGACCGGTTCCACGGCACCGGtcagca                                                                                                                                  <  8:196657/146‑6 (MQ=255)
                   cAAGGGCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGgcgc                                                                                                                               >  5:242382/1‑149 (MQ=255)
                   cAAGGGCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATAGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTTCCACGGCACCGGtcagca                                                                                                                                  >  7:196657/1‑141 (MQ=255)
                    aaGGGCACCTCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTTCCACGGCACCGGtca                                                                                                                                     <  5:103688/142‑3 (MQ=255)
                    aaGGGCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTTCCACGGCACCGGtca                                                                                                                                     >  6:103688/1‑140 (MQ=255)
                       ggCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCgg                                                                                                                                        >  8:146479/1‑136 (MQ=255)
                       ggCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCgg                                                                                                                                        <  7:146479/136‑1 (MQ=255)
                       ggCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGtt                                                                                                                           <  3:124760/149‑1 (MQ=255)
                        gCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTc                                                                                                                          <  5:143286/149‑1 (MQ=255)
                          aCCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTc                                                                                                                          >  2:149354/1‑147 (MQ=255)
                          aCCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCga                                                                                                                        <  3:196699/149‑1 (MQ=255)
                            cgcgCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTc                                                                                                                          >  1:300145/1‑145 (MQ=255)
                            cgcgCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCgaga                                                                                                                      <  5:75956/149‑1 (MQ=255)
                              cgcAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATg                                                                                                                    <  2:98393/149‑1 (MQ=255)
                                 aaCTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGc                                                                                                                 >  6:186075/1‑149 (MQ=255)
                                   cTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCAAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGa                                                                                                               <  8:40760/149‑1 (MQ=255)
                                       tGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGc                                                                                                           <  1:149354/149‑1 (MQ=255)
                                        ggAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCa                                                                                                          >  1:218514/1‑149 (MQ=255)
                                        ggAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCa                                                                                                          >  4:279039/1‑149 (MQ=255)
                                        ggAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCa                                                                                                          >  8:277313/1‑149 (MQ=255)
                                         gAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCAc                                                                                                         <  6:304005/149‑1 (MQ=255)
                                             tCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTCTTCGAGATGGGCGATCGCACGTc                                                                                                      >  4:335114/1‑148 (MQ=255)
                                                aCGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACACTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGCTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATc                                                                                                  >  8:107832/1‑149 (MQ=255)
                                                        cTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCtta                                                                                          <  5:332769/149‑1 (MQ=255)
                                                          gAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCttac                                                                                         >  2:328833/1‑148 (MQ=255)
                                                           aaCGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCttac                                                                                         <  7:277313/147‑1 (MQ=255)
                                                              gCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCttacgatta                                                                                    <  2:300145/149‑1 (MQ=255)
                                                              gCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCttacgatta                                                                                    >  1:347914/1‑149 (MQ=255)
                                                               ccGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCttacgattac                                                                                   >  8:149331/1‑149 (MQ=255)
                                                                  aTGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGaa                                                                                <  5:301723/149‑1 (MQ=255)
                                                                   tGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGtt                                                                                                                           >  2:2805/1‑105 (MQ=255)
                                                                   tGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGtt                                                                                                                           <  1:2805/105‑1 (MQ=255)
                                                                     gCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAacac                                                                             <  6:242382/149‑1 (MQ=255)
                                                                          gTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATggg                                                                                                                  >  4:40822/1‑107 (MQ=255)
                                                                               cgccTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCa                                                                   <  2:218514/149‑1 (MQ=255)
                                                                                 ccTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCAcc                                                                 <  3:335114/149‑1 (MQ=255)
                                                                                 ccTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCAc                                                                  >  2:102649/1‑148 (MQ=255)
                                                                                     tGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCAt                                                               <  1:58319/147‑1 (MQ=255)
                                                                                       aCTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCGTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATCGGCGATCGCACGTCGATCCCCCCTTACGATTACGAACACATGCGCCTCACCATCAat                                                           >  4:28405/1‑149 (MQ=255)
                                                                                       aCTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAat                                                           >  2:281874/1‑149 (MQ=255)
                                                                                         ttCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACCCATGCGCGTCACCATCAatg                                                          >  8:332819/1‑148 (MQ=255)
                                                                                         ttCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTAGAGATGGGCGATCGCACGTCGATCACGCCTTACGATTACGAACACATGCGCGTCACCATCAatg                                                          >  5:50520/1‑148 (MQ=255)
                                                                                                 cAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGc                                                 <  2:347914/149‑1 (MQ=255)
                                                                                                       cccGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAat                                                           <  8:58529/133‑1 (MQ=255)
                                                                                                       cccGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAat                                                           >  7:58529/1‑133 (MQ=255)
                                                                                                       cccGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATTACATGACGCCACTgg                                           >  6:34627/1‑149 (MQ=255)
                                                                                                       cccGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTgg                                           <  5:186075/149‑1 (MQ=255)
                                                                                                       cccGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTgg                                           >  4:158931/1‑149 (MQ=255)
                                                                                                          gTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGgtgc                                        >  2:322794/1‑149 (MQ=255)
                                                                                                           tCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGgtgcg                                       >  3:331617/1‑149 (MQ=255)
                                                                                                             cAATGGCGCGCCGACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATTGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGGGTCACCATCAATGACATGACGCCACTGgtgcgt                                      <  5:172102/148‑1 (MQ=255)
                                                                                                             cAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGgtgcgtg                                     >  8:150056/1‑149 (MQ=255)
                                                                                                               aTGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGcc                                   <  3:279039/149‑1 (MQ=255)
                                                                                                                  gcgcgcCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTTCGAACACATGCGCGTCACCATCAATGACATTACGCCACTGGTGCGTGCCTc                                 <  5:34627/148‑1 (MQ=255)
                                                                                                                   cgcgcCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGCCTCAt                               <  7:107832/149‑1 (MQ=255)
                                                                                                                   cgcgcCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGCCTCAt                               <  1:322794/149‑1 (MQ=255)
                                                                                                                    gagcCCCCCCTGGCCCTGCGGCTGACCTGTGACGTGGAAGCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCTATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGCCTCATg                              <  8:267869/147‑1 (MQ=255)
                                                                                                                            ccTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATgcgc                                                                       >  2:94531/1‑100 (MQ=255)
                                                                                                                            ccTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATgcgc                                                                       <  1:94531/100‑1 (MQ=255)
                                                                                                                            ccTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGCCTCATGGATGCggg                      <  1:328833/149‑1 (MQ=255)
                                                                                                                                gCCCTGCTGCTGACCGGTGCCGTGGAACGGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGGCGATCCCGCCTTACGTTTACGAACACATGCGCGTCACCATCAATGACATGACGCCTCTGGTGCGTGCCTCATGGATGCGGGgt                    <  3:28405/147‑1 (MQ=255)
                                                                                                                                  ccTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGCCTCATGGATGCGGGgtgtgt                >  2:304989/1‑149 (MQ=255)
                                                                                                                                      ctgctgACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGCCTCATGGATGCGGGGTGTGTCg              <  7:332819/147‑1 (MQ=255)
                                                                                                                                         ctgACCGGTTGCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGCCTCATGGATGCGGGGTGTGTCGGCCAt         <  1:304989/149‑1 (MQ=255)
                                                                                                                                          tgACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGCCTCATGGATGCGGGGTGTGTCGGCCATg        >  4:14960/1‑149 (MQ=255)
                                                                                                                                              cGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGCCTCATGGATGCGGGGTGTGTCGGCCATGccc     >  7:114633/1‑148 (MQ=255)
                                                                                                                                                gTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGCCTCATGGATGCGGGGTGTGTCGGCCATGCCCAAc  <  8:114633/149‑1 (MQ=255)
                                                                                                                                                gTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGCCTCATGGATGCGGGGTGTGTCGGCCATGCCCAAc  <  1:281874/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AGGAACACGTGTGGGAGCCCAAGGGCACCGCGCAACTGATGGAAATCGACGGCGGCCTGAACGCCGATGGCAGCGTGGCCGCCTATGACTTCCAGACCAGCTACCCGTCCAATGGCGCGCCCACCCTGGCCCTGCTGCTGACCGGTGCCGTGGAACCGGTGCCGGCGCTGTTCGAGATGGGCGATCGCACGTCGATCCCGCCTTACGATTACGAACACATGCGCGTCACCATCAATGACATGACGCCACTGGTGCGTGCCTCATGGATGCGGGGTGTGTCGGCCATGCCCAAC  >  NC_002947/4451676‑4451968

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: ‑