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

TGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACCCG  >  NC_002947/4447702‑4447996
                                                                                                                                                   |                                                                                                                                                      
tGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATg                                                                                                                                                      >  2:351569/1‑149 (MQ=255)
  aTTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGt                                                                                                                                                    >  3:396886/1‑149 (MQ=255)
   ttGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACTTTTGAACAGATCACTTGCGCTCTCTTGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACAATTATTTTGCACAGAGATGCAGCACCATGGCAGTGCCTGGTg                                                                                                                                                   <  4:415423/149‑1 (MQ=255)
     gTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgct                                                                                                                                                 >  1:53321/1‑148 (MQ=255)
      ttACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgctg                                                                                                                                                <  8:303497/149‑3 (MQ=255)
        aCCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGcc                                                                                                                                               <  6:360903/148‑1 (MQ=255)
         ccGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTgc                                                                                                                                                  <  1:106133/144‑1 (MQ=255)
         ccGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTgc                                                                                                                                                  >  2:106133/1‑144 (MQ=255)
         ccGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGccc                                                                                                                                              <  3:77602/148‑1 (MQ=255)
            gaaaTACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGGACAAGTAGGGTGTACACACTATTATTTTGCACCGAGCTGCAGCACCAGGCACTGCCATGGTGCGcccgggt                                                                                                                                          >  8:127445/2‑145 (MQ=255)
                  cTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGc                                                                                                                                    <  5:274385/149‑1 (MQ=255)
                   tGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGc                                                                                                                                    >  4:415221/1‑148 (MQ=255)
                         aTGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAg                                                                                                                             <  5:335900/149‑1 (MQ=255)
                          tgtgtgGATAGGTAACGTTTGCACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCa                                                                                                                              <  8:185884/147‑1 (MQ=255)
                          tgtgtgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAgg                                                                                                                            <  7:20679/149‑1 (MQ=255)
                            tgtgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGga                                                                                                                           <  5:327399/148‑1 (MQ=255)
                            tgtgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGga                                                                                                                           <  7:330294/148‑1 (MQ=255)
                              tgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGaa                                                                                                                        <  6:196879/149‑1 (MQ=255)
                                     gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑At                                                                                                               <  4:396886/149‑1 (MQ=255)
                                              gAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTACGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGa                                                                                                       <  2:66325/148‑1 (MQ=255)
                                               aaCAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAg                                                                                                      >  8:277108/1‑148 (MQ=255)
                                               aaCAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGGGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGc                                                                                                     >  8:309325/1‑149 (MQ=255)
                                                     tCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCGCCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACACCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGca                                                                                               >  2:125272/1‑149 (MQ=255)
                                                     tCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGca                                                                                               <  1:351569/149‑1 (MQ=255)
                                                       aCTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGAACCAGGCACTGCAATGGTGCGCCCCTGTCAACTCCACCAAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGcac                                                                                              >  1:388128/1‑148 (MQ=255)
                                                               tcgctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGc                                                                                                             >  1:83189/1‑125 (MQ=255)
                                                               tcgctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGc                                                                                                             <  2:83189/125‑1 (MQ=255)
                                                               tcgctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGAcc                                                                                      >  7:106887/1‑148 (MQ=255)
                                                                 gctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGAccc                                                                                     <  1:265480/147‑1 (MQ=255)
                                                                       aGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTc                                                                             <  8:106887/149‑1 (MQ=255)
                                                                         gATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTTCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCATATGTTCGa                                                                           <  8:218844/149‑1 (MQ=255)
                                                                             ccttttCATCTGCACAAGTAGGGTGTACACACAATTATTTTGCACCGAGATTCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCGGG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCa                                                                        <  4:42016/146‑1 (MQ=255)
                                                                               ttttCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACg                                                                      >  1:295635/1‑148 (MQ=255)
                                                                                    aTCTGCACAAGTAGGGTGTACCCACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTg                                                                 >  5:155947/1‑148 (MQ=255)
                                                                                    aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGTCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTg                                                                 >  7:219604/1‑148 (MQ=255)
                                                                                     tCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTg                                                                 >  8:212173/1‑147 (MQ=255)
                                                                                       tGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGc                                                             >  7:46580/1‑149 (MQ=255)
                                                                                         cacaAGTCGGGGGTACACACTATTATTTTGCACCGAGATGCAGCACCAGTCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGaa                                                               <  1:346567/145‑1 (MQ=255)
                                                                                            aaGTAGGGTGTACACACTATTATTTTGCTCCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGtct                                                        >  1:26433/1‑149 (MQ=255)
                                                                                            aaGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGa                                                                                                       >  4:302266/1‑102 (MQ=255)
                                                                                            aaGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGa                                                                                                       <  3:302266/102‑1 (MQ=255)
                                                                                                   gtgtACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGTGAACCCCGCCATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGt                                                   <  1:125632/149‑1 (MQ=255)
                                                                                                   gtgtACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTgg                                                 <  7:277108/149‑1 (MQ=255)
                                                                                                    tgtACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                                <  8:46580/149‑1 (MQ=255)
                                                                                                     gtACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGAc                                               >  1:69919/1‑149 (MQ=255)
                                                                                                         acacTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGTCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGaaa                                           <  8:219604/149‑1 (MQ=255)
                                                                                                          aacTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAAc                                          <  7:127445/148‑1 (MQ=255)
                                                                                                           acTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACg                                         >  5:251159/1‑149 (MQ=255)
                                                                                                                        aCCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGt                                                   <  8:355165/126‑1 (MQ=255)
                                                                                                                        aCCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGt                                                   >  7:355165/1‑126 (MQ=255)
                                                                                                                          cGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                                <  3:143110/127‑1 (MQ=255)
                                                                                                                          cGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                                >  4:143110/1‑127 (MQ=255)
                                                                                                                          cGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTg                          >  6:333858/1‑149 (MQ=255)
                                                                                                                            agaTGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCt                        >  4:104919/1‑149 (MQ=255)
                                                                                                                               tGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTc                       >  7:150929/1‑147 (MQ=255)
                                                                                                                                ggcgcaCCATGGCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTc                       >  6:418036/4‑146 (MQ=255)
                                                                                                                                    caccatggCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGc                 >  8:29892/7‑148 (MQ=255)
                                                                                                                                      ccAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCac              <  3:415221/149‑1 (MQ=255)
                                                                                                                                          ggCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCactac           >  5:221049/1‑149 (MQ=255)
                                                                                                                                           gCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCactact          >  3:161805/1‑149 (MQ=255)
                                                                                                                                            cACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTc         >  7:188224/1‑149 (MQ=255)
                                                                                                                                             aCTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTcc        >  8:362704/1‑149 (MQ=255)
                                                                                                                                                  cATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGAccc   >  5:161421/1‑149 (MQ=255)
                                                                                                                                                   aTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACCcg  >  3:253502/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACCCG  >  NC_002947/4447702‑4447996

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: