Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,447,843 C→G 7.2% intergenic (‑279/‑43) 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,8430CG 7.2% 218.7 / 6.3 69intergenic (‑279/‑43)nicC/nicX6‑hydroxynicotinate 3‑monooxygenase/2,5‑dihydroxypyridine 5,6‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base G (3/2);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACT  >  NC_002947/4447696‑4447988
                                                                                                                                                        |                                                                                                                                                 
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                                                                                                       <  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                                        <  6:126164/149‑1 (MQ=255)
                                                                                                          tgtACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                        <  7:206479/149‑1 (MQ=255)
                                                                                                            tACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGAc                                       <  2:32592/148‑1 (MQ=255)
                                                                                                               acacTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGca                                                                                       >  5:185509/1‑97 (MQ=255)
                                                                                                               acacTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGca                                                                                       <  6:185509/97‑1 (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                                                                  <  3:18536/103‑1 (MQ=255)
                                                                                                                              aCCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAg                                                                  >  4:18536/1‑103 (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)
                                                                                                                                                        |                                                                                                                                                 
TCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGA‑‑‑TGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACT  >  NC_002947/4447696‑4447988

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: