Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,447,843 C→G 8.1% 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 8.1% 193.7 / 6.3 62intergenic (‑279/‑43)nicC/nicX6‑hydroxynicotinate 3‑monooxygenase/2,5‑dihydroxypyridine 5,6‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base G (3/2);  total (32/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.56e-01

ATCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCC  >  NC_002947/4447695‑4447990
                                                                                                                                                     |                                                                                                                                                     
aTCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCAc                                                                                                                                                       <  5:422451/149‑1 (MQ=255)
       tGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATg                                                                                                                                                >  2:351569/1‑149 (MQ=255)
         aTTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGt                                                                                                                                              >  3:396886/1‑149 (MQ=255)
          ttGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACTTTTGAACAGATCACTTGCGCTCTCTTGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACAATTATTTTGCACAGAGATGCAGCACCATGGCAGTGCC‑TGGTg                                                                                                                                             <  4:415423/149‑1 (MQ=255)
            gTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCC‑TGGTgct                                                                                                                                           >  1:53321/1‑148 (MQ=255)
             ttACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCC‑TGGTgctg                                                                                                                                          <  8:303497/149‑3 (MQ=255)
               aCCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGcc                                                                                                                                         <  6:360903/148‑1 (MQ=255)
                ccGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTgc                                                                                                                                            <  1:106133/144‑1 (MQ=255)
                ccGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTgc                                                                                                                                            >  2:106133/1‑144 (MQ=255)
                ccGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGccc                                                                                                                                        <  3:77602/148‑1 (MQ=255)
                   gaaaTACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGGACAAGTAGGGTGTACACACTATTATTTTGCACCGAGCTGCAGCACCA‑GGCACTGCCATGGTGCGcccgggt                                                                                                                                    >  8:127445/2‑145 (MQ=255)
                         cTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGc                                                                                                                              <  5:274385/149‑1 (MQ=255)
                          tGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGc                                                                                                                              >  4:415221/1‑148 (MQ=255)
                                aTGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAg                                                                                                                       <  5:335900/149‑1 (MQ=255)
                                 tgtgtgGATAGGTAACGTTTGCACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCa                                                                                                                        <  8:185884/147‑1 (MQ=255)
                                 tgtgtgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAgg                                                                                                                      <  7:20679/149‑1 (MQ=255)
                                   tgtgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGga                                                                                                                     <  5:327399/148‑1 (MQ=255)
                                   tgtgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGga                                                                                                                     <  7:330294/148‑1 (MQ=255)
                                     tgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGaa                                                                                                                  <  6:196879/149‑1 (MQ=255)
                                            gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑At                                                                                                         <  4:396886/149‑1 (MQ=255)
                                                     gAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTACGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGa                                                                                                 <  2:66325/148‑1 (MQ=255)
                                                      aaCAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAg                                                                                                >  8:277108/1‑148 (MQ=255)
                                                      aaCAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGGGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGc                                                                                               >  8:309325/1‑149 (MQ=255)
                                                            tCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCGCCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACACCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGca                                                                                         >  2:125272/1‑149 (MQ=255)
                                                            tCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGca                                                                                         <  1:351569/149‑1 (MQ=255)
                                                              aCTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGAACCA‑GGCACTGCAATGGTGCGCCCCTGTCAACTCCACCAAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGcac                                                                                        >  1:388128/1‑148 (MQ=255)
                                                                      tcgctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGc                                                                                                       <  2:83189/125‑1 (MQ=255)
                                                                      tcgctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGc                                                                                                       >  1:83189/1‑125 (MQ=255)
                                                                      tcgctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGAcc                                                                                >  7:106887/1‑148 (MQ=255)
                                                                        gctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGAccc                                                                               <  1:265480/147‑1 (MQ=255)
                                                                              aGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTc                                                                       <  8:106887/149‑1 (MQ=255)
                                                                                gATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTTCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCATATGTTCGa                                                                     <  8:218844/149‑1 (MQ=255)
                                                                                    ccttttCATCTGCACAAGTAGGGTGTACACACAATTATTTTGCACCGAGATTCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCGGG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCa                                                                  <  4:42016/146‑1 (MQ=255)
                                                                                      ttttCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACg                                                                >  1:295635/1‑148 (MQ=255)
                                                                                           aTCTGCACAAGTAGGGTGTACCCACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTg                                                           >  5:155947/1‑148 (MQ=255)
                                                                                           aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGTCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTg                                                           >  7:219604/1‑148 (MQ=255)
                                                                                            tCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTg                                                           >  8:212173/1‑147 (MQ=255)
                                                                                              tGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGc                                                       >  7:46580/1‑149 (MQ=255)
                                                                                                cacaAGTCGGGGGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GTCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGaa                                                         <  1:346567/145‑1 (MQ=255)
                                                                                                   aaGTAGGGTGTACACACTATTATTTTGCTCCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGtct                                                  >  1:26433/1‑149 (MQ=255)
                                                                                                   aaGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGa                                                                                                 >  4:302266/1‑102 (MQ=255)
                                                                                                   aaGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGa                                                                                                 <  3:302266/102‑1 (MQ=255)
                                                                                                          gtgtACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGTGAACCCCGCCATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGt                                             <  1:125632/149‑1 (MQ=255)
                                                                                                          gtgtACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTgg                                           <  7:277108/149‑1 (MQ=255)
                                                                                                           tgtACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                          <  8:46580/149‑1 (MQ=255)
                                                                                                            gtACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGAc                                         >  1:69919/1‑149 (MQ=255)
                                                                                                                acacTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGTCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGaaa                                     <  8:219604/149‑1 (MQ=255)
                                                                                                                 aacTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAAc                                    <  7:127445/148‑1 (MQ=255)
                                                                                                                  acTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACg                                   >  5:251159/1‑149 (MQ=255)
                                                                                                                               aCCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGt                                             >  7:355165/1‑126 (MQ=255)
                                                                                                                               aCCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGt                                             <  8:355165/126‑1 (MQ=255)
                                                                                                                                 cGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                          <  3:143110/127‑1 (MQ=255)
                                                                                                                                 cGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                          >  4:143110/1‑127 (MQ=255)
                                                                                                                                 cGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTg                    >  6:333858/1‑149 (MQ=255)
                                                                                                                                   agaTGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCt                  >  4:104919/1‑149 (MQ=255)
                                                                                                                                      tGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTc                 >  7:150929/1‑147 (MQ=255)
                                                                                                                                       ggcgcaCCATGGCAGTGCC‑TGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTc                 >  6:418036/4‑146 (MQ=255)
                                                                                                                                           caccatggCAGTGCC‑TGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGc           >  8:29892/7‑148 (MQ=255)
                                                                                                                                             ccA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑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)
                                                                                                                                                     |                                                                                                                                                     
ATCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCA‑GGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCG‑‑ATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCC  >  NC_002947/4447695‑4447990

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: ‑