Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,502 C→G 16.3% intergenic (+19/‑113) tnaB → / → mdtL tryptophan transporter of low affinity/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,891,5020CG16.3% 521.3 / 87.4 192intergenic (+19/‑113)tnaB/mdtLtryptophan transporter of low affinity/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (89/70);  new base G (13/18);  total (104/88)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

Only 100 of 208 total aligned reads displayed.

CGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAG  >  NC_000913/3891404‑3891601
                                                                                                  |                                                                                                   
cGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAgccagcca                                                                                                    >  1:293440/1‑100 (MQ=255)
     gTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTACCGTGTCACCTACCTTTGGCCAAATCCTTCAAGAAGCCCGCCAt                                                                                                   >  2:1726641/1‑96 (MQ=255)
     gTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGc                                                                                               >  1:1584264/1‑100 (MQ=255)
       cATTCTCTTTTGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATGCGCTg                                                                                             <  2:1134549/100‑1 (MQ=255)
       cATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                             >  1:3514104/1‑100 (MQ=255)
         ttCTCTTTGGTATAATCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                           <  1:4139344/100‑1 (MQ=255)
         ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTACATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                           <  2:1234536/100‑1 (MQ=255)
         ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                           <  1:744589/100‑1 (MQ=255)
         ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                           <  1:1632663/100‑1 (MQ=255)
         ttCACTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCCACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                           <  2:3800295/100‑1 (MQ=255)
          tctcTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGct                                                                                          <  2:1120666/100‑1 (MQ=255)
           ctctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTGCAAGAAGCCAGCGAATGGCTGGctt                                                                                         <  2:3978896/100‑1 (MQ=38)
             ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                       <  1:589733/100‑1 (MQ=39)
             ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttct                                                                                       >  1:3202569/1‑100 (MQ=255)
              tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttctt                                                                                      <  1:2329341/100‑1 (MQ=255)
               ttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTg                                                                                     <  1:2544052/100‑1 (MQ=38)
               ttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTg                                                                                     <  2:3344117/100‑1 (MQ=255)
               ttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATACTTCAAGAAGACAGCCATTAGCTgccacattg                                                                                     >  1:3922453/1‑92 (MQ=255)
               ttGGGATAACCGTGATTTTGTGCTGGTTCGGGAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTGCAAGAAGCCAGCGAATGGCTGGcttct                                                                                       <  1:1223740/98‑1 (MQ=37)
                 ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaa                                                                                   <  1:2478870/100‑3 (MQ=38)
                  gTATAGCCGTGATTTTGTGCTGGGTAGGCAGCGTCGTTAACGTGTTACCTTAATTTTGCTAAAGTCTTCAAGAAGCCAGCGAATGGCTGGcttctt                                                                                      <  2:3849853/96‑1 (MQ=19)
                  gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaag                                                                                  <  1:3373527/100‑4 (MQ=38)
                  gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaag                                                                                  <  1:1416058/100‑4 (MQ=38)
                   tataACCGTGATTTTGTGCTGGTTCGGCAATGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCCTTCGCTGGCTTCTTGCctc                                                                                 >  1:550056/1‑100 (MQ=255)
                   tataACCGTGATTTTGTGCTGGTTCGGCAATGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCctc                                                                                 >  1:1000666/1‑100 (MQ=255)
                   tataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCctc                                                                                 >  2:2393255/1‑100 (MQ=255)
                   tataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCctc                                                                                 <  1:2682459/100‑1 (MQ=255)
                    ataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCctct                                                                                >  2:1069709/1‑100 (MQ=255)
                       aCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAg                                                                             <  1:46346/100‑1 (MQ=255)
                        ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAgg                                                                            >  1:3197086/1‑100 (MQ=255)
                        ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAgg                                                                            <  1:1662437/100‑1 (MQ=255)
                          gTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaa                                                                          <  1:2905978/100‑1 (MQ=255)
                           tGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaaa                                                                         <  1:1733735/100‑1 (MQ=255)
                              ttttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCa                                                                      >  1:3255732/1‑100 (MQ=255)
                               tttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCAc                                                                     <  1:771320/100‑1 (MQ=255)
                               tttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCAc                                                                     <  1:3953315/100‑1 (MQ=255)
                                 tgtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCa                                                                      >  1:3601664/1‑97 (MQ=255)
                                  gtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTa                                                                  <  2:860275/100‑1 (MQ=255)
                                  gtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTa                                                                  <  2:3593258/100‑1 (MQ=255)
                                     cTGGTTCGTCAACGTCTTGAACGTGTTACCTAAATTTGGCTAAATCCTTGAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGt                                                               <  2:1313724/100‑1 (MQ=255)
                                     cTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGt                                                               <  1:2257798/100‑1 (MQ=255)
                                      tGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTc                                                              >  1:1435072/1‑100 (MQ=255)
                                       ggTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTcc                                                             >  1:1604776/1‑100 (MQ=255)
                                       ggTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTcc                                                             >  1:406234/1‑100 (MQ=255)
                                           cGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAAt                                                         >  2:213336/1‑100 (MQ=255)
                                            ggCAACGCCGTTAACGAGTTACCTAAATTTGGCGAAAGCCTTCAAGAAGCCAGACTGTCGCTGGCTACTTGCCTCTCAGGAAATCACTTATGTCCaaa                                                          >  2:671312/1‑98 (MQ=255)
                                             gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                       >  1:3800756/1‑100 (MQ=255)
                                              cAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                      >  2:3893931/1‑100 (MQ=255)
                                              cAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                      <  2:4048329/100‑1 (MQ=255)
                                               aaCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCa                                                     >  1:2424440/1‑100 (MQ=255)
                                                   tCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTc                                                 >  1:3408269/1‑100 (MQ=255)
                                                   tCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTc                                                 >  2:3706215/1‑100 (MQ=255)
                                                      ttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGcc                                              <  2:3005718/100‑1 (MQ=255)
                                                       tAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCt                                             >  2:3525036/1‑100 (MQ=255)
                                                       tAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCt                                             >  2:1168144/1‑100 (MQ=255)
                                                        aaCGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTg                                            <  2:1907202/100‑1 (MQ=255)
                                                            tgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcc                                        <  1:4023642/100‑1 (MQ=255)
                                                            tgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcc                                        <  1:2042384/100‑1 (MQ=255)
                                                              ttACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCAGAtcctc                                      >  2:2249656/1‑100 (MQ=255)
                                                               tACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                     >  1:821561/1‑100 (MQ=255)
                                                                aCCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATcctcct                                    >  2:1953006/1‑100 (MQ=255)
                                                                aCCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATcctcct                                    >  1:1331212/1‑100 (MQ=255)
                                                                aCCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATcctcct                                    <  1:2215331/100‑1 (MQ=255)
                                                                 ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                   >  1:2749356/1‑100 (MQ=255)
                                                                 ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                   >  1:3291720/1‑100 (MQ=255)
                                                                  cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTc                                  >  1:1783531/1‑100 (MQ=255)
                                                                   tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTc                                  <  2:827089/99‑1 (MQ=255)
                                                                     aaTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcacc                               <  2:3276454/100‑1 (MQ=255)
                                                                     aaTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcacc                               >  1:2724660/1‑100 (MQ=255)
                                                                        ttGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACg                            >  1:1939646/1‑100 (MQ=255)
                                                                        ttGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACg                            >  1:1474624/1‑100 (MQ=255)
                                                                          ggCTAAATCCTTCAAGAAGCCAGCTATTCGCTGGCTTCTTTCCTCTCAGGAAATTACTTATGTCCAACTGGCAACTCGCCTGATCCTCCTTCACCACGTa                          >  2:1866768/1‑100 (MQ=255)
                                                                          ggCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTa                          <  2:499176/100‑1 (MQ=255)
                                                                          ggCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTa                          <  1:3893956/100‑1 (MQ=255)
                                                                           gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTAt                         >  1:308048/1‑100 (MQ=255)
                                                                            cTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATg                        >  1:2164000/1‑100 (MQ=255)
                                                                            cTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATg                        <  2:2838882/100‑1 (MQ=255)
                                                                              aaaTCCTTCAAGAAGCTAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCt                      <  2:2447378/100‑1 (MQ=255)
                                                                              aaaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCCACTCGCCTGATCCTCCTTCACCACGTAt                         >  2:2605659/1‑97 (MQ=255)
                                                                              aaaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGCTCCCCACTCACCACGTATGc                       >  2:3224736/1‑99 (MQ=255)
                                                                               aaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCtt                     >  1:3324572/1‑100 (MQ=255)
                                                                                aTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCttt                    >  1:3996845/1‑100 (MQ=255)
                                                                                  ccTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                  <  2:2905718/100‑1 (MQ=255)
                                                                                   cTTCAAGAAGCCAGCGATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCg                 <  1:3783554/100‑1 (MQ=255)
                                                                                    ggcAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGt                >  1:540010/3‑100 (MQ=38)
                                                                                    ggcAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGt                >  2:2821717/3‑100 (MQ=38)
                                                                                     tCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc               <  1:2958743/100‑1 (MQ=255)
                                                                                      cAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATTGCGACTCGCATGATCCTCCTTCACCACGTATGCTTTGCGTCa              >  2:2939841/1‑100 (MQ=37)
                                                                                      cAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa              >  2:1592108/1‑100 (MQ=255)
                                                                                      cAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa              >  1:3394654/1‑100 (MQ=255)
                                                                                       aagaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAc             >  2:4006705/1‑100 (MQ=39)
                                                                                       aagaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAc             >  1:2161097/1‑100 (MQ=39)
                                                                                        agaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCGAATGGCAACTCGCCTGATCCTCCTTCACTACGTATGCTTTGCGTCAc             >  2:2706975/1‑99 (MQ=255)
                                                                                         gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt           >  2:1398727/1‑100 (MQ=255)
                                                                                         gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt           >  1:997172/1‑100 (MQ=255)
                                                                                           agccagccATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTa         >  2:1770363/1‑100 (MQ=255)
                                                                                            gccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTAc        >  2:3061974/1‑100 (MQ=255)
                                                                                            gccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTAc        >  2:2334043/1‑100 (MQ=255)
                                                                                              cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACt       <  2:1418191/99‑1 (MQ=255)
                                                                                                  caTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAg  <  1:3738333/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAG  >  NC_000913/3891404‑3891601

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: