Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,814,174 T→A 21.5% intergenic (+20/+44) emrB → / ← luxS multidrug efflux system protein/S‑ribosylhomocysteine lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,814,1740TA21.5% 269.8 / 68.4 125intergenic (+20/+44)emrB/luxSmultidrug efflux system protein/S‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base T (39/59);  new base A (12/15);  total (51/74)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 128 total aligned reads displayed.

ATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGT  >  NC_000913/2814076‑2814273
                                                                                                  |                                                                                                   
atTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAtt                                                                                                    <  2:504321/100‑1 (MQ=255)
  tCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTg                                                                                                  <  2:2578050/100‑1 (MQ=255)
  tCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTg                                                                                                  >  2:887476/1‑100 (MQ=255)
     tCGTTCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGGCAGTTCAAATGAAc                                                                                               <  2:987433/100‑1 (MQ=255)
     tCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAAc                                                                                               <  1:2042408/100‑1 (MQ=255)
     tCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAAc                                                                                               <  2:2614041/100‑1 (MQ=255)
         ccTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGc                                                                                           <  2:906697/100‑1 (MQ=255)
         ccTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGc                                                                                           <  2:1442625/100‑1 (MQ=255)
         ccTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTgcc                                                                                           <  2:1970882/100‑3 (MQ=255)
         ccTGCTGAGGCTGGTGTGGGTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTACGCGAGTTCAAATGAACTGGc                                                                                           <  1:23086/100‑1 (MQ=255)
          ctgctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCt                                                                                          <  1:2258612/100‑1 (MQ=255)
          ctgctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCt                                                                                          <  1:2595324/100‑1 (MQ=255)
          ctgctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCACACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCt                                                                                          <  1:1559501/100‑1 (MQ=255)
          ctgctgGGGCTGGTGTGGGTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCt                                                                                          <  1:1670465/100‑1 (MQ=255)
           tgctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCtt                                                                                         <  1:2813892/100‑1 (MQ=255)
           tgctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCtt                                                                                         <  2:684824/100‑1 (MQ=255)
           tgctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCtt                                                                                         <  2:1149420/100‑1 (MQ=255)
           tgctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCtt                                                                                         <  2:1000509/100‑1 (MQ=255)
            gctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCttt                                                                                        <  2:1431375/100‑1 (MQ=255)
              tgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGTCGCAGGTGGCGGCGGATGCTGTGCTCACTAAGTTCAACTAAGCCAGTTCAAATGAACTGGCttagt                                                                                      <  2:2771124/100‑4 (MQ=37)
               ggggCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCtttttt                                                                                     <  1:2302733/100‑1 (MQ=255)
               ggggCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCtttttt                                                                                     >  2:475497/1‑100 (MQ=255)
               ggggCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCtttttt                                                                                     <  2:1106782/100‑1 (MQ=255)
                gggCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCttagttg                                                                                    <  2:371607/100‑6 (MQ=38)
                  gCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCGAGTTCATTTGAACTGGCGTTTTTCa                                                                                   >  1:1906979/1‑99 (MQ=255)
                   cTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaat                                                                                 <  2:2708556/100‑1 (MQ=255)
                   cTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaat                                                                                 <  1:2781291/100‑1 (MQ=255)
                   cTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCttagttgtac                                                                                 <  2:2095884/100‑9 (MQ=38)
                    tGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaatt                                                                                <  2:2537445/100‑1 (MQ=255)
                     ggTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaatta                                                                               >  2:1825061/1‑100 (MQ=255)
                       tgtgGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaattaat                                                                             >  2:851929/1‑100 (MQ=255)
                        gtgGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaattaatt                                                                            <  2:2692320/100‑1 (MQ=255)
                          ggTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATtgt                                                                          <  1:1184211/100‑1 (MQ=255)
                            tttGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGa                                                                        <  1:2226305/100‑1 (MQ=255)
                             ttGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGaa                                                                       <  2:1480110/100‑1 (MQ=255)
                              tGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAg                                                                      <  2:501921/100‑1 (MQ=255)
                               gCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGa                                                                     >  1:1055449/1‑100 (MQ=255)
                                cTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGAt                                                                    >  2:2482374/1‑100 (MQ=255)
                                  aaaCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAg                                                                  <  1:1494244/100‑1 (MQ=255)
                                    aCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGtt                                                                >  1:908653/1‑100 (MQ=255)
                                    aCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGtt                                                                <  2:2840431/100‑1 (MQ=255)
                                     ccGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGttt                                                               >  2:127108/1‑100 (MQ=255)
                                      cGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTa                                                              >  2:214018/1‑100 (MQ=255)
                                       gCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTAc                                                             >  1:1689103/1‑100 (MQ=255)
                                        ccATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACt                                                            >  2:110195/1‑100 (MQ=255)
                                         cATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTg                                                           <  1:1890178/100‑1 (MQ=255)
                                          aTTTGGCGCAGGTGGCGGCTGAGGCGGTGCGCATTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTg                                                           <  2:1679776/99‑1 (MQ=255)
                                          aTTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGa                                                          >  1:57311/1‑100 (MQ=255)
                                            ttGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACt                                                        <  2:777463/100‑1 (MQ=255)
                                            ttGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACt                                                        <  1:1326100/100‑1 (MQ=255)
                                             tGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTa                                                       >  1:1503116/1‑100 (MQ=255)
                                              ggCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAg                                                      >  2:864127/1‑100 (MQ=255)
                                              ggCGAAGGTGGCGGCGGAGGCGGTGAGCACTAATTACAACTAAGTCAGTTCATTTGACCTCGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTa                                                       <  1:2280979/99‑1 (MQ=255)
                                               gcgcAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGa                                                     <  2:2196204/100‑1 (MQ=255)
                                                    ggTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGc                                                >  2:2066820/1‑100 (MQ=255)
                                                      tGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAg                                              >  2:2204429/1‑100 (MQ=255)
                                                         cggcggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTCTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTc                                           >  1:1863470/1‑100 (MQ=255)
                                                          ggcggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTcc                                          >  1:2143174/1‑100 (MQ=255)
                                                           gcggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCt                                         <  2:2617804/100‑1 (MQ=255)
                                                           gcggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCt                                         <  1:1613138/100‑1 (MQ=255)
                                                            cggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTg                                        <  1:440672/100‑1 (MQ=255)
                                                              gAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCa                                      <  1:1163263/100‑1 (MQ=255)
                                                               aGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCaa                                     <  1:1373195/100‑1 (MQ=255)
                                                                ggCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAAc                                    <  2:2201164/100‑1 (MQ=255)
                                                                gcccgttcgcACTAAGTACAACTAAGACAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTCACTAGTTGTGCAGTTCCTGCAAc                                    <  1:1299942/93‑1 (MQ=255)
                                                                 gCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATGAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACt                                   <  1:972527/100‑1 (MQ=255)
                                                                    gTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtct                                >  1:880909/1‑100 (MQ=255)
                                                                     tGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtctc                               >  2:981610/1‑100 (MQ=255)
                                                                       cgcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCtt                             >  2:2022620/1‑100 (MQ=255)
                                                                       cgcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCtt                             >  2:2817490/1‑100 (MQ=255)
                                                                        gcACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCttt                            >  2:2620905/1‑100 (MQ=255)
                                                                         cACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTc                           >  2:756073/1‑100 (MQ=255)
                                                                             aaGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTTGATGTGTATTTCCTGCAACTTCTCTTTCGGc                        >  1:1937458/1‑99 (MQ=255)
                                                                             aaGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCa                       <  2:2415672/100‑1 (MQ=255)
                                                                               gTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGt                     >  1:1267144/1‑100 (MQ=255)
                                                                               attgaaaaaaGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGt                     >  2:1136462/9‑100 (MQ=38)
                                                                               attgaaaaaaGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGt                     >  2:51371/9‑100 (MQ=38)
                                                                                 tgaaaaaaGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGc                   >  1:2286202/7‑100 (MQ=38)
                                                                                 aCAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGc                   <  2:1538259/100‑1 (MQ=255)
                                                                                  gaaaaaaGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGcc                  >  2:1953340/6‑100 (MQ=38)
                                                                                  cAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTTACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGcc                  <  1:1593778/100‑1 (MQ=255)
                                                                                  cAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGcc                  <  1:2167262/100‑1 (MQ=255)
                                                                                   aaCTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCa                 <  1:475494/100‑1 (MQ=255)
                                                                                    aaaaaGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAg                >  1:1127507/4‑100 (MQ=255)
                                                                                     cTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAAATGTGCAGTTCCTGCAATTTCTCTTTCGGCAATGCCAGt               >  1:2663237/1‑100 (MQ=255)
                                                                                      tAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGtt              >  2:1948397/1‑100 (MQ=255)
                                                                                      aaaGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGtt              >  1:803223/2‑100 (MQ=255)
                                                                                       naGCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttc             >  1:2621037/2‑100 (MQ=255)
                                                                                       aaGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttc             <  2:2768178/100‑1 (MQ=255)
                                                                                       aaGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttc             <  2:100646/100‑1 (MQ=255)
                                                                                         gCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTTTGAAGATAGTTTACTTACTCGATGTGCAGTTCCTGCAACTTCTCTTTTGGCAGTGCTAgttgtt           >  1:2263913/1‑100 (MQ=255)
                                                                                         gCCAGTTCAATTGAACTGGCTTTTTTCAATTAATTGTGAAGACAGTTTACTGACTAGACGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttct            <  1:15711/99‑1 (MQ=255)
                                                                                         gCCAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGAAGTTCCTACAACTTCTCTTTCGCCAGTGCCAGttc             >  1:1682542/1‑98 (MQ=255)
                                                                                           cAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCg         >  1:948443/1‑100 (MQ=255)
                                                                                               tCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCGACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGc     <  1:1918206/100‑1 (MQ=255)
                                                                                                cATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCt    <  1:294769/100‑1 (MQ=255)
                                                                                                 aTTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTg   >  1:1478466/1‑100 (MQ=255)
                                                                                                 aTTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTg   <  2:1635798/100‑1 (MQ=255)
                                                                                                  tttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGt  >  2:2210691/1‑100 (MQ=255)
                                                                                                  tttGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGt  >  2:921256/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
ATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGT  >  NC_000913/2814076‑2814273

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