Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,294,120 C→T 5.4% intergenic (+15/+25) pdxH → / ← PP_1130 pyridoxine/pyridoxamine 5'‑phosphate oxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,294,1200CT 5.4% 279.4 / 4.5 93intergenic (+15/+25)pdxH/PP_1130pyridoxine/pyridoxamine 5'‑phosphate oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (46/42);  new base T (3/2);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

GAGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCT  >  NC_002947/1293974‑1294268
                                                                                                                                                  |                                                                                                                                                    
gAGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAg                                                                                                                                                   >  2:58598/1‑150 (MQ=255)
 aGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGTGTGCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGc                                                                                                                                                  <  1:201719/150‑1 (MQ=255)
 aGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGc                                                                                                                                                  <  2:2855928/150‑1 (MQ=255)
 aGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGc                                                                                                                                                  <  2:2750119/150‑1 (MQ=255)
    aTTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGc                                                                                                                                               >  1:2557366/1‑150 (MQ=255)
     gggggggtgATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGCGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTCGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGcc                                                                                                                                              <  2:184877/148‑1 (MQ=255)
      tGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGccc                                                                                                                                             <  1:1763928/150‑1 (MQ=255)
           gggATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGcc                                                                                                                                        >  1:524971/1‑150 (MQ=255)
              aTACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCAtt                                                                                                                                     <  2:949028/150‑1 (MQ=255)
                  cGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGacc                                                                                                                                 >  1:1077803/1‑148 (MQ=255)
                  cGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATc                                                                                                                                 >  1:1156723/1‑150 (MQ=255)
                   gTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGc                                                                                                                                                  >  2:959300/1‑132 (MQ=255)
                   gTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGc                                                                                                                                                  <  1:959300/132‑1 (MQ=255)
                    tttGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTTGCTGCGAGATCGCCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCaa                                                                                                                               >  1:695031/1‑150 (MQ=255)
                    tttGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCaa                                                                                                                               <  2:329791/150‑1 (MQ=255)
                      tactgCCCGAGCGCCTCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCAGCCAAGCGGCCCCAGCCATTGATCAAGc                                                                                                                             <  1:2282600/147‑1 (MQ=255)
                           cccGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTAccc                                                                                                                        <  1:1205419/150‑1 (MQ=255)
                             cGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCg                                                                                                                      >  1:1867572/1‑150 (MQ=255)
                                  gcATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGc                                                                                                                                               >  2:464964/1‑120 (MQ=255)
                                  gcATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGc                                                                                                                                               <  1:464964/120‑1 (MQ=255)
                                    aTCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAAcc                                                                                                               >  1:2662205/1‑150 (MQ=255)
                                     tCGAGTTCTGGCAGGGGCGGGCAAGTCTTCTGCAGGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGAGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAAccc                                                                                                              <  1:1199595/150‑2 (MQ=255)
                                       gAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGCCCCAgt                                                                                                                                         >  1:2450036/1‑120 (MQ=255)
                                       gAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCTTGGCGGCCCCAgt                                                                                                                                         <  2:2450036/121‑2 (MQ=255)
                                        aGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTgcc                                                                                                           >  1:1894026/1‑150 (MQ=255)
                                        aGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCACTGGCTGGGAGAGCGGCTGGCGCCCTGAGCAACTGGGCCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTgcc                                                                                                           >  2:1036753/1‑150 (MQ=255)
                                             tGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTAcc                                                                                                                         <  2:545313/131‑1 (MQ=255)
                                             tGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTAcc                                                                                                                         >  1:545313/1‑131 (MQ=255)
                                             tGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTc                                                                                                      >  2:1547671/1‑150 (MQ=255)
                                               gCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAATGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAACGGCCCCAGCCATTGATCAAGCTACCCCGTATAACCTGCCGCCTctc                                                                                                    >  2:378442/1‑148 (MQ=255)
                                                   gggcgggcAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTc                                                                                                >  1:2623489/1‑150 (MQ=255)
                                                          cAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAaccac                                                                                         >  1:1889585/1‑150 (MQ=255)
                                                           aaGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCaa                                                                                                                               >  1:2964402/1‑111 (MQ=255)
                                                           aaGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCaa                                                                                                                               <  2:2964402/111‑1 (MQ=255)
                                                            aGTCGTCTGCATGATAGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCAccc                                                                                       >  1:69963/1‑150 (MQ=255)
                                                                 tcTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCccacca                                                                                  <  2:2557366/150‑1 (MQ=255)
                                                                 tcTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCa                                                                                  <  1:50637/150‑1 (MQ=255)
                                                                    gCATGCTCGTTTGCACTCTTCGCAGGTGCACGGGGAGTTGCTGCTAGTTCTGCCGTTGCTCTTGGGAACTTGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGtt                                                                               >  2:32220/1‑150 (MQ=255)
                                                                     cATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTg                                                                              >  1:1686822/1‑150 (MQ=255)
                                                                     cATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTg                                                                              <  2:1686822/150‑1 (MQ=255)
                                                                     cATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGAGGCCGCCAAGCGGCCCCAGCCAATGTTCCAGCTACCCCAGATAACCTGCGGCCTCGCGCTCCACCCACCCCGCGAGTTg                                                                              >  1:2934985/1‑150 (MQ=255)
                                                                     cATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCCAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGttt                                                                              >  2:2227028/1‑149 (MQ=255)
                                                                      aTGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGc                                                                             >  1:2267163/1‑150 (MQ=255)
                                                                       tGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGcc                                                                            >  2:2065373/1‑149 (MQ=255)
                                                                            gTTTTAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGc                                                                       <  2:1867572/150‑1 (MQ=255)
                                                                             tttGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGcc                                                                      >  1:2896735/1‑150 (MQ=255)
                                                                             tttGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGcc                                                                      >  2:789566/1‑150 (MQ=255)
                                                                                 aaCTATCGGCTGGTTGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATaa                                                                                                                 >  2:119437/1‑103 (MQ=255)
                                                                                 aaCTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATaa                                                                                                                 <  1:119437/103‑1 (MQ=255)
                                                                                 aaCTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCAc                                                                  >  2:2627243/1‑150 (MQ=255)
                                                                                      tCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTtct                                                             >  1:1359261/1‑150 (MQ=255)
                                                                                         gCTGGTGGCCGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCt                                                          <  1:2310920/150‑1 (MQ=255)
                                                                                         gCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGCTCTCGCTACTCCGGATAACCGTTCGCCGCGCTCTGTAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCt                                                          <  1:58598/150‑1 (MQ=255)
                                                                                         gCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGACCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCt                                                          <  2:1077803/150‑1 (MQ=255)
                                                                                           tggtggACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGc                                                        >  1:1139135/1‑150 (MQ=255)
                                                                                            ggtggACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGcc                                                       <  2:1139135/150‑1 (MQ=255)
                                                                                                 aCGGGCAGTGGCTGCGAGAGCGGCTGGCGCTCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGcc                                                  >  1:2387190/1‑150 (MQ=255)
                                                                                                 aCGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGc                                                                                                                             <  2:605134/75‑1 (MQ=255)
                                                                                                 aCGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGc                                                                                                                             >  1:605134/1‑75 (MQ=255)
                                                                                                    ggCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTgc                                               <  1:378442/150‑1 (MQ=255)
                                                                                                      cAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTgcgc                                             <  2:1156723/150‑1 (MQ=255)
                                                                                                         tGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCaa                                          <  2:2662205/150‑1 (MQ=255)
                                                                                                          ggCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAAc                                         <  1:2627243/150‑1 (MQ=255)
                                                                                                             tGCGAGAGCGGCTGGCGCGCTGAGCAACTGGGGCCGCCAAGCCGTGCCAGCCATTGAAGAAGATATCCCCGGTAACATACCGCCTCGCGCTCCAACCACTCCGGAAATAGCGTAAGTCACAGTGTCTCTTTCGGCGCCTGCGCCAACAtt                                      >  2:1306751/1‑150 (MQ=255)
                                                                                                                gagaGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTcc                                   >  1:672435/1‑150 (MQ=255)
                                                                                                                    gCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCa                               <  2:524971/150‑1 (MQ=255)
                                                                                                                             gcCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAcgcgc                      >  2:279066/1‑150 (MQ=255)
                                                                                                                             gcCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCAACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAcgcgc                      <  2:2896735/150‑1 (MQ=255)
                                                                                                                               ccTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAac                                                                                            >  1:1224938/1‑78 (MQ=255)
                                                                                                                               ccTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAac                                                                                            <  2:1224938/78‑1 (MQ=255)
                                                                                                                               ccTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTa                                                                          >  2:2791116/1‑96 (MQ=255)
                                                                                                                               ccTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTa                                                                          <  1:2791116/96‑1 (MQ=255)
                                                                                                                                 tGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAaccac                                                                                         >  2:2774837/1‑79 (MQ=255)
                                                                                                                                 tGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAaccac                                                                                         <  1:2774837/79‑1 (MQ=255)
                                                                                                                                  gAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTcc                                                           <  1:1810072/108‑1 (MQ=255)
                                                                                                                                  gAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTcc                                                           >  2:1810072/1‑108 (MQ=255)
                                                                                                                                    gCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTc               >  2:537441/1‑150 (MQ=255)
                                                                                                                                    gCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTc               <  1:789566/150‑1 (MQ=255)
                                                                                                                                     cAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTcc                                                           >  1:477402/1‑105 (MQ=255)
                                                                                                                                     cAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTcc                                                           <  2:477402/105‑1 (MQ=255)
                                                                                                                                     cAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCAAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCa              <  1:1689805/150‑1 (MQ=255)
                                                                                                                                      aaCTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAg             <  2:69963/150‑1 (MQ=255)
                                                                                                                                       gcTGGGGCCGCTTGGCGGCCCCAGTCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAg             >  1:3051108/2‑149 (MQ=255)
                                                                                                                                       gcTGGGGCCGCTTGGCGGCCCCAGTCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAg             <  2:3051108/148‑1 (MQ=255)
                                                                                                                                       aCTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAgg            >  2:1237345/1‑150 (MQ=255)
                                                                                                                                        cTGGGGCCGCTTGGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGa           >  1:783943/1‑150 (MQ=255)
                                                                                                                                         tGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGTGCGTTTTTCAGGAt          >  2:1604470/1‑150 (MQ=255)
                                                                                                                                          ggggCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGTCGTGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATc         <  2:2134791/150‑1 (MQ=255)
                                                                                                                                           gggCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAc                          <  2:2940378/132‑1 (MQ=255)
                                                                                                                                           gggCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAc                          <  2:2940234/132‑1 (MQ=255)
                                                                                                                                           gggCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAc                          >  1:2940378/1‑132 (MQ=255)
                                                                                                                                           gggCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAAc                          >  1:2940234/1‑132 (MQ=255)
                                                                                                                                            ggCCGCCAAGCGGCCCCAGCCATTGATCAAGGTACCCGGCATAACCTGCCGCCTCGCGCTCCAACGGCCTGGCCAGTGGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCgg       <  1:2065373/150‑1 (MQ=255)
                                                                                                                                               cgccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGccc    <  2:2623489/150‑1 (MQ=255)
                                                                                                                                                gccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGcccc   >  2:1051772/1‑150 (MQ=255)
                                                                                                                                                gccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGcccc   >  2:779141/1‑150 (MQ=255)
                                                                                                                                                 ccAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCt  <  2:1889585/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GAGCATTGGGGGGGATACCGTTTGTTGCCCGAGCGCATCGAGTTCTGGCAGGGGCGGGCAAGTCGTCTGCATGATCGTTTGAACTATCGGCTGGTGGACGGGCAGTGGCTGCGAGAGCGGCTGGCGCCCTGAGCAACTGGGGCCGCCAAGCGGCCCCAGCCATTGATCAAGCTACCCCGGATAACCTGCCGCCTCGCGCTCCAACCACCCCGCCAGTTGCTTACGCCCCACCTTCTCCTGCCGCGCCTGCGCCAACCTTTCCAGCATGAACGCGCGTTTTTCAGGATCGGCCCCT  >  NC_002947/1293974‑1294268

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