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% 1020.4 / NA 264G168D (GGC→GAC) rhaDrhamnulose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (119/144);  minor base G (0/1);  total (119/145)
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 282 total aligned reads displayed.

TCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATC  >  NC_000913/4091695‑4091890
                                                                                                  |                                                                                                 
tctcTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCt                                                                                                 >  1:5767089/1‑101 (MQ=255)
 ctctTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCtt                                                                                                <  1:6700337/101‑1 (MQ=255)
 ctctTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCtt                                                                                                <  1:360046/101‑1 (MQ=255)
  tctTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTc                                                                                               >  1:624188/1‑101 (MQ=255)
   ttagtgCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTcc                                                                                              <  2:1296952/98‑1 (MQ=255)
    ttGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTccc                                                                                             >  1:5882817/1‑101 (MQ=255)
      gtgCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCcac                                                                                           >  2:2234236/1‑101 (MQ=255)
      gtgCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCcac                                                                                           >  2:5705642/1‑101 (MQ=255)
        gCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAg                                                                                         <  2:1398981/101‑1 (MQ=255)
        gCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAg                                                                                         <  1:2955796/101‑1 (MQ=255)
          ggTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                                       >  2:2757369/1‑101 (MQ=255)
          ggTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                                       >  1:2649239/1‑101 (MQ=255)
          ggTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                                                       >  1:1958715/1‑101 (MQ=255)
           gTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTg                                                                                      >  2:2139591/1‑101 (MQ=255)
             ggCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGACTACCCCCAGACCCTCGGTGCCGTCTTCCCACAGTTGGc                                                                                    <  2:7022142/100‑1 (MQ=255)
             cGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGc                                                                                    <  1:4108349/101‑1 (MQ=255)
              gcttGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGTAATACCACCAGACACTCGGTGCTGTCTTCCCCCAGTTGGCg                                                                                   <  1:6787816/98‑1 (MQ=255)
              gCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCg                                                                                   <  2:658358/101‑1 (MQ=255)
              gCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCg                                                                                   >  2:4674270/1‑101 (MQ=255)
               ccTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                                                  >  2:1593574/1‑101 (MQ=255)
                cTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAg                                                                                 >  1:4903441/1‑101 (MQ=255)
                 tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGt                                                                                >  1:3787208/1‑101 (MQ=255)
                 tGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGt                                                                                <  1:5296229/101‑1 (MQ=255)
                  ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTg                                                                               >  2:6033915/1‑101 (MQ=255)
                  ggCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCAGGGAATACCACCAGACACTCGTTGCTGTCTTCCCACAGTTGGCGAGTg                                                                               >  1:2163929/1‑101 (MQ=255)
                    ccGATTTCGCCCGTCCCGGCCACCATCCACGGCAAAATGCCAACACCATCCGGGAATACCACCAGACACTCGGTGCTGGCTGCCCACAGTTGGCGAGTGaa                                                                             <  2:6751372/101‑1 (MQ=255)
                     cGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAg                                                                            <  1:6539211/101‑1 (MQ=255)
                      gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGa                                                                           >  2:4492581/1‑101 (MQ=255)
                      gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGa                                                                           >  2:4925938/1‑101 (MQ=255)
                      gATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGa                                                                           <  1:2504893/101‑1 (MQ=255)
                       aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                          <  2:6195427/101‑1 (MQ=255)
                       aTTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                                          <  2:240586/101‑1 (MQ=255)
                        tttCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAcc                                                                         <  2:364482/101‑1 (MQ=255)
                        tttCGTCCGTCCCGGGCACCATCCACGGCAAAATGACAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAcc                                                                         <  2:3607283/101‑1 (MQ=255)
                         ttCGTCCGTGCCGGGCACCATCCCCGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcg                                                                        <  1:5202243/101‑1 (MQ=255)
                         ttCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGACTAACACCAGACACTCGGCGCTGTCTTCCCACAGTTGGCGAGTGAAGACcg                                                                        >  2:6406054/1‑101 (MQ=255)
                         ttCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcg                                                                        <  2:4853542/101‑1 (MQ=255)
                         ttCGGCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcg                                                                        <  2:5403106/101‑1 (MQ=255)
                          tCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACGCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgc                                                                       <  1:2443875/101‑1 (MQ=255)
                            gTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCCCCAGACACTCGGTGCTGTCTTCCCACAGCTGGCGAGTGAAGGCcgctg                                                                     <  2:4167142/101‑3 (MQ=255)
                            gTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACAATCGGTGCTGTCTTCCCACAGTTGGCGAGTGACGACCGCgg                                                                     >  1:6636723/1‑101 (MQ=255)
                             tCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgt                                                                    <  1:5302292/101‑1 (MQ=255)
                                 tGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCg                                                                <  1:4960884/101‑1 (MQ=255)
                                 tGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCg                                                                <  1:790746/101‑1 (MQ=255)
                                   ccGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtt                                                              <  2:2949179/101‑1 (MQ=255)
                                   ccGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtt                                                              >  2:3813337/1‑101 (MQ=255)
                                    cGGGCACCATCCACGGCAAAATGCCACCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGAGAGTGAAGACCGCGGTGTCGttt                                                             >  1:5919589/1‑101 (MQ=255)
                                    cGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                                             >  2:2534224/1‑101 (MQ=255)
                                    cGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                                             <  2:3131504/101‑1 (MQ=255)
                                    cGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                                             >  1:6395780/1‑101 (MQ=255)
                                    cGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                                             >  1:421875/1‑101 (MQ=255)
                                     gggCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtttt                                                            >  2:204388/1‑101 (MQ=255)
                                      ggCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTc                                                           <  1:3322803/101‑1 (MQ=255)
                                         aCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAg                                                        <  1:6972284/101‑1 (MQ=255)
                                         aCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAg                                                        <  1:5318014/101‑1 (MQ=255)
                                           cATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGCGAAGACCGCGGTGTCGTTTTCAAGTa                                                      >  1:5516789/1‑101 (MQ=255)
                                            aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTAc                                                     >  2:623363/1‑101 (MQ=255)
                                            aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGAGTCGTTTTCAAGTAc                                                     >  2:5217390/1‑101 (MQ=255)
                                             tCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACa                                                    >  2:2913592/1‑101 (MQ=255)
                                              ccACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACAt                                                   >  2:5370721/1‑101 (MQ=255)
                                               cACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                                                  >  2:1171119/1‑101 (MQ=255)
                                                tcGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATan                                                 <  2:1029507/100‑2 (MQ=255)
                                                aCGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAg                                                 >  2:5866487/1‑101 (MQ=255)
                                                  ggCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGt                                               <  1:6897759/101‑1 (MQ=255)
                                                   gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTg                                              <  2:1872315/101‑1 (MQ=255)
                                                        ntGCCAACGCCATCCGGGAATGCCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGc                                         <  2:5622147/100‑1 (MQ=255)
                                                         tGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCg                                        >  2:6109730/1‑101 (MQ=255)
                                                           ccAACGTCATCCGGGATTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGCTTGAAGCACGAGGTGCAGTTTTCAAGTACATAGGTGAGGGCGAt                                      <  2:6933461/101‑1 (MQ=255)
                                                            ctaagCCATCCGGGAATACCACCAGACACTCGGTGCTGGCTTCCCACAGTTGGCGAGTGAAGAGCGCGGTGTCGGTTTCAAGTACATAGGTGAGGGCGATc                                     >  2:3466573/5‑101 (MQ=255)
                                                            cAAGGCCACCCGGGAATACCACCAGACCCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGGTTTCAAGTACATAGGTGAGGGCGATc                                     <  2:646108/101‑1 (MQ=255)
                                                              aCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTCCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAg                                   <  1:100175/101‑1 (MQ=255)
                                                                 ccATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGtt                                <  1:1544752/101‑1 (MQ=255)
                                                                   aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtgg                              >  2:5140897/1‑101 (MQ=255)
                                                                     ccGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtg                            <  1:1144636/101‑1 (MQ=255)
                                                                      cGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtgg                           >  2:4739828/1‑101 (MQ=255)
                                                                       gggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggc                          >  2:4818795/1‑101 (MQ=255)
                                                                          aaTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtg                       <  2:885835/101‑1 (MQ=255)
                                                                            tACCACCAGACTCTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtggc                     <  1:3193714/101‑1 (MQ=255)
                                                                            tACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGgtggcgtggc                     >  2:4334128/1‑101 (MQ=255)
                                                                             accaccAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCa                    <  1:5369084/101‑1 (MQ=255)
                                                                             accaccAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCa                    <  1:4594519/101‑1 (MQ=255)
                                                                               caccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGt                  <  1:902614/101‑1 (MQ=255)
                                                                                  caGACACTCGGTGCTGTCTTCCCACGGTTGGCGAGTGAAGACCGCGGGGTCGTTTTCAAGTACATAGGCGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCa               >  1:2895775/1‑101 (MQ=255)
                                                                                   aGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCAt              <  2:6102407/101‑1 (MQ=255)
                                                                                     acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGa            <  2:6738885/101‑1 (MQ=255)
                                                                                     acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGa            <  2:5767089/101‑1 (MQ=255)
                                                                                     acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGa            <  1:4739588/101‑1 (MQ=255)
                                                                                      cacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGAt           <  1:2037895/101‑1 (MQ=255)
                                                                                       acTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATc          <  1:4691975/101‑1 (MQ=255)
                                                                                       acTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATc          <  1:6033915/101‑1 (MQ=255)
                                                                                        cTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCa         <  1:6481589/101‑1 (MQ=255)
                                                                                          cGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCAcc       >  1:4442935/1‑101 (MQ=255)
                                                                                          cGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCAcc       <  1:6593419/101‑1 (MQ=255)
                                                                                           ggTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCAccc      <  2:805413/101‑1 (MQ=255)
                                                                                           ggTGCTGTCTTCCCACAGTTGGCGAGGGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGTAGTGCATGATCcccc      <  2:1991781/101‑1 (MQ=255)
                                                                                            gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCg     <  2:5131419/101‑1 (MQ=255)
                                                                                            gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCg     <  1:2561337/101‑1 (MQ=255)
                                                                                            gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCg     <  1:1406547/101‑1 (MQ=255)
                                                                                            gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCg     >  2:3230366/1‑101 (MQ=255)
                                                                                             tgctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGa    >  1:1822360/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
TCTCTTGTGCGGTCGCCTGGCCGATTTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGGCGTGGCAGTGCATGATCACCCGATC  >  NC_000913/4091695‑4091890

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