Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 481,968 A→G 43.1% G812G (GGT→GGC acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913481,9680AG43.1% 59.5 / 190.8 158G812G (GGT→GGCacrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (52/38);  new base G (34/34);  total (86/72)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.56e-01

Only 100 of 161 total aligned reads displayed.

TCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTCGCCGATATCATCCGGCAGCATACGG  >  NC_000913/481870‑482067
                                                                                                  |                                                                                                   
 ccATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAAcc                                                                                                   >  1:2427422/1‑100 (MQ=255)
  cATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCg                                                                                                  <  1:1640273/100‑1 (MQ=255)
  cATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCg                                                                                                  >  2:2169608/1‑100 (MQ=255)
  cATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCg                                                                                                  <  2:791221/100‑1 (MQ=255)
     ccgccccATTGCTTCACCGGTACTTTTACCCGGCGCCGCCTCTCCTACGATTTCCATGGATGGCACGCCGATGTAACGTTCCATACGCGGCGAACCGTAc                                                                                               >  1:455025/6‑100 (MQ=255)
     cAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTAc                                                                                               <  2:3728467/100‑1 (MQ=255)
     cAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTAc                                                                                               <  1:1949160/100‑1 (MQ=255)
        cggcATTGCTTCACCGGTACTTTTACCCGGGGCCGCCTGACCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTc                                                                                             <  1:3095431/96‑1 (MQ=255)
        cTCCCTTGCTTCACCGGTACTTTGACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTcc                                                                                            <  1:3813862/100‑1 (MQ=255)
        cTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTcc                                                                                            <  1:2169603/100‑1 (MQ=255)
        cTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTcc                                                                                            >  1:1077231/1‑100 (MQ=255)
         tCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTccc                                                                                           >  2:532677/1‑100 (MQ=255)
         tCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTccc                                                                                           >  2:405629/1‑100 (MQ=255)
           cATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCaa                                                                                         >  2:808094/1‑100 (MQ=255)
           aaTTGCATCACCGGTACTTTTACCCGGTGCCGCCTTGCATAAGCTTTCCATGGATGGCAGGCAGTTGTGACGGTCCAGACGCGACGAACCGTACTCCCa                                                                                          >  1:2881132/2‑99 (MQ=255)
               gCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACgag                                                                                     >  1:1555095/1‑100 (MQ=255)
               gCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACgag                                                                                     >  2:3412119/1‑100 (MQ=255)
                cTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACgaga                                                                                    >  1:3747797/1‑100 (MQ=255)
                cTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACgaga                                                                                    >  1:631127/1‑100 (MQ=255)
                cTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACgaga                                                                                    <  1:3119263/100‑1 (MQ=255)
                cTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACgaga                                                                                    <  1:2832250/100‑1 (MQ=255)
                    aCCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAagag                                                                                <  2:2511660/100‑1 (MQ=255)
                    aCCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAagag                                                                                <  2:3253695/100‑1 (MQ=255)
                    aCCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTCCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAagag                                                                                <  1:3452580/100‑1 (MQ=255)
                     ccGGTACTTTTACCCGGTGCCGCCTGGCCTAAGGTTTCCATGGATGGCACGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAAGGAGAAgagg                                                                               >  2:1048987/1‑100 (MQ=255)
                     ccGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAgagg                                                                               <  2:1768558/100‑1 (MQ=255)
                     ccGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAgagg                                                                               >  1:1797465/1‑100 (MQ=255)
                     ccGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAgagg                                                                               <  1:421028/100‑1 (MQ=255)
                     ccGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAgagg                                                                               >  1:3842092/1‑100 (MQ=255)
                     ccGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAgagg                                                                               <  1:1322673/100‑1 (MQ=255)
                      cGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAgagga                                                                              <  2:4206037/100‑1 (MQ=255)
                      cGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAgagg                                                                               >  1:1354322/1‑99 (MQ=255)
                        gTAGTTATACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCAGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGgaga                                                                            <  1:11074/100‑1 (MQ=255)
                        gTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGgaga                                                                            >  2:2533396/1‑100 (MQ=255)
                           cTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAAc                                                                          >  2:3471595/1‑99 (MQ=255)
                           cTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACg                                                                         >  2:3470851/1‑100 (MQ=255)
                             tttACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGcc                                                                       >  1:578342/1‑100 (MQ=255)
                               tccccGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCga                                                                     <  2:3799504/98‑1 (MQ=255)
                               tACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCga                                                                     >  2:2552834/1‑100 (MQ=255)
                                 cccGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCgaga                                                                   >  2:1080116/1‑100 (MQ=255)
                                 cccGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCgaga                                                                   >  1:2169849/1‑100 (MQ=255)
                                   cGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAAt                                                                 >  2:16754/1‑100 (MQ=255)
                                     gTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATgg                                                               >  2:2987394/1‑100 (MQ=255)
                                     gTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATgg                                                               <  1:290893/100‑1 (MQ=255)
                                      tGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGc                                                              <  1:3687046/100‑1 (MQ=255)
                                       gccgccTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGc                                                              >  1:1002531/1‑99 (MQ=255)
                                        ccgccTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCAc                                                            >  2:1489588/1‑100 (MQ=255)
                                        ccgccTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCAc                                                            <  2:1077234/100‑1 (MQ=255)
                                        ccgccTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACAGTCCTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCa                                                             <  1:2588345/99‑1 (MQ=255)
                                         cgccTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGCACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCAcc                                                           >  2:1278173/1‑100 (MQ=255)
                                         cgccTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCAcc                                                           >  2:1825496/1‑100 (MQ=255)
                                          gccTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCa                                                             <  1:3014397/97‑1 (MQ=255)
                                              ggCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCa                                                          >  1:2327672/1‑96 (MQ=255)
                                                ccTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGAc                                                    >  2:1957859/1‑100 (MQ=255)
                                                ccTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGAc                                                    <  2:2957387/100‑1 (MQ=255)
                                                ccTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGAc                                                    >  1:100506/1‑100 (MQ=255)
                                                 cTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGAcc                                                   <  2:167747/100‑1 (MQ=255)
                                                 cTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGAc                                                    <  1:1411448/99‑1 (MQ=255)
                                                  tAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGAcc                                                   >  1:2906497/1‑99 (MQ=255)
                                                   cttccgatcttggatggCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCAt                                                 <  2:3902476/90‑1 (MQ=255)
                                                   aaGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCa                                                  <  2:3293281/99‑1 (MQ=255)
                                                   aaGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCAt                                                 >  2:12539/1‑100 (MQ=255)
                                                    aGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATc                                                <  2:2423499/100‑1 (MQ=255)
                                                       tttCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATcagc                                             <  2:1706282/100‑1 (MQ=255)
                                                       tttCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATcagc                                             <  1:2987382/100‑1 (MQ=255)
                                                       tttCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATcagc                                             >  2:1822034/1‑100 (MQ=255)
                                                       tttCCATGGATGGCAGGCCGTTGTAACGTGCCACACGCGGCGAACCGTACTCCCAACGAGGCGGCGAGACCGACGAGACTGGCACCATCTCAcc                                                   >  1:3624035/1‑94 (MQ=255)
                                                        ttCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGCACGCCGAGAATGGCCCCATCTGACCATcagaa                                            >  1:1836946/1‑98 (MQ=255)
                                                         tCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATcagcag                                           >  2:2275574/1‑100 (MQ=255)
                                                          ccATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCagcagc                                          >  2:3057802/1‑100 (MQ=255)
                                                           cATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGACGAGGAGAACGCCGAGAATGGCACCATCTGACCATCCgcagca                                         >  1:709201/1‑100 (MQ=255)
                                                             tggatggCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCAtagatcggaag                                       >  1:3902458/1‑90 (MQ=255)
                                                             tggatggCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACg                                       <  2:3842109/100‑1 (MQ=255)
                                                              ggatggCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGa                                      >  2:2038443/1‑100 (MQ=255)
                                                                atggCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAAc                                    <  2:2844684/100‑1 (MQ=255)
                                                                 tggCAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACa                                   <  2:617075/100‑1 (MQ=255)
                                                                  ggCAGCCCGTTGTAACGTTCCCGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACAt                                  <  1:3843493/100‑1 (MQ=255)
                                                                    cAGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATAc                                <  1:1896924/100‑1 (MQ=255)
                                                                     aGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATAcc                               <  1:3077640/100‑1 (MQ=255)
                                                                     aGGCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATAcc                               <  2:631129/100‑1 (MQ=255)
                                                                       gCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAg                             <  1:4012383/100‑1 (MQ=255)
                                                                       gCCGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAg                             <  2:3951155/100‑1 (MQ=255)
                                                                        ccGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGt                            >  1:4232700/1‑100 (MQ=255)
                                                                        ccGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGt                            >  2:3248674/1‑100 (MQ=255)
                                                                         cGTTGTAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTc                           <  1:1080113/100‑1 (MQ=255)
                                                                         cGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTc                           <  2:2446902/100‑1 (MQ=255)
                                                                           ttGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGACGAGAACGCCGAGAATGCCACCACCTGCCCATCATCAGCACGAACAt                                  >  1:1742559/1‑91 (MQ=255)
                                                                             gTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTCGCCg                       >  2:1774200/1‑100 (MQ=255)
                                                                              tAACGTTCCAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTCGCCGa                      >  2:2227615/1‑100 (MQ=255)
                                                                              tAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGATCATACCAGTCGCCGa                      >  2:4213806/1‑100 (MQ=255)
                                                                               aaCGTTCCAGACGCGGCGAACCGTACTCCCAACGCGAAGAGGAGAACGCCGCGAATGGCACCCTCTGCCCATCAGCAGCACGCACATACCAGTc                           >  1:1483724/1‑94 (MQ=255)
                                                                               aaCGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTCGCCGat                     >  2:903837/1‑100 (MQ=255)
                                                                               aaCGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTCGCCGat                     >  1:2928738/1‑100 (MQ=255)
                                                                                 cGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTCGCCGatat                   >  2:3140288/1‑100 (MQ=255)
                                                                                    tCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTCGCCGATatcat                <  1:405627/100‑1 (MQ=255)
                                                                                      cAGACGCGGCGAGCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTCGCCGATATCATcc              <  2:1002534/100‑1 (MQ=255)
                                                                                      cAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTCGCCGATATCATcc              <  2:2018426/100‑1 (MQ=255)
                                                                                        gACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTCGCCGATATCATcc              >  1:1509497/1‑98 (MQ=255)
                                                                                                  aCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTCGCCGATATCATCCGGCAGCATACgg  >  2:3411070/1‑100 (MQ=255)
                                                                                                  aCCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTCGCCGATATCATCCGGCAGCATACg   >  1:239778/1‑99 (MQ=255)
                                                                                                  |                                                                                                   
TCCATCAGCTCCATTGCTTCACCGGTACTTTTACCCGGTGCCGCCTGGCCTAAGATTTCCATGGATGGCAGGCCGTTGTAACGTTCCAGACGCGGCGAACCGTACTCCCAACGAGAAGAGGAGAACGCCGAGAATGGCACCATCTGACCATCAGCAGCACGAACATACCAGTCGCCGATATCATCCGGCAGCATACGG  >  NC_000913/481870‑482067

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: