Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,049,346 G→T 5.8% intergenic (‑32/+28) PP_1822 ← / ← folEA‑I hypothetical protein/GTP cyclohydrolase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,049,3460GT 5.8% 301.8 / 4.6 96intergenic (‑32/+28)PP_1822/folEA‑Ihypothetical protein/GTP cyclohydrolase I
Reads supporting (aligned to +/- strand):  ref base G (42/47);  new base T (4/2);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

GCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACATG  >  NC_002947/2049199‑2049494
                                                                                                                                                   |                                                                                                                                                     
gCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGaa                                                                                                                                                     <  1:633232/150‑1 (MQ=255)
 cTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGCAGTTGAGCTGGCCGAGGCTTGCCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAg                                                                                                                                                    <  1:886423/150‑1 (MQ=255)
  ttGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAATTTCTGAAATGCAAAAGCCGACCCTGAAgg                                                                                                                                                   <  1:665863/150‑1 (MQ=255)
  ttGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAgg                                                                                                                                                   <  1:589269/150‑1 (MQ=255)
  ttGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAgg                                                                                                                                                   <  1:915899/150‑1 (MQ=255)
     cTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTc                                                                                                                                                <  2:1081208/150‑1 (MQ=255)
        gATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGc                                                                                                                                             >  1:105697/1‑150 (MQ=255)
            cGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCgg                                                                                                                                         <  1:2058401/150‑1 (MQ=255)
                 cTGGGCGGCGGCGTCGCGTTTGCGCTTGGCAGGGCGGCTGTTCCACTCACCGAAGACGATGAGCTGGCCGAGGGCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCa                                                                                                                                    <  1:1115094/150‑1 (MQ=255)
                  tGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGaa                                                                                                                                                     >  1:1959728/1‑132 (MQ=255)
                  tGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGaa                                                                                                                                                     <  2:1959728/132‑1 (MQ=255)
                  tGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAgg                                                                                                                                                   <  2:1267021/134‑1 (MQ=255)
                  tGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAgg                                                                                                                                                   >  1:1267021/1‑134 (MQ=255)
                  tGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAg                                                                                                                                   <  2:800424/150‑1 (MQ=255)
                       ggcggcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCgg                                                                                                                                         >  1:1513288/1‑139 (MQ=255)
                       ggcggcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCgg                                                                                                                                         <  2:1513288/139‑1 (MQ=255)
                         cggcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCaa                                                                                                                           <  1:1710547/150‑1 (MQ=255)
                           gcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCaa                                                                                                                           >  1:204580/1‑148 (MQ=255)
                           gcgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCaa                                                                                                                           <  2:204580/148‑1 (MQ=255)
                            cgTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCACCCCGCAGGGCTGTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCa                                                                                                                        <  2:1381415/150‑1 (MQ=255)
                             gTCGCGTTTGCGCTTGGCCGGGCGGCTTATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCAc                                                                                                                       >  1:2277526/1‑150 (MQ=255)
                              tCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACt                                                                                                                      <  2:644454/150‑1 (MQ=255)
                                  gTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCCGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTCAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGAt                                                                                                                  <  1:1044288/150‑1 (MQ=255)
                                         cTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGATCGGCCCGGGGCCAGCAACCAATCACTTGATCAGGCTg                                                                                                           >  2:1145307/1‑150 (MQ=255)
                                          ttGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAgg                                                                                                                                                   <  2:41182/110‑1 (MQ=255)
                                          ttGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAgg                                                                                                                                                   >  1:41182/1‑110 (MQ=255)
                                          ttGGCCGGGCGGCTGATCCAGGCACCTAAGGGTATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGTTCATGGGCAACTCCTGAAAGGCAAAAGCCGACCCTTAAGGGTCTGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGa                                                                                                          <  2:989046/150‑1 (MQ=255)
                                             gCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGTCCAGCAACCAATCACtt                                                                                                                     >  2:472172/1‑136 (MQ=255)
                                             gCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGTCCAGCAACCAATCACtt                                                                                                                     <  1:472172/136‑1 (MQ=255)
                                                   cggcTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGc                                                                                                 <  2:316034/150‑1 (MQ=255)
                                                      cTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGc                                                                                              <  2:105697/150‑1 (MQ=255)
                                                      cTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGc                                                                                              <  1:1059366/150‑1 (MQ=255)
                                                         aTCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCAc                                                                                                                       >  2:915718/1‑122 (MQ=255)
                                                         aTCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCAc                                                                                                                       <  1:915718/122‑1 (MQ=255)
                                                         aTCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCggg                                                                                           >  1:1470727/1‑150 (MQ=255)
                                                            cAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTgg                                                                                        >  1:1030282/1‑150 (MQ=255)
                                                                   cGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGtt                                                                                 >  2:147671/1‑150 (MQ=255)
                                                                          gATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTcgc                                                                            >  2:452312/1‑148 (MQ=255)
                                                                          gATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTcgc                                                                            <  1:452312/148‑1 (MQ=255)
                                                                               gCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTc                                                                     >  1:2155711/1‑150 (MQ=255)
                                                                               gCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGCCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTc                                                                     <  1:221515/150‑1 (MQ=255)
                                                                                    ccGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCa                                                                >  1:969332/1‑150 (MQ=255)
                                                                                       aGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCa                                                             <  1:667721/150‑1 (MQ=255)
                                                                                         gCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTcgc                                                                            <  1:2412231/133‑1 (MQ=255)
                                                                                         gCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGACAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTcgc                                                                            >  2:2412231/1‑133 (MQ=255)
                                                                                           cAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAc                                                         >  2:902694/1‑150 (MQ=255)
                                                                                            aaCCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAcc                                                        >  2:65673/1‑150 (MQ=255)
                                                                                                  cAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCAcc                                                                  <  1:2065741/134‑1 (MQ=255)
                                                                                                  cAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCAcc                                                                  >  2:2065741/1‑134 (MQ=255)
                                                                                                     ggCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATc                                               >  2:680399/1‑150 (MQ=255)
                                                                                                        tttGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAATCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGa                                                      <  2:1233892/140‑1 (MQ=255)
                                                                                                        tttGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAATCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGa                                                      >  1:1233892/1‑140 (MQ=255)
                                                                                                        tttGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCg                                                                                                  >  1:1077594/1‑96 (MQ=255)
                                                                                                        tttGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCg                                                                                                  <  2:1077594/96‑1 (MQ=255)
                                                                                                         ttGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCTCGGTGAtgctgg                                           <  2:2277526/150‑6 (MQ=255)
                                                                                                         ttGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGACGAGCCTCGGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCACCACCGAGGTGATCATgg                                           >  1:413827/1‑150 (MQ=255)
                                                                                                           gACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTc                                         <  1:65673/150‑1 (MQ=255)
                                                                                                              gATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCgg                                                                                            >  1:2384220/1‑96 (MQ=255)
                                                                                                              gATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCgg                                                                                            <  2:2384220/96‑1 (MQ=255)
                                                                                                              gATCATGCGCAACTCCTTAAATGTAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAGATAACATATCAGGCTGAGGAACTCGCTGCGGGTGTCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGaa                                      <  1:147671/150‑1 (MQ=255)
                                                                                                                tCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGc                                                                                                                                      >  1:543176/1‑53 (MQ=255)
                                                                                                                tCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGc                                                                                                                                      <  2:543176/53‑1 (MQ=255)
                                                                                                                   tGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTg                                 <  2:1002463/150‑1 (MQ=255)
                                                                                                                      gCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCtt                              >  2:280464/1‑150 (MQ=255)
                                                                                                                         aCTCCTGAAATGCAAAGGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                           <  2:1221945/150‑1 (MQ=255)
                                                                                                                         aCTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                           <  1:649598/150‑1 (MQ=255)
                                                                                                                         aCTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACGCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                           <  2:940786/150‑1 (MQ=255)
                                                                                                                           tCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                           >  2:2099705/1‑148 (MQ=255)
                                                                                                                           tCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtctc                           <  1:2099705/148‑1 (MQ=255)
                                                                                                                           tCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGa                         >  1:1560231/1‑150 (MQ=255)
                                                                                                                           tCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGa                         >  1:1582196/1‑150 (MQ=255)
                                                                                                                            ccTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGACCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGt                                          >  2:2293854/1‑132 (MQ=255)
                                                                                                                            ccTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGACCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGt                                          <  1:2293854/132‑1 (MQ=255)
                                                                                                                              tGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGacac                      >  1:1209938/1‑150 (MQ=255)
                                                                                                                              tGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGacac                      >  1:164538/1‑150 (MQ=255)
                                                                                                                                 aaTGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACcgc                   >  1:981883/1‑150 (MQ=255)
                                                                                                                                   tGCAAAAGAAGGCCCTGGAAGGTCGGAACTACGCCGACAACCAAGCACTTGATCAGGCTGACGCACAGGCCACCGTGGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCgcgc                 <  1:2253611/150‑1 (MQ=255)
                                                                                                                                     cccgaGCCGACCCTTCAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtct                            <  2:893538/133‑1 (MQ=255)
                                                                                                                                     cccgaGCCGACCCTTCAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTtct                            >  1:893538/5‑137 (MQ=255)
                                                                                                                                      ccgaGCCGACCCTTCAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatc              >  1:1753317/4‑150 (MQ=255)
                                                                                                                                       cgaGCCGACCCTTCAGGGTCGGCTTGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatca             >  2:2396963/3‑150 (MQ=255)
                                                                                                                                           ccGACCCTTCAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatcat         >  1:2191824/1‑150 (MQ=255)
                                                                                                                                           ccGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCTCTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatcat         <  2:1102621/150‑1 (MQ=255)
                                                                                                                                           ccGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAcc                                                        >  1:1166700/1‑103 (MQ=255)
                                                                                                                                           ccGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCAcc                                                        <  2:1166700/103‑1 (MQ=255)
                                                                                                                                           ccGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGcatcatcat         >  1:1160315/1‑150 (MQ=255)
                                                                                                                                              aCCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGca      <  2:1470727/150‑1 (MQ=255)
                                                                                                                                                ccTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAAc                                                                       <  2:99459/83‑1 (MQ=255)
                                                                                                                                                ccTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAAc                                                                       >  1:99459/1‑83 (MQ=255)
                                                                                                                                                ccTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGac                        >  1:2374430/1‑130 (MQ=255)
                                                                                                                                                ccTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGac                        <  2:2374430/130‑1 (MQ=255)
                                                                                                                                                ccTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGcaca    >  1:2097137/1‑150 (MQ=255)
                                                                                                                                                 cTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACAt   >  2:352229/1‑150 (MQ=255)
                                                                                                                                                  tGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAg                                                                                                         <  1:27947/47‑1 (MQ=255)
                                                                                                                                                  tGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAg                                                                                                         >  2:27947/1‑47 (MQ=255)
                                                                                                                                                  tGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACAtg  >  1:1132815/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTTGCTCGATGCGGGCCTGGGCGGCGGCGTCGCGTTTGCGCTTGGCCGGGCGGCTGATCCAGTCACCGAAGACGATGAGCTGGCCGAGGCCAACCCGCAGGGCTTTGACGATCATGGGCAACTCCTGAAATGCAAAAGCCGACCCTGAAGGGTCGGCTCGGGGCCAGCAACCAATCACTTGATCAGGCTGAGGAACTCGCTGCGGGTGGCGGCGTTGTCGCGGAACTCACCCAGCATCACCGAGGTGATCATGGTCGAATTCTGCTTCTCGACACCGCGCATCATCATGCACATG  >  NC_002947/2049199‑2049494

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