Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,447,882:1 +T 11.5% intergenic (‑274/‑48) 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,8821.T11.5% 297.5 / 37.3 87intergenic (‑274/‑48)nicC/nicX6‑hydroxynicotinate 3‑monooxygenase/2,5‑dihydroxypyridine 5,6‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base . (41/36);  new base T (4/6);  total (45/42)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.79e-01

GAAATCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTAC  >  NC_002947/4447736‑4448031
                                                                                                                                                   |                                                                                                                                                        
gaaaTCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGc                                                                                                                                                       <  1:689711/150‑1 (MQ=255)
gaaaTCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGc                                                                                                                                                       >  1:2034483/1‑150 (MQ=255)
  aaTCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCAc                                                                                                                                                     <  1:752519/150‑1 (MQ=255)
        gTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCAt                                                                                                                                               >  2:716115/1‑150 (MQ=255)
         ttGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATg                                                                                                                                              <  2:128243/150‑1 (MQ=255)
          tGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGcctgg                                                                                                                                              >  2:1063656/1‑147 (MQ=255)
          tGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGcctggt                                                                                                                                             >  1:714087/1‑146 (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:1303191/1‑150 (MQ=255)
                     gTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGt                                                                                                                                  >  2:1913478/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                                                                                                                                          <  2:851139/105‑1 (MQ=255)
                                                          aCAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCATGGCAGTGCCTGGTgc                                                                                                                                          >  1:851139/1‑105 (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                                                      <  2:1785660/150‑1 (MQ=255)
                                                                                              aTCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGaa                                                      <  1:716115/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                                                                    >  2:1216989/1‑132 (MQ=255)
                                                                                                  gCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTc                                                                    <  1:1216989/132‑1 (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)
                                                                                                                                                   |                                                                                                                                                        
GAAATCCGGTTGATTGTTACCGTAAATACTGGTCAATGTGTGGATAGGTAACGTTTGAACAGATCACTTGCGCTCGCTGGCAGGATGCGTTTTCATCTGCACAAGTAGGGTGTACACACTATTATTTTGCACCGAGATGCAGCACCAGGCACTGCCATGGTGCGCCCCTGTCAACGCCACCCAGGAGAACCCCGATGCCGGTGAGCAATGCACAACTG‑‑ACCCAGATGTTCGAGCACGTACTGAAGCTGTCTCGAGTGGACGAAACGCAGAGCGTGGCCGTGCTCAAGAGCCACTAC  >  NC_002947/4447736‑4448031

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: