Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,937,185 T→C 5.0% intergenic (+17/‑109) 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,1850TC 5.0% 367.0 / 6.8 119intergenic (+17/‑109)rbsB/rbsKD‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor/ribokinase
Reads supporting (aligned to +/- strand):  ref base T (59/54);  new base C (3/3);  total (62/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.98e-01

Only 100 of 119 total aligned reads displayed.

CTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAG  >  NC_000913/3937037‑3937328
                                                                                                                                                    |                                                                                                                                                
atAGCAGCGACTATCGCTCAGCTACCAGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATg                                                                                                                                                 <  1:62682/149‑1 (MQ=255)
   gcagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGAcc                                                                                                                                              <  1:1557401/150‑1 (MQ=255)
   gcagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGAcc                                                                                                                                              <  1:1710307/150‑1 (MQ=255)
   gcagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGAcc                                                                                                                                              <  1:600262/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)
         aCTCTCGCTCAGCTACCCGATCAGATTGGCGCGAAACGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCTAGAAAGTTCATGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGGTTTAATCAGGTTGTATGACCTGATgg                                                                                                                                        <  1:1292357/150‑1 (MQ=255)
         aCTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTCATACAACCTGAtta                                                                                                                                        <  2:709322/150‑3 (MQ=255)
           tATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                                                      <  1:1261212/150‑1 (MQ=255)
           tATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                                                      <  2:168437/150‑1 (MQ=255)
           tATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                                                      <  2:408899/150‑1 (MQ=255)
           tATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                                                      >  1:912825/1‑150 (MQ=255)
           tATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTg                                                                                                                                      <  2:1473143/150‑1 (MQ=255)
            aTCGCTCAGCTACGCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGa                                                                                                                                     <  1:625209/150‑1 (MQ=255)
            aTCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGTTGa                                                                                                                                     >  1:1122707/1‑150 (MQ=255)
            aTCGCACAGGTACCCGAGCAGATTTGCGCGAAAGGCGGCGAAACCGCAGTTAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGTTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGa                                                                                                                                     <  2:820191/150‑1 (MQ=255)
             tCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTCATACAACCTGAttaaaac                                                                                                                                    <  2:1208932/150‑7 (MQ=255)
              cGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACa                                                                                                                                   <  1:794416/150‑1 (MQ=255)
              cGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACa                                                                                                                                   >  1:30740/1‑150 (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)
                     cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGAAAAAGTGCTGAAAGGCGAGAAAGTTCAGCCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGAAGCTGACATAAATAc                                                                                                                            >  2:841990/1‑150 (MQ=255)
                        cccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTATAATCAGGTTGTATGACCTGATGGTGACATAAATACGTc                                                                                                                         >  1:1646622/1‑150 (MQ=255)
                              cAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGAc                                                                                                                   >  2:1332945/1‑150 (MQ=255)
                                gATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAg                                                                                                                 >  1:1725411/1‑150 (MQ=255)
                                     gcgcGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGaa                                                                                                            >  1:1195060/1‑150 (MQ=255)
                                      cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAAc                                                                                                           >  1:1113162/1‑150 (MQ=255)
                                         gAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                                                        <  2:1374515/150‑1 (MQ=255)
                                           aaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTa                                                                                                      <  2:798205/150‑1 (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)
                                                             aTAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAggg                                                                                   <  2:470016/150‑1 (MQ=255)
                                                              tAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGc                                                                                  >  1:1259063/1‑150 (MQ=255)
                                                                aaGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCAc                                                                                >  2:839306/1‑150 (MQ=255)
                                                                  gTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgc                                                                              >  1:1050009/1‑150 (MQ=255)
                                                                  gTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgc                                                                              >  1:1430572/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                                                          >  1:1018437/1‑150 (MQ=255)
                                                                                      tCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcg                                                          >  2:70657/1‑150 (MQ=255)
                                                                                      tCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcg                                                          <  2:526703/150‑1 (MQ=255)
                                                                                       cAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcgc                                                         <  2:883261/150‑1 (MQ=255)
                                                                                         ggCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAt                                                       <  1:814639/150‑1 (MQ=255)
                                                                                         ggCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAt                                                       >  1:143945/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)
                                                                                                tATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATgg                                                <  2:1626859/150‑1 (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                                   >  2:644383/1‑150 (MQ=255)
                                                                                                             tGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAatat                                   >  1:326147/1‑150 (MQ=255)
                                                                                                              gAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATg                                  >  1:521759/1‑150 (MQ=255)
                                                                                                              gAAACTGGATGTTAACCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATCAACGTCTAATATAAAGAAAAGCATAGCACGCGCCACCCAAACACGGTGGCGCATTTAATGGACATCCCGAatata                                  >  1:15529/1‑149 (MQ=255)
                                                                                                               aaaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGc                                 <  2:1252404/150‑1 (MQ=255)
                                                                                                               aaaCTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGc                                 >  1:283390/1‑150 (MQ=255)
                                                                                                                 aCTGGTTGTTAAGCAGTAGTTTTAATCAGTTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCTGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGATTATGCaa                               <  1:1508019/150‑1 (MQ=255)
                                                                                                                    ggTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAAc                            <  2:30740/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)
                                                                                                                               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)
                                                                                                                                       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)
                                                                                                                                              gTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGgcag  <  1:986573/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAATATGCAAAACGCAGGCAGCCTCGTTGTTCTTGGCAG  >  NC_000913/3937037‑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: