Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,568,688 C→A 20.3% intergenic (‑169/‑310) mdtM ← / → yjiP multidrug efflux system protein/pseudogene, transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,568,6880CA20.3% 373.5 / 85.1 164intergenic (‑169/‑310)mdtM/yjiPmultidrug efflux system protein/pseudogene, transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base C (63/66);  new base A (19/14);  total (83/81)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.04e-01

Only 100 of 166 total aligned reads displayed.

AGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAATTAAAAT  >  NC_000913/4568589‑4568787
                                                                                                   |                                                                                                   
 ggCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCg                                                                                                    >  2:2039433/1‑100 (MQ=255)
  gCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCgg                                                                                                   >  2:1679856/1‑100 (MQ=255)
   cAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTACTCAACCACACTGCCTGTGAACACACAGACAGAGAGATACTTTCTGACTTCggg                                                                                                  <  1:64666/100‑1 (MQ=255)
    aGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAgg                                                                                                   >  1:2564277/1‑98 (MQ=255)
     ggTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCgg                                                                                                   >  1:2583884/1‑97 (MQ=255)
        gagaAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGa                                                                                             >  1:2015015/1‑100 (MQ=255)
        gagaAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGa                                                                                             >  1:178505/1‑100 (MQ=255)
           aaGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATg                                                                                          <  1:1261715/100‑1 (MQ=255)
           aaGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATg                                                                                          >  2:441478/1‑100 (MQ=255)
            aGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGt                                                                                         >  2:1233045/1‑100 (MQ=255)
            aGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGt                                                                                         >  1:593007/1‑100 (MQ=255)
            aGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGt                                                                                         >  1:350344/1‑100 (MQ=255)
             gggCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGta                                                                                        <  2:2909176/100‑1 (MQ=255)
             gggCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGta                                                                                        <  2:336984/100‑1 (MQ=255)
             gggCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGta                                                                                        >  1:548059/1‑100 (MQ=255)
             gggCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGta                                                                                        >  2:3088105/1‑100 (MQ=255)
                cAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATaa                                                                                     >  2:400511/1‑100 (MQ=255)
                    gagaCTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGAt                                                                                 <  2:813696/100‑1 (MQ=255)
                      gaCTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTa                                                                               <  2:2988819/100‑1 (MQ=255)
                        cTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACa                                                                             >  2:2544065/1‑100 (MQ=255)
                        cTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACa                                                                             >  2:3306092/1‑100 (MQ=255)
                        cTGCAAATGTACTGAATGTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTAc                                                                              >  1:2394540/1‑99 (MQ=255)
                         tGCAAATGTACTGAATTTCATGTAATCGACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACaa                                                                            <  1:1923594/100‑1 (MQ=255)
                           cAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAAt                                                                          <  1:36909/100‑1 (MQ=255)
                           cAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAAt                                                                          >  1:1654286/1‑100 (MQ=255)
                           cAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAAt                                                                          >  2:1896650/1‑100 (MQ=255)
                            aaaTGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATa                                                                         <  2:1584775/100‑1 (MQ=255)
                            aaaTGTACTGAATTTCATGGAATCAACCACACTTCCTGTTAACACACAGACAGAGAAATACTTTCTGACTGCGGGTCAGAATGTATAACGATTACAAATa                                                                         <  1:141227/100‑1 (MQ=255)
                             aaTGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAg                                                                        >  1:2788411/1‑100 (MQ=255)
                             aaTGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAg                                                                        <  2:2167129/100‑1 (MQ=255)
                             aaTGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAg                                                                        <  2:2216558/100‑1 (MQ=255)
                              aTGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGt                                                                       <  1:657030/100‑1 (MQ=255)
                               tGTACTGAATTACATGTAATCACCCACACTGCCTGTAAACACACAAACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATACtc                                                                      <  2:747561/100‑1 (MQ=255)
                                gTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGAGTCAGAATGTATAACGATTACAAATAGtca                                                                     >  1:539702/1‑100 (MQ=255)
                                   cTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGtcatca                                                                  <  1:2887982/100‑1 (MQ=255)
                                    tGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAg                                                                 <  1:2792791/100‑1 (MQ=255)
                                    tGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAg                                                                 <  1:868383/100‑1 (MQ=255)
                                     gAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGc                                                                >  2:581550/1‑100 (MQ=255)
                                      aaTTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGca                                                               <  1:1523627/100‑1 (MQ=255)
                                          tCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACAtt                                                           <  1:2465001/100‑1 (MQ=255)
                                          tCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACAtt                                                           >  2:2705905/1‑100 (MQ=255)
                                            aTGTAATCAACCACACTGCCTGTGAACACACAGACAAAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTac                                                         >  2:2392695/1‑100 (MQ=255)
                                               tAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACAcc                                                      <  2:2942483/100‑1 (MQ=255)
                                                   cAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCtt                                                  >  2:1003662/1‑100 (MQ=255)
                                                   cAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGCCAGAATGTATAAACATTACAAATAGTAATCAGCACATTaca                                                        >  1:326954/1‑94 (MQ=255)
                                                    aaCCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTg                                                 <  1:784064/100‑1 (MQ=255)
                                                    aaCCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTg                                                 <  1:146271/100‑1 (MQ=255)
                                                    aaCCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTg                                                 <  2:2803104/100‑1 (MQ=255)
                                                     aCCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGa                                                >  1:198892/1‑100 (MQ=255)
                                                       cacaCTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAAt                                              >  2:2273019/1‑100 (MQ=255)
                                                       cacaCTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAAt                                              <  2:1574950/100‑1 (MQ=255)
                                                       cacaCTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAAt                                              <  1:104780/100‑1 (MQ=255)
                                                        acacTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAACTAGTCATCAGCACATTACACCTCTTGAAtt                                             <  1:1975167/100‑1 (MQ=255)
                                                        acacTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAAtt                                             <  1:2518358/100‑1 (MQ=255)
                                                         cacTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTg                                            <  1:2744931/100‑1 (MQ=255)
                                                           cTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGAt                                          >  1:3229478/1‑100 (MQ=255)
                                                            tGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATg                                         <  1:2239124/100‑1 (MQ=255)
                                                             gCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATgg                                        >  1:465252/1‑100 (MQ=255)
                                                             gCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATgg                                        <  2:143196/100‑1 (MQ=255)
                                                             gCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATgg                                        >  2:1522764/1‑100 (MQ=255)
                                                              ccTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATggg                                       <  2:601443/100‑1 (MQ=255)
                                                              ccTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATggg                                       >  2:1268325/1‑100 (MQ=255)
                                                                 gtgAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTg                                    >  1:2947600/1‑100 (MQ=255)
                                                                 gtgAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTg                                    <  2:863828/100‑1 (MQ=255)
                                                                 gtgAACACACAGACAGAGAAATACTTTCTGACTTCGGGACAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTg                                    >  1:1891288/1‑100 (MQ=255)
                                                                  tgAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGc                                   <  1:1793241/100‑1 (MQ=255)
                                                                  tgAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGc                                   >  1:13596/1‑100 (MQ=255)
                                                                  tgAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGc                                   <  1:3305149/100‑1 (MQ=255)
                                                                   gAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCa                                  <  2:2788434/100‑1 (MQ=255)
                                                                    aaCACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCAt                                 <  1:1325043/100‑1 (MQ=255)
                                                                        cacaGACAGAGAAGTACTTTCTGCCTTAGGGTCAGAATGTATAACGATTCCAAATAGTCATCAGCACATTACACCTCTTGAATTGCTGGGTTACATCGa                              >  1:2096039/1‑99 (MQ=255)
                                                                        cacaGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGaa                             >  2:1872799/1‑100 (MQ=255)
                                                                           agacagaGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAat                          <  1:3335727/100‑1 (MQ=255)
                                                                            gacagaGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAata                         <  1:1159769/100‑1 (MQ=255)
                                                                            gacagaGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAata                         >  1:1955895/1‑100 (MQ=255)
                                                                            gacagaGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAata                         <  1:914437/100‑1 (MQ=255)
                                                                            gacagaGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAata                         <  1:1691004/100‑1 (MQ=255)
                                                                              cagaGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAAtata                       >  2:2593302/1‑100 (MQ=255)
                                                                                gagaAATACTTTCTGACTTCGGGTCAGAATGTGTAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATaca                     >  2:3212889/1‑100 (MQ=255)
                                                                                gagaAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATaca                     >  2:2430793/1‑100 (MQ=255)
                                                                                  gaAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACAcc                   >  1:2210834/1‑100 (MQ=255)
                                                                                     aTACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTa                <  2:1551306/100‑1 (MQ=255)
                                                                                       aCTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAg              <  2:256198/100‑1 (MQ=255)
                                                                                         tttCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAgtg            <  2:1469061/100‑1 (MQ=255)
                                                                                         tttCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAgtg            <  1:939088/100‑1 (MQ=255)
                                                                                         ttgttGACTTGGGATCATAATGTATAACGATTACAAATATTCATCAGCACATTACACCTTGTGAATTGATGGGTTGCATAGCAATTATACACCTTAAgtg            <  1:2081629/96‑1 (MQ=255)
                                                                                          ttCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAgtgt           <  1:572477/100‑1 (MQ=255)
                                                                                          ttCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAgtgt           >  1:1592183/1‑100 (MQ=255)
                                                                                            cTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTaa         >  1:970903/1‑100 (MQ=255)
                                                                                            cTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTaa         >  2:1205534/1‑100 (MQ=255)
                                                                                            cTGACTTAGGGTCAGAATGTATAACTATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTa          <  1:1364852/99‑1 (MQ=255)
                                                                                            cTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTaa         <  1:373800/100‑1 (MQ=255)
                                                                                             tGACTTCGGGTCAGAATGTATAACGATTACAAATATTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAAt        <  2:623762/100‑1 (MQ=255)
                                                                                             tGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAAt        >  1:54098/1‑100 (MQ=255)
                                                                                             tGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAAt        <  1:2966685/100‑1 (MQ=255)
                                                                                              gACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAAtt       >  2:2877321/1‑100 (MQ=255)
                                                                                               aCTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAAtt       >  1:1508448/1‑99 (MQ=255)
                                                                                                cTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAATTGTAATTaa     >  1:3242167/1‑100 (MQ=255)
                                                                                                  tCGGGTCAGAATGTATAACGATTACAAATAGTCGTCAGCACATTACACCTCTTTAATTGATGTGTTGCATAGAAAATATACACCTTAAGTGTAATTaaaa   <  1:3125387/100‑1 (MQ=255)
                                                                                                   cGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAATTAAAAt  >  1:1233373/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
AGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAATTAAAAT  >  NC_000913/4568589‑4568787

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