Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,447,848 Δ1 bp 6.9% 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. 6.9% 224.6 / 10.4 72intergenic (‑284/‑38)nicC/nicX6‑hydroxynicotinate 3‑monooxygenase/2,5‑dihydroxypyridine 5,6‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base . (2/3);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACC  >  NC_002947/4447703‑4447994
                                                                                                                                                  |                                                                                                                                                  
gATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATg                                                                                                                                                  <  8:284683/148‑1 (MQ=255)
   tGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTgc                                                                                                                                              >  2:347244/1‑149 (MQ=255)
    gTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTgcg                                                                                                                                             >  6:31536/1‑149 (MQ=255)
          gTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTg                                                                                                                                       <  2:205305/149‑1 (MQ=255)
               taatgttCAATGTTTGGATAGGTAAAGTTTTTACAGATCTCTTTCGCTCGCTGTCAGGATGCGTTTTCATATGCACAAGTAGGGTGTACACTCTATTATTTTGCACCTTGATGCAGCACCATGATCTGTCATGGTGGGGCGCTGTCaag                                                                                                                                  <  8:287511/143‑2 (MQ=255)
               tACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAAc                                                                                                                                  >  3:345828/1‑149 (MQ=255)
                  tGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGc                                                                                                                                >  1:266000/1‑148 (MQ=255)
                    gTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCAc                                                                                                                             >  3:9671/1‑149 (MQ=255)
                      cAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCAccc                                                                                                                           <  3:192538/149‑1 (MQ=255)
                         tgtgtgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAg                                                                                                                         >  4:191312/1‑148 (MQ=255)
                           tgtgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAggcc                                                                                                                      >  8:146052/1‑147 (MQ=255)
                                aTAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAAccc                                                                                                                 <  1:347244/149‑1 (MQ=255)
                                  aGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGcc                                                                                                                               <  5:281475/133‑1 (MQ=255)
                                  aGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGcc                                                                                                                               >  6:281475/1‑133 (MQ=255)
                                    gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTgcgc                                                                                                                                            <  2:140915/118‑1 (MQ=255)
                                         gTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCgg                                                                                                        >  1:204506/1‑149 (MQ=255)
                                         gTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCACGAGAACCCCGATGCCgg                                                                                                        >  5:39312/1‑149 (MQ=255)
                                           ttGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgct                                                                                                                                             >  5:120308/1‑109 (MQ=255)
                                           ttGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgct                                                                                                                                             <  6:120308/110‑2 (MQ=255)
                                                cAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCaa                                                                                                 >  8:306555/1‑149 (MQ=255)
                                                 aGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACg                                                                                                                                 <  8:374410/116‑1 (MQ=255)
                                                 aGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACg                                                                                                                                 >  7:374410/1‑116 (MQ=255)
                                                  gATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATg                                                                                               <  7:306555/149‑1 (MQ=255)
                                                  gATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATg                                                                                               <  2:266000/149‑1 (MQ=255)
                                                       cTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACaa                                                                                          <  1:25494/149‑1 (MQ=255)
                                                        ttGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAAc                                                                                         >  7:388724/1‑149 (MQ=255)
                                                            gctcgctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGAc                                                                                     >  1:382680/1‑149 (MQ=255)
                                                              tcgctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGAcc                                                                                    >  5:392366/1‑148 (MQ=255)
                                                                gctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTAATATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGCTCAGCAATGCACAACTGAcc                                                                                    >  8:353694/1‑146 (MQ=255)
                                                                      aGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGa                                                                                      <  2:4350/138‑1 (MQ=255)
                                                                      aGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGa                                                                                      >  1:4350/1‑138 (MQ=255)
                                                                        gATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTc                                                                           >  6:334139/1‑147 (MQ=255)
                                                                         aTGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGTCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATTTTCGAg                                                                        <  2:204506/149‑1 (MQ=255)
                                                                         aTGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAg                                                                        <  6:39312/149‑1 (MQ=255)
                                                                             gTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACg                                                                    >  6:270358/1‑149 (MQ=255)
                                                                              ttttCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGt                                                                   <  5:334139/149‑1 (MQ=255)
                                                                                  cATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACt                                                                <  2:382680/148‑1 (MQ=255)
                                                                                   aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTg                                                               >  5:91740/1‑148 (MQ=255)
                                                                                   aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTg                                                               <  3:90362/148‑1 (MQ=255)
                                                                                      tGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGc                                                           >  8:226532/1‑149 (MQ=255)
                                                                                        cacaAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGaa                                                             <  6:192359/145‑1 (MQ=255)
                                                                                        cacaAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGaa                                                             <  1:395043/145‑1 (MQ=255)
                                                                                        cacaAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAg                                                            <  3:191312/146‑1 (MQ=255)
                                                                                           aaGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCt                                                                                                                                        >  8:166601/1‑67 (MQ=255)
                                                                                           aaGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCt                                                                                                                                        <  7:166601/67‑1 (MQ=255)
                                                                                           aaGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGtct                                                      <  6:392366/149‑1 (MQ=255)
                                                                                             gTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGtctc                                                     >  3:236961/1‑148 (MQ=255)
                                                                                                    gtACACAGTATTATTTTGCACCGATATGCAGCAACAGTCACTGCCATGGTGCGCCCCTTTCAACGGCACCCAGGAGAAGCCCGATGACGGTTAGCAATGCACAACTGACCCATTTGTTCGATCACGTACTGAAGCTGTCTCGTGTGGa                                              <  4:345828/148‑1 (MQ=255)
                                                                                                       cacacTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCCATCCCGGTGAGCAATGCACAACTGACCCAGATGATCGAGCACGTACTGAAGCTGTCTCGAGTGGACGa                                           >  4:300263/1‑148 (MQ=255)
                                                                                                        acacTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGaaa                                         <  4:9671/149‑1 (MQ=255)
                                                                                                            tattatTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGaa                                                             <  4:340572/125‑1 (MQ=255)
                                                                                                            tattatTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGaa                                                             >  3:340572/1‑125 (MQ=255)
                                                                                                               tattTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACGGAAGCTGTCTCGAGTGGACGAAACGCagag                                  >  5:322849/1‑149 (MQ=255)
                                                                                                                      cACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTg                                                                                       <  5:104644/89‑1 (MQ=255)
                                                                                                                      cACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTg                                                                                       >  6:104644/1‑89 (MQ=255)
                                                                                                                       aCCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGAccc                                                                                   <  1:140232/92‑1 (MQ=255)
                                                                                                                       aCCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGAccc                                                                                   >  2:140232/1‑92 (MQ=255)
                                                                                                                       aCCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCg                          >  7:169287/1‑149 (MQ=255)
                                                                                                                       aCCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCg                          >  5:268101/1‑149 (MQ=255)
                                                                                                                        ccGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGt                         >  3:261861/1‑149 (MQ=255)
                                                                                                                              gggcgcaCCATGGCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGaa                                                             <  4:246913/103‑1 (MQ=255)
                                                                                                                                    accatggCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAagag                >  2:360658/6‑146 (MQ=255)
                                                                                                                                    accatggCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAagag                <  1:360658/141‑1 (MQ=255)
                                                                                                                                     cccctgCACTGCCATGGTGCGCCCCTGTCAACGCCACCCACGAGAAACCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCa             <  8:388724/144‑1 (MQ=255)
                                                                                                                                          gCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGt                                                                   <  7:24764/90‑1 (MQ=255)
                                                                                                                                          gCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGt                                                                   >  8:24764/1‑90 (MQ=255)
                                                                                                                                                ccATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGaa                                                             <  4:256307/90‑1 (MQ=255)
                                                                                                                                                ccATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGaa                                                             >  3:256307/1‑90 (MQ=255)
                                                                                                                                                ccATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCactac         <  2:124564/142‑1 (MQ=255)
                                                                                                                                                ccATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCactac         >  1:124564/1‑142 (MQ=255)
                                                                                                                                                ccATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGAcc  >  8:279491/1‑149 (MQ=255)
                                                                                                                                                ccATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGAc   <  5:31536/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACC  >  NC_002947/4447703‑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: