Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,937,181 T→C 5.0% intergenic (+13/‑113) rbsB → / → rbsK D‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor/ribokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,937,1810TC 5.0% 360.1 / 6.7 120intergenic (+13/‑113)rbsB/rbsKD‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor/ribokinase
Reads supporting (aligned to +/- strand):  ref base T (59/52);  new base C (3/3);  total (62/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 121 total aligned reads displayed.

GCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAG  >  NC_000913/3937032‑3937328
                                                                                                                                                     |                                                                                                                                                    
 cAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTg                                                                                                                                                     <  2:1644856/150‑1 (MQ=255)
     atAGCAGCGACTATCGCTCAGCTACCAGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATg                                                                                                                                                 <  1:62682/149‑1 (MQ=255)
        gcagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGAcc                                                                                                                                              <  1:1710307/150‑1 (MQ=255)
        gcagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGAcc                                                                                                                                              <  1:1557401/150‑1 (MQ=255)
         aagcGACTATCGCTCTTCTACCCTATCAGATTGGCGAGAAAGTCGTCGAAACCGCAGATAAATTGCTGAAATGAGAGAAAGTTCATGGTAAGTATCCGGTTGATCTGAATCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCt                                                                                                                                             <  1:236298/149‑1 (MQ=255)
          agcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTg                                                                                                                                            <  2:177966/150‑1 (MQ=255)
           gcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGa                                                                                                                                           >  1:758128/1‑150 (MQ=255)
              aCTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTCATACAACCTGAtta                                                                                                                                        <  2:709322/150‑3 (MQ=255)
               cTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGGCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGt                                                                                                                                       <  2:1286083/150‑1 (MQ=255)
                tATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                                                      <  2:408899/150‑1 (MQ=255)
                tATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                                                      <  2:1473143/150‑1 (MQ=255)
                tATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                                                      >  1:912825/1‑150 (MQ=255)
                tATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                                                      <  1:1261212/150‑1 (MQ=255)
                 aTCGCTCAGCTACGCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGa                                                                                                                                     <  1:625209/150‑1 (MQ=255)
                 aTCGCACAGGTACCCGAGCAGATTTGCGCGAAAGGCGGCGAAACCGCAGTTAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGTTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGa                                                                                                                                     <  2:820191/150‑1 (MQ=255)
                  tCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTCATACAACCTGAttaaaac                                                                                                                                    <  2:1208932/150‑7 (MQ=255)
                   cGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACa                                                                                                                                   <  1:794416/150‑1 (MQ=255)
                   cGCTCAGCTAACCGATCAGATTGGCGCGACAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACa                                                                                                                                   <  1:241159/150‑1 (MQ=255)
                    gcgcAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGATGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACAt                                                                                                                                  <  1:1334008/147‑1 (MQ=255)
                    gCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACAt                                                                                                                                  <  1:1444032/150‑1 (MQ=255)
                    gCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTATAATCAGGTTGTATGACCTGATGGTGACAt                                                                                                                                  <  1:1562823/150‑1 (MQ=255)
                     cTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATa                                                                                                                                 >  1:853454/1‑150 (MQ=255)
                     cTCAGCTACCCGATCAGAATTGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATa                                                                                                                                 <  1:799516/150‑1 (MQ=255)
                        aGCTACCCGATCTGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAAt                                                                                                                              <  1:1366259/150‑1 (MQ=255)
                        aGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAAt                                                                                                                              >  2:1594780/1‑150 (MQ=255)
                        aGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAAt                                                                                                                              >  1:671150/1‑150 (MQ=255)
                        aGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAAt                                                                                                                              >  1:1326783/1‑150 (MQ=255)
                          cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGAAAAAGTGCTGAAAGGCGAGAAAGTTCAGCCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGAAGCTGACATAAATAc                                                                                                                            >  2:841990/1‑150 (MQ=255)
                             cccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTATAATCAGGTTGTATGACCTGATGGTGACATAAATACGTc                                                                                                                         >  1:1646622/1‑150 (MQ=255)
                                   cAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGAc                                                                                                                   >  2:1332945/1‑150 (MQ=255)
                                    cgATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGATCTGAAAGGCGAGAAAGTTCAGCCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACa                                                                                                                  <  2:865958/149‑1 (MQ=255)
                                        tGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATg                                                                                                              <  1:1150342/150‑1 (MQ=255)
                                          gcgcGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGaa                                                                                                            >  1:1195060/1‑150 (MQ=255)
                                           cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAAc                                                                                                           >  1:1113162/1‑150 (MQ=255)
                                                aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATAAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTa                                                                                                      <  2:929080/150‑1 (MQ=255)
                                                 aGGCGTCGAAACCGCAGATAAAGGGCTGAAAGGCGAGAAAGTTCAGGCTATGAATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTGGTATGACCTGATGGTGACATAAATACGGCATCGACAGATGAACGTGTaa                                                                                                     <  2:976167/150‑1 (MQ=255)
                                                  ggCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCTGTTGATCTGAAACTGGGTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATCCGTCATCGACAGATGAACGTGTAat                                                                                                    <  1:612864/150‑1 (MQ=255)
                                                        gAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAg                                                                                              <  1:1469122/150‑1 (MQ=255)
                                                              gCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGc                                                                                        >  1:1691005/1‑150 (MQ=255)
                                                                 gATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAgg                                                                                    >  1:1432858/1‑150 (MQ=255)
                                                                  aTAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAggg                                                                                   <  2:470016/150‑1 (MQ=255)
                                                                   tAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGc                                                                                  >  1:1259063/1‑150 (MQ=255)
                                                                     aaGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCAc                                                                                >  2:839306/1‑150 (MQ=255)
                                                                      aGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcg                                                                               >  1:571584/1‑150 (MQ=255)
                                                                       gTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgc                                                                              >  1:1050009/1‑150 (MQ=255)
                                                                         gCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTATTCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACgcgc                                                                            >  2:1417076/1‑150 (MQ=255)
                                                                            gAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCAc                                                                         >  1:1285430/1‑150 (MQ=255)
                                                                            gAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCAc                                                                         <  2:1279120/150‑1 (MQ=255)
                                                                             aaaGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTATTCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCAcc                                                                        <  1:1417076/150‑1 (MQ=255)
                                                                             aaaGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCAcc                                                                        <  1:1673928/150‑1 (MQ=255)
                                                                             aaaGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCAcc                                                                        <  2:536923/150‑1 (MQ=255)
                                                                                  cGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAac                                                                   >  2:1112445/1‑150 (MQ=255)
                                                                                  cGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCAACCAAac                                                                   >  2:1668370/1‑150 (MQ=255)
                                                                                     gaAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGAAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACg                                                                >  2:921788/1‑150 (MQ=255)
                                                                                         gTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGTGCACGCGCCACCCTAACACGGTgg                                                            >  1:692853/1‑150 (MQ=255)
                                                                                          ttcttgCTAAGTTTGCTGTTGATCTGAATCTGGTTGTTAAGCTGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGGCACGCGCCACCCTAACACGGTgg                                                            <  2:1144848/145‑1 (MQ=255)
                                                                                           tCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcg                                                          >  2:70657/1‑150 (MQ=255)
                                                                                           tCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcg                                                          <  2:526703/150‑1 (MQ=255)
                                                                                           tCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcg                                                          >  1:1018437/1‑150 (MQ=255)
                                                                                              ggCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAt                                                       <  1:814639/150‑1 (MQ=255)
                                                                                              ggCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAt                                                       >  1:143945/1‑150 (MQ=255)
                                                                                              ggCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAt                                                       >  1:1560922/1‑150 (MQ=255)
                                                                                                cTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAttt                                                     <  2:1499847/150‑1 (MQ=255)
                                                                                                   aGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTAt                                                  >  1:427455/1‑150 (MQ=255)
                                                                                                     tATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATgg                                                >  1:606401/1‑150 (MQ=255)
                                                                                                      caccGCTTGATCCGAAACTGGTGGCTAAGCAGTAGTTATACTCAGGTGGTATGACATGAAGGTGACATAACTACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACCCGGTGGCGCAGTTTATGGa                                               <  2:306379/148‑1 (MQ=255)
                                                                                                      aTCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGa                                               >  1:166713/1‑150 (MQ=255)
                                                                                                       tCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGAc                                              <  2:853454/150‑1 (MQ=255)
                                                                                                         cGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACAt                                            >  2:1339595/1‑150 (MQ=255)
                                                                                                           gTTGATCTGAAACTGGTTGTTAAGCAGAAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGGCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATcc                                          <  2:1681466/150‑1 (MQ=255)
                                                                                                                 cTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCTCGCGCCCCCCTAACACGGTGGCGCATTTTATGGACAGCCCGAata                                    >  2:1055823/1‑150 (MQ=255)
                                                                                                                  tGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGTCACGCGCCACCCTAACACGGTGGCGCATGTTATGGACATACCGAataa                                   >  2:692389/1‑149 (MQ=255)
                                                                                                                  tGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAatat                                   >  1:687118/1‑150 (MQ=255)
                                                                                                                  tGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAatat                                   >  1:326147/1‑150 (MQ=255)
                                                                                                                   gAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATg                                  >  1:521759/1‑150 (MQ=255)
                                                                                                                   gAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATg                                  >  2:743669/1‑150 (MQ=255)
                                                                                                                    aaaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGc                                 >  1:283390/1‑150 (MQ=255)
                                                                                                                      aCTGGTTGTTAAGCAGTAGTTTTAATCAGTTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCTGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGATTATGCaa                               <  1:1508019/150‑1 (MQ=255)
                                                                                                                         ggTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAAc                            <  2:593628/150‑1 (MQ=255)
                                                                                                                          gttgttAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACg                           <  2:1725411/150‑1 (MQ=255)
                                                                                                                               tAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACgcaggc                      >  1:86990/1‑150 (MQ=255)
                                                                                                                                aaGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCAGCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACgcaggca                     >  1:95506/1‑150 (MQ=255)
                                                                                                                                 aGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACgcaggcag                    >  2:1446696/1‑150 (MQ=255)
                                                                                                                                 aGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACgcaggcag                    >  2:163219/1‑150 (MQ=255)
                                                                                                                                    agtagtTTTAATCAGGTGGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCt                 <  2:163058/150‑1 (MQ=255)
                                                                                                                                       agtTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGTCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCgt              <  1:680172/150‑1 (MQ=255)
                                                                                                                                       agtTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCgg              >  1:1568933/1‑149 (MQ=255)
                                                                                                                                        gtcaccaTCAGGTCATACAACCTGATGGTGCCATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCCCCCTAACACGGTGGCGCATTTTATGGACATCCAGAAAATGCAAAACGCAGGCAGACTCgtt             >  1:965293/7‑150 (MQ=255)
                                                                                                                                        gtcaccaTCAGGTCATACAACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAGTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCgtt             >  1:201242/7‑150 (MQ=255)
                                                                                                                                         ttttAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCgttg            <  2:912825/150‑1 (MQ=255)
                                                                                                                                         ttttAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCACATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCgttg            >  1:257367/1‑150 (MQ=255)
                                                                                                                                            tAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCTGGCCGCCTCgttgtta         >  1:1398855/1‑149 (MQ=255)
                                                                                                                                            tAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTc         >  2:739884/1‑150 (MQ=255)
                                                                                                                                            ccaTCAGGTCATACAACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTc         >  1:971347/3‑150 (MQ=255)
                                                                                                                                              aTCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACAAGGTGGCGCATTTTATGGACATCCCTAATATGCAAAAAGCATGCAGCCTAGTTGTTCtt       >  1:578827/1‑150 (MQ=255)
                                                                                                                                                cAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCCCCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTgg     >  2:1289581/1‑150 (MQ=255)
                                                                                                                                                 aGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGCTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGAAGCATCGTTGTTCTTGgc    >  1:103907/1‑150 (MQ=255)
                                                                                                                                                  ggTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGgca   >  1:1051294/1‑150 (MQ=255)
                                                                                                                                                   gTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAACCGCAGGCAGCCTCGTTGTTCTTGgcag  >  2:260697/1‑150 (MQ=255)
                                                                                                                                                   gTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGgcag  <  1:986573/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAG  >  NC_000913/3937032‑3937328

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