Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,447,892 Δ1 bp 9.7% 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,8920A. 9.7% 179.8 / 13.2 62intergenic (‑284/‑38)nicC/nicX6‑hydroxynicotinate 3‑monooxygenase/2,5‑dihydroxypyridine 5,6‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (25/31);  new base . (4/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

GATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACCCGC  >  NC_002947/4447747‑4448041
                                                                                                                                                  |                                                                                                                                                     
gATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGt                                                                                                                                                   <  2:1618855/150‑1 (MQ=255)
gATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGt                                                                                                                                                   <  2:1101571/150‑1 (MQ=255)
   tGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgct                                                                                                                                                <  1:559237/150‑2 (MQ=255)
   tGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTgcg                                                                                                                                                >  1:1877613/1‑150 (MQ=255)
    gTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgctg                                                                                                                                               <  1:84310/150‑3 (MQ=255)
          gTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGt                                                                                                                                         <  2:1394543/150‑1 (MQ=255)
             aaTACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCaa                                                                                                                                      <  1:867496/150‑1 (MQ=255)
              aTACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAAc                                                                                                                                     <  1:9201/150‑1 (MQ=255)
                aCTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGc                                                                                                                                   <  1:1567439/150‑1 (MQ=255)
                  tGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCa                                                                                                                                 <  1:464279/150‑1 (MQ=255)
                         tgtgtgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGga                                                                                                                          <  2:79971/150‑1 (MQ=255)
                         tgtgtgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCACCACCAGGCACTGCCATGGTGCTCCCCTGTCATCGCTACCCAGga                                                                                                                          <  2:1755865/150‑1 (MQ=255)
                          gtgtgGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGgag                                                                                                                         <  2:418792/150‑1 (MQ=255)
                                    gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAAc                                                                                                                                     <  2:1622595/128‑1 (MQ=255)
                                    gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAAc                                                                                                                                     >  1:1622595/1‑128 (MQ=255)
                                    gTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATg                                                                                                               >  1:155682/1‑150 (MQ=255)
                                       aCGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCg                                                                                                            <  2:155682/150‑1 (MQ=255)
                                       aCGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCg                                                                                                            <  2:291029/150‑1 (MQ=255)
                                         gTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGt                                                                                                          >  2:1320324/1‑150 (MQ=255)
                                           ttGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGa                                                                                                        <  2:1436796/150‑1 (MQ=255)
                                            tGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAg                                                                                                       <  1:21692/150‑1 (MQ=255)
                                                   aTCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGca                                                                                                >  1:202598/1‑150 (MQ=255)
                                                    tCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGcac                                                                                               >  2:1416351/1‑150 (MQ=255)
                                                      aCTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACaa                                                                                             <  2:202598/150‑1 (MQ=255)
                                                      aCTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACaa                                                                                             <  2:1877613/150‑1 (MQ=255)
                                                                  tGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGCGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGa                                                                                                                 >  1:1411043/1‑118 (MQ=255)
                                                                  tGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGCGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGa                                                                                                                 <  2:1411043/118‑1 (MQ=255)
                                                                     cAGGATGCGTTTTGATCTGCACAAGTCGGGTGTACACATTATTATTCTGCACCGGGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCCATGCACAACTGACCCAGATGTTc                                                                              <  2:815628/150‑1 (MQ=255)
                                                                     cAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTc                                                                              >  2:974104/1‑150 (MQ=255)
                                                                        gATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAg                                                                           <  1:1888817/150‑1 (MQ=255)
                                                                                   aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGTACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGa                                                                                                        >  2:441531/1‑110 (MQ=255)
                                                                                   aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGTACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGa                                                                                                        <  1:441531/110‑1 (MQ=255)
                                                                                    tCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAg                                                               >  2:934159/1‑150 (MQ=255)
                                                                                    tCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCAAAACTGACCCAGATGTTCGAGCACGTACTGAAg                                                               >  2:1447459/1‑150 (MQ=255)
                                                                                         acaAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGAGCCAGATGTTCGAGCACGTACTGAAGCTGtc                                                          >  2:1604331/1‑150 (MQ=255)
                                                                                               aGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACACCTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGt                                                    <  1:240550/150‑1 (MQ=255)
                                                                                                       cacacTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGa                                                                                                        >  1:575946/1‑90 (MQ=255)
                                                                                                       cacacTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGa                                                                                                        <  2:575946/90‑1 (MQ=255)
                                                                                                          acTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAgg                                                                                                                           >  2:201039/1‑68 (MQ=255)
                                                                                                          acTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAgg                                                                                                                           <  1:201039/68‑1 (MQ=255)
                                                                                                           cTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAGCCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCa                                        <  1:974104/150‑1 (MQ=255)
                                                                                                           cTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCa                                        <  1:1320324/150‑1 (MQ=255)
                                                                                                           cTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCAACCACGAGTCCCATGGTATCGGTGAGCAATGGTGTTCTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCa                                        <  1:1178322/150‑1 (MQ=255)
                                                                                                                  tttGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTATTGAAGCTGTCTCGAGTGGACGAAAAGCAGAGCGTg                                 <  2:1014065/150‑1 (MQ=255)
                                                                                                                     gCACCGAGATGCAGCACCAGCCAGTGCCATGGTGCGCCCCTGTCAACGCCACCCAAGAAAACAACATTGCCGTTGAGCAATGCTCAACTGACACAGATGGTCGAGCACGTACTGATGCTGTCTCCAGTTGACGAAACGCAGAGCGTGGcc                              <  2:826758/150‑1 (MQ=255)
                                                                                                                       aCCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGt                            >  1:251281/1‑150 (MQ=255)
                                                                                                                        ccGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAGGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTg                           >  2:40547/1‑150 (MQ=255)
                                                                                                                         cGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGc                          >  1:860371/1‑150 (MQ=255)
                                                                                                                         cGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACGGAAGCTGTCTCGAGAGGACGAAACGCAGAGCGTGGCCGtcc                          >  2:36335/1‑148 (MQ=255)
                                                                                                                             ggggcgcaCCATGGCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTAGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCaa                      >  2:916604/6‑150 (MQ=255)
                                                                                                                              gggcgcaCCATGGCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCaaa                     >  2:1728328/5‑149 (MQ=255)
                                                                                                                                  gcaCCATGGCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGcc                 >  2:1790979/1‑150 (MQ=255)
                                                                                                                                  gcaCCATGGCAGTGCCTGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGcc                 >  2:621720/1‑150 (MQ=255)
                                                                                                                                         ggCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCactacc           >  2:1509747/1‑149 (MQ=255)
                                                                                                                                            aCTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCg        >  1:1645407/1‑150 (MQ=255)
                                                                                                                                            aCTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCg        >  1:1645214/1‑150 (MQ=255)
                                                                                                                                              tGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCaa                      <  2:1389332/134‑1 (MQ=255)
                                                                                                                                              tGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCaa                      >  1:1389332/1‑134 (MQ=255)
                                                                                                                                                ccATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGAccc    <  1:1416351/150‑1 (MQ=255)
                                                                                                                                                ccATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGAccc    >  1:1399039/1‑150 (MQ=255)
                                                                                                                                                 cATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACCcg   >  2:1412614/1‑150 (MQ=255)
                                                                                                                                                 cATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACCcg   >  1:1036348/1‑150 (MQ=255)
                                                                                                                                                  aTGGTGCGCCCCTGTCAACGCCACCCAGTAGAACCCCGATGCCGGTGAGCAATGCACAACTGACTCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACCcgc  <  2:722975/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTGACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACCCGC  >  NC_002947/4447747‑4448041

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: