Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,248,432 (G)5→6 coding (660/663 nt) yqjA → general envelope maintenance protein; DedA family inner membrane protein; putative multidrug efflux transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,248,4271.G99.3% 531.4 / ‑2.8 143coding (655/663 nt)yqjAgeneral envelope maintenance protein; DedA family inner membrane protein; putative multidrug efflux transporter
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (85/57);  total (85/58)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 153 total aligned reads displayed.

GACCAGCTGATGTCATGCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGTTGCTGGCCTG  >  NC_000913/3248328‑3248521
                                                                                                    |                                                                                              
gACCAGCTGATGTCATGCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATc                                                                                                 <  1:235193/100‑1 (MQ=255)
 aCCAGCTGATGTCATGCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCg                                                                                               >  2:3016546/1‑100 (MQ=255)
    aacTGATGTCATGACTGATGCTGCTCCCGGTAGTGCTGCTGGGGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCgggg                                                                                            >  1:2610326/3‑100 (MQ=255)
    aGCTGATGTCATGCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCgggg                                                                                            <  1:1536499/100‑1 (MQ=255)
    aGCTGATGTCATGCCTGATGCAGCTCCCGGTGGTGCTGCTGGTGGTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCgggg                                                                                            <  1:178967/100‑1 (MQ=255)
     gCTGATGTCATGCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATAGGGAAATCggggg                                                                                           <  1:2535088/100‑1 (MQ=255)
           gTCATGCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGACTGGCAGGTTCTCTGGTCGTGTTATGAAAAAAGAAATATGGAAATCGGGGGGTAAgg                                                                                      <  2:112407/100‑1 (MQ=255)
            tCATGCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAggg                                                                                     <  1:1338336/100‑1 (MQ=255)
            tCATGCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAggg                                                                                     >  2:272544/1‑100 (MQ=255)
             cATGCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGa                                                                                    >  2:3786163/1‑100 (MQ=255)
             cATGCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGa                                                                                    >  2:119558/1‑100 (MQ=255)
              aTGCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTGTGGCCTGGCAGGTGCTCTGGTCGTGTTATTGAAAAAGAAATATGGAAATCGGGGGGTAAggg                                                                                     <  2:949291/98‑1 (MQ=255)
               tGCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATg                                                                                  <  2:2450754/100‑1 (MQ=255)
                gCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGc                                                                                 >  1:435569/1‑100 (MQ=255)
                 ccTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCa                                                                                >  2:2854496/1‑100 (MQ=255)
                 ccTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCa                                                                                <  1:74422/100‑1 (MQ=255)
                 ccTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCa                                                                                >  2:958990/1‑100 (MQ=255)
                   tGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCaaa                                                                              >  2:3841376/1‑100 (MQ=255)
                    gATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAAt                                                                             >  2:983025/1‑100 (MQ=255)
                     aTGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATa                                                                            <  1:3787276/100‑1 (MQ=255)
                       gctgctCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATAcc                                                                          <  1:3301606/100‑1 (MQ=255)
                         tgctCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTc                                                                        >  1:697404/1‑100 (MQ=255)
                           gtcacagtggtgCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGc                                                                      >  1:3689042/7‑100 (MQ=255)
                           ctCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATTCAAATACCTCGc                                                                      <  1:2547301/100‑1 (MQ=255)
                           ctCCCGGTGGTGATGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTCTGGCAAAAGAAATATGGCAATCGGGGGGTAAGGGATGCAAATACCTCGc                                                                      >  2:1735084/1‑100 (MQ=255)
                            tCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGACGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCa                                                                     <  2:2722366/100‑1 (MQ=255)
                              ccGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATg                                                                   >  2:569908/1‑100 (MQ=255)
                                ggtggtGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTc                                                                 >  2:1625310/1‑100 (MQ=255)
                                ggtggtGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTc                                                                 >  1:1180759/1‑100 (MQ=255)
                                 gtggtgCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCg                                                                <  1:3007967/100‑1 (MQ=255)
                                  tggtgCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGc                                                               <  1:133211/100‑1 (MQ=255)
                                     tgctgcTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCagatcggaag                                                            >  2:510681/1‑91 (MQ=255)
                                      gctgctGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCg                                                           >  2:3090416/1‑100 (MQ=255)
                                        tgctgGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGcc                                                         >  2:3502448/1‑100 (MQ=255)
                                        tgattgtgtTTGTACTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCg                                                           <  2:203129/93‑1 (MQ=255)
                                         gctgGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGcc                                                         >  2:1069102/1‑99 (MQ=255)
                                         gctgGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCa                                                        <  1:3500910/100‑1 (MQ=255)
                                         gctgGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCa                                                        >  1:711969/1‑100 (MQ=255)
                                         gctgGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCa                                                        <  2:1529153/100‑1 (MQ=255)
                                           tgGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGc                                                      <  1:2880109/100‑1 (MQ=255)
                                            ggTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCt                                                     >  1:2832689/1‑100 (MQ=255)
                                             gtgtTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTa                                                    <  1:857337/100‑1 (MQ=255)
                                              tgtTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAg                                                   >  1:526131/1‑100 (MQ=255)
                                               gtTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGc                                                  >  1:2219578/1‑100 (MQ=255)
                                               gtTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGc                                                  >  1:2473071/1‑100 (MQ=255)
                                               gtTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGc                                                  >  2:3313849/1‑100 (MQ=255)
                                               gtTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGc                                                  >  2:2224752/1‑100 (MQ=255)
                                                 ttGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTACCCt                                                <  1:694966/100‑1 (MQ=255)
                                                  tGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTg                                               >  2:2903845/1‑100 (MQ=255)
                                                   ggCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTgg                                              >  1:3159067/1‑100 (MQ=255)
                                                     ccTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTc                                            >  1:398238/1‑100 (MQ=255)
                                                     ccTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTc                                            >  2:1175131/1‑100 (MQ=255)
                                                      cTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTcc                                           >  1:2593644/1‑100 (MQ=255)
                                                      cTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTc                                            <  2:503955/99‑1 (MQ=255)
                                                       tGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCg                                          >  1:336670/1‑100 (MQ=255)
                                                         gCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGgc                                        >  1:2494532/1‑100 (MQ=255)
                                                         gCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGgc                                        <  1:2500545/100‑1 (MQ=255)
                                                           aGGTTCTCTGGGCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGgc                                        <  2:3291197/98‑1 (MQ=255)
                                                           aGGTGCTCTGGTCGTGTTATGGGAAAAGAAATATGGAACTCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGgc                                        <  2:3291059/98‑1 (MQ=255)
                                                            ggTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCt                                     >  1:2044052/1‑100 (MQ=255)
                                                             gTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTg                                    <  2:2189901/100‑1 (MQ=255)
                                                                ctctGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGttt                                 <  2:398237/100‑1 (MQ=255)
                                                                ctctGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGttt                                 <  1:2321013/100‑1 (MQ=255)
                                                                ctctGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGt                                   >  2:1101902/1‑98 (MQ=255)
                                                                 tctGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGtttt                                >  1:2019856/1‑100 (MQ=255)
                                                                  ctGGTCGTGTTGTGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTa                               >  1:3684211/1‑100 (MQ=255)
                                                                  ctGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTa                               >  2:749405/1‑100 (MQ=255)
                                                                    ggTCGTGTGATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACt                             >  2:1417208/1‑100 (MQ=255)
                                                                     gTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACtt                            >  1:2686863/1‑100 (MQ=255)
                                                                     gTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACtt                            <  2:2686848/100‑1 (MQ=255)
                                                                        gtgtTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTg                         <  1:1940649/100‑1 (MQ=255)
                                                                         tgtTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTgg                        >  1:1675325/1‑100 (MQ=255)
                                                                         tgtTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTgg                        >  2:2469424/1‑100 (MQ=255)
                                                                          gtTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGt                       >  2:236806/1‑100 (MQ=255)
                                                                          gtTATGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGt                       >  2:1760500/1‑100 (MQ=255)
                                                                             aTGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCgg                    <  2:2017658/100‑1 (MQ=255)
                                                                             aTGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCgg                    <  2:697403/100‑1 (MQ=255)
                                                                             aTGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCgg                    <  1:2392699/100‑1 (MQ=255)
                                                                              tGGAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGc                   >  1:2466188/1‑100 (MQ=255)
                                                                               ggAAAAAGAAATATGGCAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTTGTCGGCa                  <  2:2645577/100‑1 (MQ=255)
                                                                               ggAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCg                     >  2:3394564/1‑97 (MQ=255)
                                                                               ggAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCa                  >  1:2197574/1‑100 (MQ=255)
                                                                               ggAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCa                  >  2:1504790/1‑100 (MQ=255)
                                                                                gAAAAAGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCAc                 >  1:548770/1‑100 (MQ=255)
                                                                                  aaaaGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGc               >  2:3038337/1‑100 (MQ=255)
                                                                                   aaaGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCt              <  1:1488114/100‑1 (MQ=255)
                                                                                     aGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGt            >  2:3404620/1‑100 (MQ=255)
                                                                                     aGAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGt            >  2:1040001/1‑100 (MQ=255)
                                                                                      gAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGtt           <  1:1043748/100‑1 (MQ=255)
                                                                                      gAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGtt           >  1:2881482/1‑100 (MQ=255)
                                                                                      gAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGtt           <  1:2903860/100‑1 (MQ=255)
                                                                                      gAAATATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGtt           >  2:2551931/1‑100 (MQ=255)
                                                                                       aaaTATGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGTTg          <  1:3175083/100‑1 (MQ=255)
                                                                                          tatGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGTTGCTg       <  1:3797952/100‑1 (MQ=255)
                                                                                           atGGAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCTAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGTTGCTgg      >  2:255954/1‑100 (MQ=255)
                                                                                             ggAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGTTGCTGGcc    >  2:589754/1‑100 (MQ=255)
                                                                                              gAAATCGGGGGGTAAGGGATGCAAATATCTCGCATGTCTCGTCGCGAGCGAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGTTGCTGGCCt   >  2:318717/1‑100 (MQ=255)
                                                                                              gAAATCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGTTGCTGGCCt   >  1:898333/1‑100 (MQ=255)
                                                                                               aaaTCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGTTGCTGGCCTg  >  1:1805568/1‑100 (MQ=255)
                                                                                               aaaTCGGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGTTGCTGGCCTg  >  2:1057190/1‑100 (MQ=255)
                                                                                                    |                                                                                              
GACCAGCTGATGTCATGCCTGATGCTGCTCCCGGTGGTGCTGCTGGTGTTTGGCCTGGCAGGTTCTCTGGTCGTGTTATGGAAAAAGAAATATGGAAATCGGGGGTAAGGGATGCAAATACCTCGCATGTCGCTTCGCCAGCTAGCCTGGTCCGGCGCTGTTTTACTTCTGGTCGGCACGCTGTTGCTGGCCTG  >  NC_000913/3248328‑3248521

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: