Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,447,838:1 +T 7.5% intergenic (‑274/‑48) 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,8381.T 7.5% 262.4 / 15.5 67intergenic (‑274/‑48)nicC/nicX6‑hydroxynicotinate 3‑monooxygenase/2,5‑dihydroxypyridine 5,6‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base . (28/34);  new base T (3/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

CGAAATCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTA  >  NC_002947/4447691‑4447986
                                                                                                                                                        |                                                                                                                                                    
cgaaaTCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAg                                                                                                                                                     <  6:275387/149‑1 (MQ=255)
 gaaaTCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAgg                                                                                                                                                    <  6:28786/149‑1 (MQ=255)
     tCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCTCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACt                                                                                                                                                <  2:286618/149‑1 (MQ=255)
              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                                                                                                     <  1:26246/149‑1 (MQ=255)
                                                gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGAt                                                                                                     >  5:47010/1‑149 (MQ=255)
                                                gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGAt                                                                                                     <  6:223137/149‑1 (MQ=255)
                                                gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGAt                                                                                                     <  2:137813/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                                                                <  6:186312/131‑1 (MQ=255)
                                                                                                       aaGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAg                                                                >  5:186312/1‑131 (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                                      <  6:126164/149‑1 (MQ=255)
                                                                                                               tgtACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                      <  7:206479/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                                                   <  3:240706/124‑1 (MQ=255)
                                                                                                                           tattTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGc                                                   >  4:240706/1‑124 (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)
                                                                                                                                                        |                                                                                                                                                    
CGAAATCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTA  >  NC_002947/4447691‑4447986

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: