Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,568,688 C→A 52.6% 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,6880CA52.6% 4.5 / 170.9 121intergenic (‑169/‑310)mdtM/yjiPmultidrug efflux system protein/pseudogene, transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base A (42/21);  total (69/52)
Fisher's exact test for biased strand distribution p-value = 4.22e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

GGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAATTAAAA  >  NC_000913/4568590‑4568786
                                                                                                  |                                                                                                   
gacacttGATAAGGGCAACGAGACTGCAAATGAACTGAAGTTCATGTCATCAACCACCCTGCCTGTGAACACGCAGACATATAAATACTTTCTGAGTTCg                                                                                                    <  2:2247831/94‑1 (MQ=255)
  cAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCggg                                                                                                  <  2:1212658/100‑1 (MQ=255)
  cAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAggg                                                                                                  >  1:1281127/1‑100 (MQ=255)
  cAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAggg                                                                                                  >  2:1909601/1‑100 (MQ=255)
   aGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGt                                                                                                 <  1:1367790/100‑1 (MQ=255)
     gTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCa                                                                                               >  2:1734358/1‑100 (MQ=255)
       gagaAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGa                                                                                             >  1:915465/1‑100 (MQ=255)
         gaAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAAt                                                                                           >  2:2312378/1‑100 (MQ=255)
          taGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATg                                                                                          <  1:2352916/99‑1 (MQ=255)
          aaGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATg                                                                                          <  1:913828/100‑1 (MQ=255)
          aaGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATg                                                                                          >  1:2387874/1‑100 (MQ=255)
          aaGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATg                                                                                          <  2:2280699/100‑1 (MQ=255)
           aGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATga                                                                                         >  2:2352888/1‑99 (MQ=255)
            gggCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGta                                                                                        <  2:278271/100‑1 (MQ=255)
            gggCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGta                                                                                        <  1:1780929/100‑1 (MQ=255)
             ggCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGta                                                                                        <  1:1273433/99‑1 (MQ=255)
             cgCAACGAGACTGCAAATGAACTGAATATCATGAAATCAACAACACAGACAGTGAACACACAGATAGAGAAATACTTTGTGACTTCTTATCAGAATGtat                                                                                       >  2:1254986/2‑100 (MQ=255)
               cAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATaa                                                                                     >  2:1851105/1‑100 (MQ=255)
                aaCGAGACTGCAAATGTACTGACTTTCATGTCATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAAc                                                                                    <  2:1618165/100‑1 (MQ=255)
                aaCGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAAc                                                                                    >  2:1553766/1‑100 (MQ=255)
                 aCGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACg                                                                                   <  2:547635/100‑1 (MQ=255)
                   gagaCTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGAt                                                                                 <  2:2218580/100‑1 (MQ=255)
                    agaCTGTAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGAtt                                                                                <  1:582378/100‑1 (MQ=255)
                      aCTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGtatttctctgtc                                                                              >  1:1376944/1‑91 (MQ=255)
                          cAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAAt                                                                         >  2:830730/1‑100 (MQ=255)
                          cAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAAt                                                                         >  1:1269977/1‑100 (MQ=255)
                          cAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAAt                                                                         >  2:1747906/1‑100 (MQ=255)
                          cAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAAt                                                                         >  2:422499/1‑100 (MQ=255)
                            aaTGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAg                                                                       <  2:1655664/100‑1 (MQ=255)
                             aTGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGt                                                                      <  2:2048748/100‑1 (MQ=255)
                                 aCTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGtcatc                                                                  <  2:862208/100‑1 (MQ=255)
                                   tGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAg                                                                <  2:2011330/100‑1 (MQ=255)
                                    tccgatcttGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGc                                                               <  2:2021291/92‑1 (MQ=255)
                                     aaTTTCATGTAATCAACCACACTGCATGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGca                                                              <  1:598223/100‑1 (MQ=255)
                                      aTTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGcac                                                             >  2:36547/1‑100 (MQ=255)
                                      aTTTCATGTAATCAACCACACTACCTGTGAACACACAGACAGAGAAATAATTTCTAACTTCGGGTCAGAATGTATAACGATTAGAAATAGTCATCAacac                                                             >  2:1541383/1‑100 (MQ=255)
                                         tCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACAtt                                                          <  2:779698/100‑1 (MQ=255)
                                          cATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTa                                                         >  2:1032162/1‑100 (MQ=255)
                                          cATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTa                                                         >  1:674164/1‑100 (MQ=255)
                                          cATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTa                                                         >  2:523053/1‑100 (MQ=255)
                                           aTGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTac                                                        >  2:1995235/1‑100 (MQ=255)
                                            tGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGcagatcgga                                                       >  1:2021313/1‑93 (MQ=255)
                                              tAATCCACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACAcc                                                     >  2:1485177/1‑100 (MQ=255)
                                              tAATCAACCACACTGCCTGTGAACACACAGACAGAGAATTACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACAcc                                                     >  2:1473671/1‑100 (MQ=255)
                                              tAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACAcc                                                     >  2:1917671/1‑100 (MQ=255)
                                              tAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACAcc                                                     >  1:537131/1‑100 (MQ=255)
                                                aTCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACct                                                    >  1:661725/1‑99 (MQ=255)
                                                aTCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACct                                                    >  1:1879154/1‑99 (MQ=255)
                                                  cAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCtt                                                 <  1:1344846/100‑1 (MQ=255)
                                                    aCCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGa                                               <  1:1032173/100‑1 (MQ=255)
                                                     ccACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGaa                                              <  1:1611366/100‑1 (MQ=255)
                                                     ccACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGaa                                              <  1:2182386/100‑1 (MQ=255)
                                                       acacTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATATTCATCAGCACATTACACCTCTTGAAtt                                            >  2:872559/1‑100 (MQ=255)
                                                         acTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGa                                          >  2:1061550/1‑100 (MQ=255)
                                                          cTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACATCTTGAATTGAt                                         >  1:1164661/1‑100 (MQ=255)
                                                           tGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATg                                        <  2:707841/100‑1 (MQ=255)
                                                           tGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATg                                        <  1:2249283/100‑1 (MQ=255)
                                                            gCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATgg                                       >  2:482880/1‑100 (MQ=255)
                                                             ccTGTGAACCCAAAGAAAGAGAAAAAATTTCTGAGGTAGGCACAGAATGTATAACGATTACAAATAGTCATCAACACATTACAGCTCTGGAATTGATggg                                      >  2:245237/1‑100 (MQ=255)
                                                             ccTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATgg                                       >  1:219027/1‑99 (MQ=255)
                                                              cTGTGAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAAAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTAATggg                                      >  2:2180265/1‑99 (MQ=255)
                                                               tgtgAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGtt                                    >  2:915659/1‑100 (MQ=255)
                                                               tgtgAACACACAGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGtt                                    >  1:230761/1‑100 (MQ=255)
                                                                gtgAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTg                                   >  1:5492/1‑100 (MQ=255)
                                                                 tgAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGc                                  <  2:272508/100‑1 (MQ=255)
                                                                  gAACAAACAGACAGAGAAATTCTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGc                                  <  1:204340/100‑1 (MQ=255)
                                                                     cacacaGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAg                              >  2:1333464/1‑100 (MQ=255)
                                                                     cacacaGACAGAGAAATACTTTCTGACGTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAg                              <  1:278759/100‑1 (MQ=255)
                                                                       cacaGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGaa                            >  2:601691/1‑100 (MQ=255)
                                                                       cacaGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGTTTGCATAGaa                            >  1:707314/1‑100 (MQ=255)
                                                                       cacaGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGaa                            >  1:849540/1‑100 (MQ=255)
                                                                       cacaGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGaa                            >  2:1111877/1‑100 (MQ=255)
                                                                       cacaGACAGAGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGaa                            >  1:1547523/1‑100 (MQ=255)
                                                                         cagacagaGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGaaaa                          >  2:70446/1‑100 (MQ=255)
                                                                          agacagaGGAATACTTTCTGACTTCGGGTCAGAATATATAACGATTACAAATAGTCATCATCATATTACACCTCTTGATTTGATGGGTTGCATAGaaan                          <  2:88761/99‑2 (MQ=255)
                                                                          agacagaGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAat                         <  1:673167/100‑1 (MQ=255)
                                                                            acagaGAAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAata                        <  1:338426/99‑1 (MQ=255)
                                                                             cagaGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAAtata                      >  2:689031/1‑100 (MQ=255)
                                                                             cagaGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAAtata                      >  1:2204815/1‑100 (MQ=255)
                                                                                agaAATACTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAACATATacac                   >  2:1519616/1‑100 (MQ=255)
                                                                                   aaTACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATTGGTTGCATAGAAAATATACACCtt                >  2:1981364/1‑100 (MQ=255)
                                                                                    aTACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTa               <  2:63563/100‑1 (MQ=255)
                                                                                      aCTTTCTGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTaan             <  2:1596177/100‑2 (MQ=255)
                                                                                       cTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTTACGAGCACATTACACCTCTTGAATTAATGGGTTGGAGAGATAATATACACCtt                >  2:2098273/1‑96 (MQ=255)
                                                                                        tttCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTaa              >  1:705623/1‑97 (MQ=255)
                                                                                         ttCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAgtgt          <  1:447465/100‑1 (MQ=255)
                                                                                         ttCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAgtgt          >  1:1800513/1‑100 (MQ=255)
                                                                                           cTGACTTAGGGTCAGAATTTATAACGATTACAAATAGTCATCAGCAAATGACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAgtgt          <  2:189155/98‑1 (MQ=255)
                                                                                            tGATTTCGTGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTTAATTGATAGGTTGTATAGAAAATATACACCTTAAGTGTAAt       <  2:1609980/100‑1 (MQ=255)
                                                                                            tGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAAt       <  1:915670/100‑1 (MQ=255)
                                                                                            tGACTTCGGATCAGAATGTATAACGATTACAAATAGTCATCAGAACATTACACCTCTTGAATTGATGGGTTGGATAGAAAATATACACCTTAAGTGTAAt       <  2:235991/100‑1 (MQ=255)
                                                                                            tGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAAt       >  2:1661858/1‑100 (MQ=255)
                                                                                            tGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAAt       <  1:1789996/100‑1 (MQ=255)
                                                                                            tGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAAt       <  1:1384880/100‑1 (MQ=255)
                                                                                            tGACTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAAt       <  1:1251006/100‑1 (MQ=255)
                                                                                              aCTTTGGGTCGGAATGTATAACGATTAGTAATAGGCATCAGCATATTACGCCTCTTGAATTGTTGCGTTGCATAGAACACATACGTCTTTGGTGTAAt       <  2:1904373/98‑1 (MQ=255)
                                                                                              aCTTCGGGTCATAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAATTa     <  2:2040650/100‑1 (MQ=255)
                                                                                              aCTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAATTa     >  2:136165/1‑100 (MQ=255)
                                                                                              aCTTAGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAATTa     >  1:47010/1‑100 (MQ=255)
                                                                                                 tCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAATTaaaa  >  1:515455/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAACACACAGACAGAGAAATACTTTCTGACTTCGGGTCAGAATGTATAACGATTACAAATAGTCATCAGCACATTACACCTCTTGAATTGATGGGTTGCATAGAAAATATACACCTTAAGTGTAATTAAAA  >  NC_000913/4568590‑4568786

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