Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,447,892 Δ1 bp 9.0% 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.0% 261.3 / 20.8 89intergenic (‑284/‑38)nicC/nicX6‑hydroxynicotinate 3‑monooxygenase/2,5‑dihydroxypyridine 5,6‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (46/35);  new base . (2/6);  total (48/41)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACCC  >  NC_002947/4447744‑4448039
                                                                                                                                                     |                                                                                                                                                      
gTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCAt                                                                                                                                                       >  2:716115/1‑150 (MQ=255)
 ttGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATg                                                                                                                                                      <  2:128243/150‑1 (MQ=255)
   gATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAATAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGt                                                                                                                                                    <  2:2008464/150‑1 (MQ=255)
   gATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGt                                                                                                                                                    <  2:1428962/150‑1 (MQ=255)
    aTTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgc                                                                                                                                                  <  2:714087/150‑1 (MQ=255)
     ttGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTgc                                                                                                                                                  >  1:1920413/1‑150 (MQ=255)
      tGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgct                                                                                                                                                 <  1:1063656/150‑2 (MQ=255)
      tGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgct                                                                                                                                                 <  1:1165664/150‑2 (MQ=255)
           ccGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTCTACAGAGCACAAGTAAGTTTTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgctgcatc                                                                                                                                            <  1:2099839/150‑7 (MQ=255)
             gTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGt                                                                                                                                          >  2:1913478/1‑150 (MQ=255)
             gTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGt                                                                                                                                          >  2:1303191/1‑150 (MQ=255)
                 tccgatctCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAAc                                                                                                                                      <  1:571444/143‑1 (MQ=255)
                         cAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTgcgc                                                                                                                                                <  2:345062/132‑1 (MQ=255)
                         cAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTgcgc                                                                                                                                                >  1:345062/1‑132 (MQ=255)
                          aaTGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAg                                                                                                                             <  1:1183311/150‑1 (MQ=255)
                                        tAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGc                                                                                                               >  2:1671177/1‑150 (MQ=255)
                                          aCGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCg                                                                                                             <  2:2223658/150‑1 (MQ=255)
                                             tttGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTg                                                                                                         >  1:1794783/1‑150 (MQ=255)
                                                  aCAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgc                                                                                                                                                  >  1:851139/1‑105 (MQ=255)
                                                  aCAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgc                                                                                                                                                  <  2:851139/105‑1 (MQ=255)
                                                   cAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAAt                                                                                                   <  1:1303191/150‑1 (MQ=255)
                                                    aGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATg                                                                                                  <  2:1920413/150‑1 (MQ=255)
                                                         aCTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACaa                                                                                             <  1:920558/150‑1 (MQ=255)
                                                          cTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGa                                                                                                        <  2:1618527/138‑1 (MQ=255)
                                                          cTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGa                                                                                                        >  1:1618527/1‑138 (MQ=255)
                                                           ttGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACt                                                                                           >  1:1941814/1‑150 (MQ=255)
                                                           ttGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACt                                                                                           >  1:2155773/1‑150 (MQ=255)
                                                           ttGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACt                                                                                           >  1:504190/1‑150 (MQ=255)
                                                                 tcgctGTCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGATATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACACCGATGCCGGTGAGCAATGCACAACTG‑‑Accgt                                                                                   >  2:1189167/1‑148 (MQ=255)
                                                                   gctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgctgc                                                                                                                                               <  1:27996/91‑4 (MQ=38)
                                                                   gctGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgctgc                                                                                                                                               >  2:27996/1‑88 (MQ=38)
                                                                      ggCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGt                                                                              >  1:1099254/1‑150 (MQ=255)
                                                                        cAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTc                                                                            <  1:2014197/150‑1 (MQ=255)
                                                                           gATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑Accc                                                                                    >  1:2146170/1‑139 (MQ=255)
                                                                           gATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑Accc                                                                                    <  2:2146170/139‑1 (MQ=255)
                                                                            aTGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGAt                                                                                                                 <  2:1678215/112‑1 (MQ=255)
                                                                            aTGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGAt                                                                                                                 >  1:1678215/1‑112 (MQ=255)
                                                                              gCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCAc                                                                      >  1:1785660/1‑150 (MQ=255)
                                                                               cGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACg                                                                     <  1:1189167/150‑1 (MQ=255)
                                                                                      aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGaa                                                              <  1:716115/150‑1 (MQ=255)
                                                                                      aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGaa                                                              <  2:1785660/150‑1 (MQ=255)
                                                                                      aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGaa                                                              >  2:730038/1‑150 (MQ=255)
                                                                                       tCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGt                                                                                                                                                    <  1:774611/66‑1 (MQ=255)
                                                                                       tCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGt                                                                                                                                                    >  2:774611/1‑66 (MQ=255)
                                                                                       tCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAg                                                             <  2:2155773/150‑1 (MQ=255)
                                                                                       tCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAg                                                             >  1:2155658/1‑150 (MQ=255)
                                                                                       ccttcacaAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCAAAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAg                                                             >  1:1891849/5‑150 (MQ=255)
                                                                                         tGCACAAGTAGGGGGTACACACTGTTATTTTGCACCGAGAGGCAGCACCAGGCACTGCCATGGTGCGACCCTGTCAACGCAACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCCCGTACTCAAGcc                                                           >  2:486892/1‑149 (MQ=255)
                                                                                          gCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTc                                                                            <  1:1216989/132‑1 (MQ=255)
                                                                                          gCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTc                                                                            >  2:1216989/1‑132 (MQ=255)
                                                                                              aaGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCAATGGTTCGCCCCTGTCAACGCCACACAAGCGAACCCCGATGCCGCTGAACAATACACAACTG‑‑ACCCAGAAGTTCGAGCACGTACTGAAGCTGtctc                                                      >  2:1027196/1‑150 (MQ=255)
                                                                                                gTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGa                                                    >  2:2065233/1‑150 (MQ=255)
                                                                                                 tAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAg                                                   <  2:1099254/150‑1 (MQ=255)
                                                                                                 tAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAg                                                   <  1:730038/150‑1 (MQ=255)
                                                                                                  aGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGt                                                  <  2:1891849/150‑1 (MQ=255)
                                                                                                   gggTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTg                                                 >  2:291677/1‑150 (MQ=255)
                                                                                                    ggTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTgg                                                >  2:1062012/1‑150 (MQ=255)
                                                                                                     gtgtACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                               >  2:581682/1‑150 (MQ=255)
                                                                                                     gtgtACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGa                                               <  1:2065233/150‑1 (MQ=255)
                                                                                                             acTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTTCGCCCCTGTCAACGCCAACCAGGAGAACCCCGCTGACCGTGTGCAAGGCACAACGG‑‑AGACAGATCTTCGAGCACGTACTGAAGCTGTCTCGCGTGGACGAACCGc                                       >  2:642024/1‑150 (MQ=255)
                                                                                                               tattatTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCTCCCCTGTCAACGCACACCAGGAGAAAACCAAGATCGGAAGAGCA‑‑‑CACGTCTGAAATCCAGATGTTCGAGCACGTACTGAAGCTGTCTAGAATAGACGAAACGCag                                     >  1:551063/1‑150 (MQ=255)
                                                                                                               tattatTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAACGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCag                                     >  2:1277542/1‑150 (MQ=255)
                                                                                                                attattTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCaga                                    >  2:1108640/1‑150 (MQ=255)
                                                                                                                 ttattTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGca                                                                                                >  2:558510/1‑91 (MQ=255)
                                                                                                                 ttattTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCGGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGca                                                                                                <  1:558510/91‑1 (MQ=255)
                                                                                                                     tttGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTg                               <  2:2034483/150‑1 (MQ=255)
                                                                                                                       tGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGc                             <  2:504190/150‑1 (MQ=255)
                                                                                                                         cACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGccc                           >  2:1521556/1‑149 (MQ=255)
                                                                                                                          aCCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGt                          <  1:581682/150‑1 (MQ=255)
                                                                                                                            ttataTGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGc                        <  2:551063/146‑1 (MQ=255)
                                                                                                                            cggggtGCAGCACTACGCGCTGGCATGGTGCGCCCCTGTCAACGCAACCCAGGAGAACCCCGAGGCCGGGGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGc                        >  2:1296928/6‑150 (MQ=255)
                                                                                                                                  gcagcaCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGa                                                    <  1:1105775/116‑1 (MQ=255)
                                                                                                                                  gcagcaCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGa                                                    >  2:1105775/1‑116 (MQ=255)
                                                                                                                                   cagcaCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAagag                 <  1:291677/150‑1 (MQ=255)
                                                                                                                                    agcaCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGc                >  1:1703710/1‑150 (MQ=255)
                                                                                                                                     gcaCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGcc               >  1:574244/1‑150 (MQ=255)
                                                                                                                                      caCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGAGTGGCCGTGCTCAAGAGCCa              >  1:1929660/1‑150 (MQ=255)
                                                                                                                                       cccAGGCACTGCCATGGTGCGCCCCTGTCACCGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCac             <  2:1638057/149‑1 (MQ=255)
                                                                                                                                       aCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTTGCCGTGCTCAAGAGCCac             <  1:1062012/150‑1 (MQ=255)
                                                                                                                                       aCCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCac             >  1:1638057/1‑150 (MQ=255)
                                                                                                                                         cAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGt                                                         >  2:335025/1‑104 (MQ=255)
                                                                                                                                         cAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGt                                                         <  1:335025/104‑1 (MQ=255)
                                                                                                                                          aGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCactac          >  2:2128038/1‑150 (MQ=255)
                                                                                                                                           agacACTGCCATGGTGCCCCCCTGTCAAAGCCACCCAGAAGAACTTCGATTAAACAGAATTAAGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTAACGCTGTCTCGAATTGACAAACCTCCACTCGTGCCCGTGCTAAAAAGCCaatac          >  2:329235/4‑150 (MQ=255)
                                                                                                                                            ggCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCGGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCactact         >  2:1027048/1‑150 (MQ=255)
                                                                                                                                             gCACTGCCATGGTGCGCCCCTGTCACCGCTGCCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGCTGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTc        >  1:2245749/1‑150 (MQ=255)
                                                                                                                                                   ccATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTg                                                          <  1:1157362/94‑1 (MQ=255)
                                                                                                                                                   ccATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTg                                                          >  2:1157362/1‑94 (MQ=255)
                                                                                                                                                   ccATGGAGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGTGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGAccc  >  2:1543153/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTACTCCGACCC  >  NC_002947/4447744‑4448039

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