Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,091,793 C→T G168D (GGC→GAC)  rhaD ← rhamnulose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,091,7930CT99.6% 2043.6 / NA 530G168D (GGC→GAC) rhaDrhamnulose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (278/249);  minor base G (1/1);  total (279/251)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 571 total aligned reads displayed.

CATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATC  >  NC_000913/4091693‑4091890
                                                                                                    |                                                                                                 
  tctcTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCt                                                                                                 >  2:5222138/1‑101 (MQ=255)
    tctTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTc                                                                                               >  1:6608749/1‑101 (MQ=255)
    tctTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTc                                                                                               >  1:14666861/1‑101 (MQ=255)
      ttGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTccc                                                                                             <  2:617213/101‑1 (MQ=255)
           cGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                                                        >  1:10249304/1‑101 (MQ=255)
            nnncGCATGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                                       <  2:13278431/98‑1 (MQ=255)
            ggTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                                       >  1:9687890/1‑101 (MQ=255)
            ggTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                                       >  1:8029310/1‑101 (MQ=255)
            ggTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                                       >  1:4750166/1‑101 (MQ=255)
               cGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGc                                                                                    >  2:6585715/1‑101 (MQ=255)
               cGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGc                                                                                    >  1:3320791/1‑101 (MQ=255)
                gCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCg                                                                                   >  1:6810353/1‑101 (MQ=255)
                gCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGGGCTGTCTTCCCACAGTGGGCg                                                                                   >  2:1193764/1‑101 (MQ=255)
                 ccTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTggnnn                                                                                  >  2:10772365/1‑98 (MQ=255)
                   tagCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCgnnn                                                                                >  2:11134701/3‑98 (MQ=255)
                   nnnccGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGt                                                                                <  2:7563988/98‑1 (MQ=255)
                    ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTg                                                                               >  1:12585047/1‑101 (MQ=255)
                    ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTg                                                                               <  1:55988/101‑1 (MQ=255)
                       cGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAg                                                                            >  1:7237981/1‑101 (MQ=255)
                        nnnttCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGa                                                                           <  2:14015244/98‑1 (MQ=255)
                        gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGannn                                                                           >  2:11530668/1‑98 (MQ=255)
                         gtttCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                          <  2:2099034/100‑1 (MQ=255)
                         aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACTGTTGGCGAGTGAAGAc                                                                          >  2:7101575/1‑101 (MQ=255)
                         aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                          >  2:69824/1‑101 (MQ=255)
                            tCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgc                                                                       <  1:9987225/101‑1 (MQ=255)
                                ccGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtg                                                                   <  1:6349394/101‑1 (MQ=255)
                                  gTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTc                                                                 >  2:3249776/1‑101 (MQ=255)
                                    nnngggCACCAGCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGt                                                               <  2:12852455/98‑1 (MQ=255)
                                    gCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGt                                                               <  2:1222765/101‑1 (MQ=255)
                                     nnnggCACCATCCACGGCAAAATGCCAATGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtt                                                              <  2:7821144/98‑1 (MQ=255)
                                        ggCCCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTAGCGAGTGAAGACCGCGGTGTCGTTTTc                                                           <  1:10548640/101‑1 (MQ=255)
                                        ggCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtttgc                                                           <  2:4673410/101‑3 (MQ=255)
                                         gCACAATCCACGGCAAAATGCCAACGCCATCCGGGAATCCCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCa                                                          <  2:2637164/101‑1 (MQ=255)
                                          cACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCaa                                                         >  1:6086311/1‑101 (MQ=255)
                                          cACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCaa                                                         <  2:807252/101‑1 (MQ=255)
                                           aCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAg                                                        >  2:6736147/1‑101 (MQ=255)
                                            tcATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGGGCTGTCTTCCCACAGTTGGCGAGGGAAGACCGCGGCGTCGTTTTCannn                                                       >  2:14170540/2‑98 (MQ=255)
                                            nnntCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCTGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                                                       <  2:13077396/98‑1 (MQ=255)
                                             cATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTa                                                      >  2:933197/1‑101 (MQ=255)
                                              aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTAc                                                     >  1:175189/1‑101 (MQ=255)
                                              aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTAc                                                     >  1:3958436/1‑101 (MQ=255)
                                               tacACGACAAAATCCCCACGCAGTCCGGGAATCCCTCCAGACACCAGGTGCTGTCTTCCCACACTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACa                                                    <  2:4000989/99‑1 (MQ=255)
                                                 cACGGCAAAATGCGAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                                                  >  1:1606247/1‑101 (MQ=255)
                                                 cACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                                                  <  1:4564477/101‑1 (MQ=255)
                                                  aCGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAg                                                 >  2:555835/1‑101 (MQ=255)
                                                    ggCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATannn                                               >  2:9655559/1‑98 (MQ=255)
                                                    cgCAAAATGCCACCGCCATCCGGGACTACCTCCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGt                                               <  2:2331422/100‑1 (MQ=255)
                                                     gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTg                                              <  1:8340622/101‑1 (MQ=255)
                                                     gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTg                                              >  1:589769/1‑101 (MQ=255)
                                                     gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTg                                              >  1:7793957/1‑101 (MQ=255)
                                                       nnnaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAg                                            <  2:14459494/98‑1 (MQ=255)
                                                       aaaaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTACCAGTACATAGGagag                                            >  2:1908562/1‑101 (MQ=255)
                                                         aaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAggg                                          <  2:4312156/101‑1 (MQ=255)
                                                          aTGCCAACGCCACCCGGGAATACCACCAGACGCTCAGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAGCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGc                                         <  1:11960435/101‑1 (MQ=255)
                                                            gCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGGGTCGTTTTCAAGTACATAGGTGAGGGCGa                                       >  1:9443505/1‑101 (MQ=255)
                                                            cccAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGCCGa                                       <  1:12966951/100‑1 (MQ=255)
                                                             nnnaCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGAt                                      <  2:13902523/98‑1 (MQ=255)
                                                             ccAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGAt                                      >  1:8657613/1‑101 (MQ=255)
                                                              cAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc                                     >  1:11537250/1‑101 (MQ=255)
                                                               caCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCa                                    <  2:4884834/100‑1 (MQ=255)
                                                                    nnnccGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtg                               <  2:7448871/98‑1 (MQ=255)
                                                                    cATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtg                               >  1:2520378/1‑101 (MQ=255)
                                                                     aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtgg                              >  1:7361567/1‑101 (MQ=255)
                                                                     aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGGGAGGGCGATCAGGtnnn                              >  2:9067094/1‑98 (MQ=255)
                                                                     aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGGTGGCGAGTGAAGACCGCGGTGTCGCGTTCAAGTACATAGGAGAGGGCGATCAGGTtgg                              >  1:2221121/1‑101 (MQ=255)
                                                                      tCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggt                             >  1:5089089/1‑101 (MQ=255)
                                                                       nnnggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtg                            <  2:10292392/98‑1 (MQ=255)
                                                                       ccGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtg                            <  2:1568110/101‑1 (MQ=255)
                                                                        cGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtgg                           >  1:12278650/1‑101 (MQ=255)
                                                                         nnnaaTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggc                          <  2:8854341/98‑1 (MQ=255)
                                                                         gggCATACCACCAGACCCTCGGTGCTGTCTTCCCAAAGTTGGCGAGTGACGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggc                          <  2:5919071/101‑1 (MQ=255)
                                                                         gggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAATACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggc                          <  1:7536508/101‑1 (MQ=255)
                                                                         gggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTCCATAGGTGAGGGCGATCAGGTTGgtggc                          <  1:3572285/101‑1 (MQ=255)
                                                                          ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCGGGTtggggnnn                         >  2:13559352/1‑96 (MQ=255)
                                                                          ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcg                         >  1:13588433/1‑101 (MQ=255)
                                                                          ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcg                         >  1:8176206/1‑101 (MQ=255)
                                                                             aTACCCCCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCACGTACATAGGTGAGGGCGATCAGGTTGgtggcgtag                      <  2:3433950/101‑3 (MQ=255)
                                                                             aTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtgg                      >  1:2696893/1‑101 (MQ=255)
                                                                              tACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggc                          <  1:10146625/96‑1 (MQ=255)
                                                                                ccaccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGCGAGGGCGATCAGGTTGGTGGCGTGGCAg                   <  1:11774591/101‑1 (MQ=255)
                                                                                 caccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAgg                  >  1:7011347/1‑100 (MQ=255)
                                                                                   nnngACACTCGGTGCTGTCTTCCCACAGTTGGCGACTGAAGTCCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGtac                <  2:10560501/98‑3 (MQ=255)
                                                                                   ccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGc                <  1:13875020/101‑1 (MQ=255)
                                                                                    caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCa               >  1:12806551/1‑101 (MQ=255)
                                                                                    caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCa               >  2:4513006/1‑101 (MQ=255)
                                                                                    caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCa               >  2:3702121/1‑101 (MQ=255)
                                                                                    caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCa               >  1:5369772/1‑101 (MQ=255)
                                                                                    caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCa               >  1:910077/1‑101 (MQ=255)
                                                                                       acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCannn            >  2:12588212/1‑98 (MQ=255)
                                                                                       acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCannn            >  2:11351761/1‑98 (MQ=255)
                                                                                        cacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGAt           <  1:10861898/101‑1 (MQ=255)
                                                                                         nnncGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATc          <  2:8328036/98‑1 (MQ=255)
                                                                                         acTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATc          <  1:7143815/101‑1 (MQ=255)
                                                                                           tCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCAc        <  2:5120784/101‑1 (MQ=255)
                                                                                              gTGCTGTCTTCCCACAGTTGGCGAGTGAAGTCCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCg     <  1:14194626/101‑1 (MQ=255)
                                                                                              gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCg     >  1:284496/1‑101 (MQ=255)
                                                                                               tgctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGa    <  1:11134701/101‑1 (MQ=255)
                                                                                               tgctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGa    >  1:13224195/1‑101 (MQ=255)
                                                                                                gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGAt   <  1:4434937/101‑1 (MQ=255)
                                                                                                 ctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATc  <  1:930050/101‑1 (MQ=255)
                                                                                                    |                                                                                                 
CATCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATC  >  NC_000913/4091693‑4091890

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: