Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,447,848 Δ1 bp 7.2% intergenic (‑284/‑38) nicC ← / → nicX 6‑hydroxynicotinate 3‑monooxygenase/2,5‑dihydroxypyridine 5,6‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,447,8480A. 7.2% 219.3 / 10.7 69intergenic (‑284/‑38)nicC/nicX6‑hydroxynicotinate 3‑monooxygenase/2,5‑dihydroxypyridine 5,6‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (30/34);  new base . (3/2);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACC  >  NC_002947/4447705‑4447994
                                                                                                                                                    |                                                                                                                                                  
ttGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTg                                                                                                                                               >  2:294473/1‑149 (MQ=255)
  gTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTgcg                                                                                                                                             >  5:223137/1‑149 (MQ=255)
   ttACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTgcgc                                                                                                                                            <  2:196660/149‑1 (MQ=255)
      ccGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGccc                                                                                                                                          <  5:86467/148‑1 (MQ=255)
       cGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCt                                                                                                                                        <  4:256920/149‑1 (MQ=255)
       cGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCt                                                                                                                                        >  1:247902/1‑149 (MQ=255)
        gTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTg                                                                                                                                       >  6:34717/1‑149 (MQ=255)
           aaTACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCATGGCAGTGCCTGGTgctgc                                                                                                                                           >  4:121071/1‑139 (MQ=255)
           aaTACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCATGGCAGTGCCTGGTgctgc                                                                                                                                           <  3:121071/142‑4 (MQ=255)
             tACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAAc                                                                                                                                  <  2:23875/149‑1 (MQ=255)
               cTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCATGCACCAGGCACTGCCATGGTGCGCCCCTGTCAAc                                                                                                                                  <  1:206610/148‑1 (MQ=255)
                         tgtgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGga                                                                                                                       <  2:247902/148‑1 (MQ=255)
                          gtgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGgag                                                                                                                      >  4:312362/1‑148 (MQ=255)
                              aTAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAAccc                                                                                                                 >  8:249236/1‑149 (MQ=255)
                              aTAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAAccc                                                                                                                 >  2:120653/1‑149 (MQ=255)
                                 ggTAACGTTTGAACAGATCACTTGCGCTCTCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGTCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGa                                                                                                              <  5:98269/149‑1 (MQ=255)
                                  gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCCCCCAGGAGAACCCCGAt                                                                                                             >  3:46154/1‑149 (MQ=255)
                                  gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGAt                                                                                                             <  2:137813/149‑1 (MQ=255)
                                  gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGAt                                                                                                             >  5:47010/1‑149 (MQ=255)
                                  gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGAt                                                                                                             <  6:223137/149‑1 (MQ=255)
                                  gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGAt                                                                                                             <  1:26246/149‑1 (MQ=255)
                                           gAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGa                                                                                                     <  2:223807/148‑1 (MQ=255)
                                            aaCAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAg                                                                                                    >  5:126164/1‑148 (MQ=255)
                                                gATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATg                                                                                               <  6:47010/149‑1 (MQ=255)
                                                gATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATg                                                                                               <  4:46154/149‑1 (MQ=255)
                                                 aTCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATg                                                                                               >  8:206479/1‑148 (MQ=255)
                                                   cACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGca                                                                                             >  2:154083/1‑148 (MQ=255)
                                                      ttGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAAc                                                                                         >  7:183394/1‑149 (MQ=255)
                                                      ttGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACACAGGAGAACCCCGATGCCGGTGAGCAATGCACaa                                                                                          >  8:55904/1‑148 (MQ=255)
                                                        gcgcTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATAGTGGAGCACCATGCACTTCCATTGTTCGCCCCAGTCAACCCCACCCAAGAGAACCCCGATGCCCGTGAGCAAAGcaca                                                                                           >  1:100533/1‑148 (MQ=255)
                                                        gcgcTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAAc                                                                                         <  1:294473/147‑1 (MQ=255)
                                                            tcgctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGAcc                                                                                    >  2:288455/1‑148 (MQ=255)
                                                                       aTGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAg                                                                        >  6:11913/1‑149 (MQ=255)
                                                                            ttttCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACg                                                                    >  3:63436/1‑148 (MQ=255)
                                                                               tCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTAc                                                                 <  5:34717/148‑1 (MQ=255)
                                                                               tCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACt                                                                >  7:178672/1‑149 (MQ=255)
                                                                                 aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTg                                                               >  5:232508/1‑148 (MQ=255)
                                                                                      cacaAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCt                                                          <  7:55904/148‑1 (MQ=255)
                                                                                      cacaAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCt                                                          <  7:249236/148‑1 (MQ=255)
                                                                                      cacaAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCt                                                          <  6:232508/148‑1 (MQ=255)
                                                                                       acaAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGt                                                        <  3:312362/149‑1 (MQ=255)
                                                                                         aaGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAg                                                                        >  5:186312/1‑131 (MQ=255)
                                                                                         aaGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAg                                                                        <  6:186312/131‑1 (MQ=255)
                                                                                         aaGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGCTCGAGCACGTACTGAAGCTGtct                                                      <  8:183394/149‑1 (MQ=255)
                                                                                          aGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGtctc                                                     <  1:288455/149‑1 (MQ=255)
                                                                                           gTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGtctc                                                     >  1:32592/1‑148 (MQ=255)
                                                                                              gggTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGt                                                 <  3:312316/149‑1 (MQ=255)
                                                                                                 tgtACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                              <  7:206479/149‑1 (MQ=255)
                                                                                                 tgtACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                              <  6:126164/149‑1 (MQ=255)
                                                                                                   tACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGAc                                             <  2:32592/148‑1 (MQ=255)
                                                                                                      acacTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGca                                                                                             <  6:185509/97‑1 (MQ=255)
                                                                                                      acacTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGca                                                                                             >  5:185509/1‑97 (MQ=255)
                                                                                                          tattatTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACg                                       >  2:154843/1‑147 (MQ=255)
                                                                                                             tattTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGc                                                           >  4:240706/1‑124 (MQ=255)
                                                                                                             tattTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGc                                                           <  3:240706/124‑1 (MQ=255)
                                                                                                              attTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGc                                 >  8:154706/1‑149 (MQ=255)
                                                                                                                   gCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTgg                             <  1:120653/148‑1 (MQ=255)
                                                                                                                     aCCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAg                                                                        >  4:18536/1‑103 (MQ=255)
                                                                                                                     aCCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAg                                                                        <  3:18536/103‑1 (MQ=255)
                                                                                                                      ccGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGt                         <  7:154706/149‑1 (MQ=255)
                                                                                                                        gaga‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGc                       <  8:178672/149‑1 (MQ=255)
                                                                                                                        gaga‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGc                       <  1:154083/149‑1 (MQ=255)
                                                                                                                                   cgcaCCATGGCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAAACCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAagag                >  3:23245/2‑149 (MQ=255)
                                                                                                                                    gcaCCATGGCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGa                                                                                      <  6:149056/78‑1 (MQ=38)
                                                                                                                                    gcaCCATGGCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCACGATGCCGGTGAGCAATGCACAACTGa                                                                                      >  5:149056/1‑78 (MQ=37)
                                                                                                                                     caCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGcc              >  5:270984/1‑149 (MQ=255)
                                                                                                                                         aGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCacta          <  1:154843/149‑1 (MQ=255)
                                                                                                                                            gCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCactact        >  1:32517/1‑149 (MQ=255)
                                                                                                                                                  ccATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGAcc  >  7:268813/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACC  >  NC_002947/4447705‑4447994

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: