Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 483,817 A→G F196S (TTC→TCC)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913483,8170AG98.6% 261.4 / ‑3.7 71F196S (TTC→TCC) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (37/33);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

AGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATATCCTCCTGCGTCATGGTGC  >  NC_000913/483675‑483964
                                                                                                                                               |                                                                                                                                                   
aGTAGGGGTGAGCGGCGTCTGAACAATAAGGGGGGTTTTAGGTGGTGGTCCTTTCCCCGGGGGCGCCCCCCCGATTGGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTg                                                                                                                                                <  2:832713/148‑1 (MQ=255)
aGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGtt                                                                                                                                              >  2:1280205/1‑150 (MQ=255)
 gtgggggTCAGCCGCGTCGGAGAAAAAAGGGAGGCGTTAGGCTGTTGCCCTTTCCCCGGCGGCGCCCCCCAGAGCGGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCa                                                                                                                                            <  2:23408/146‑1 (MQ=255)
   agagGTAAGCCGCGTCGGACAAATAGTAGGGGTGTTAGTGTGTTGCCCTTTCGCCGCGGGCGTCCCCCGGGGCTTACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTct                                                                                                                                            <  2:1068594/149‑2 (MQ=255)
    ggggTCAGACGCGTCGGAGAAATAATAGGGGTGTTAGGTTGTTGGCCTTCCCCCGGGGGCGCCCCCCGGAGCTGCCCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAg                                                                                                                                           <  2:808838/147‑1 (MQ=255)
    gagGTCAGCCGCGTCTAAAAAAAAAGGGAGGCTTAAGGTTTTTGCCCTTTCCCCGGCGGCGTCCCCCCGGGCTCCCCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCa                                                                                                                                            <  2:900482/148‑1 (MQ=255)
     agGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTc                                                                                                                                        >  2:143125/1‑151 (MQ=255)
            aCGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCggg                                                                                                                                 >  2:1269603/1‑151 (MQ=255)
              gcgTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGtt                                                                                                                               >  1:610002/1‑151 (MQ=255)
               cgTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGCTTTGTTCAGCTCATTCGGGTTc                                                                                                                              >  1:1157278/1‑151 (MQ=255)
               cgTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTc                                                                                                                              >  1:634133/1‑151 (MQ=255)
                gTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTc                                                                                                                              >  1:909701/1‑150 (MQ=255)
                 tCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCAt                                                                                                                            >  2:562890/1‑151 (MQ=255)
                    gAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCAt                                                                                                                            >  1:337197/1‑148 (MQ=255)
                        aataaaAGAGGGGTTAGGCTGTTGCCCTTCCCCCGGCGGCTTCCCACCGGTCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGa                                                                                                                       <  2:1228275/149‑1 (MQ=255)
                           aataGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGc                                                                                                                  >  1:19844/1‑151 (MQ=255)
                            ataGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGc                                                                                                                  >  1:1332531/1‑150 (MQ=255)
                                 ggCGTTAAGGTGTTGGCTTTCCCCGGGCGGCGTCCCCCGGGGCTACCCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATcgcg                                                                                                            <  2:1218809/151‑1 (MQ=255)
                                      taaggctggtgCCTTTCCCCCGGGGGCGTCCCCCCGGGCGGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGt                                                                                                           <  2:1275031/138‑1 (MQ=255)
                                          cTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCCGGTTCATCCAGATACGCATCGCGTACtgtg                                                                                                     >  2:1237513/1‑149 (MQ=255)
                                           tGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGaa                                                                                                   >  1:37026/1‑150 (MQ=255)
                                            gTTGCCCTTCCCCCGGCGGCGTCCCCCCGAGCTGCCCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGa                                                                                                    <  2:263618/148‑1 (MQ=255)
                                             ttGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCg                                                                                                >  1:1070621/1‑151 (MQ=255)
                                                   tttCCCGGGGGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTAc                                                                                                         <  2:1094477/136‑1 (MQ=255)
                                                   ttccACCGGCGCCGTCCCCCCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAAcc                                                                                                 <  2:582618/141‑1 (MQ=255)
                                                    ttcccccGGGGGCGTACCCCCGGGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGa                                                                                                    <  1:150332/136‑1 (MQ=255)
                                                      cccccGGGGGCGCCCCCCCGAGCGGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGaacaac                                                                                          <  1:164396/146‑1 (MQ=255)
                                                          cggcggCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAAc                                                                                     >  1:749248/1‑150 (MQ=255)
                                                           ggcggcGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTTAGCTCATTCGGGTTCATCCAGATACGCCTCGCGTACTGTGAACCGAACAACTGAACAt                                                                                   >  1:580581/1‑151 (MQ=255)
                                                           ggcggcCTCCCCCCGCGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGaaca                                                                                            <  2:1382149/142‑1 (MQ=255)
                                                                    ccaccGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGc                                                                          >  2:835632/1‑151 (MQ=255)
                                                                      accGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGccc                                                                        >  2:850307/1‑151 (MQ=255)
                                                                          aGCTGACCCGCCGCAACCTGGGGGTTCTGCGCTTTGATGGCGGTAATTACATCAACCGGCGTTATCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACg                                                                    >  2:59972/1‑151 (MQ=255)
                                                                                cccGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTa                                                                  <  2:304958/147‑1 (MQ=255)
                                                                                 ccgccgCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCt                                                             <  1:686555/151‑1 (MQ=255)
                                                                                      gcAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATgg                                                        >  2:507679/1‑151 (MQ=255)
                                                                                      gcAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATgg                                                        <  2:440036/151‑1 (MQ=255)
                                                                                             gggCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGCCGTACGGCTGACGGCctctc                                                  <  2:609929/150‑2 (MQ=255)
                                                                                             gggCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCtt                                                  <  1:12246/150‑1 (MQ=255)
                                                                                               gCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCa                                               <  2:189651/151‑1 (MQ=255)
                                                                                                cGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCat                                              >  2:47139/1‑151 (MQ=255)
                                                                                                cGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCat                                              <  1:1104321/151‑1 (MQ=255)
                                                                                                 gTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCata                                             <  2:580512/151‑1 (MQ=255)
                                                                                                     tGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTc                                          >  1:743531/1‑150 (MQ=255)
                                                                                                        gcTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTcgccg                                      <  1:73189/151‑1 (MQ=255)
                                                                                                        gcTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTcgcc                                       >  1:693299/1‑150 (MQ=255)
                                                                                                        gcTTTGATGGCGGGAATGACATCCACCGGCGTTATCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGCACAACTGAACATCACCCACGCCCGCCGTACGGCTGATGGCATCTTTCATATTcgccg                                      >  2:1324015/1‑151 (MQ=255)
                                                                                                          tttGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCgccgcc                                    >  1:831104/1‑151 (MQ=255)
                                                                                                            tGATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCgccgcc                                    <  1:1237613/149‑1 (MQ=255)
                                                                                                             gATGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACg                                 <  2:749171/151‑1 (MQ=255)
                                                                                                              aTGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGt                                <  1:437776/151‑1 (MQ=255)
                                                                                                               tGGCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGt                                >  1:623789/1‑150 (MQ=255)
                                                                                                                ggCGGTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTa                               <  1:143145/150‑1 (MQ=255)
                                                                                                                    gTAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCgg                          <  2:858030/151‑1 (MQ=255)
                                                                                                                     tAATGACATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCgg                          >  1:1448367/1‑150 (MQ=255)
                                                                                                                          aCATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGatat                    >  1:941524/1‑151 (MQ=255)
                                                                                                                          aCATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGatat                    >  2:631863/1‑151 (MQ=255)
                                                                                                                           cATCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAAGAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATAtc                   >  1:280770/1‑151 (MQ=255)
                                                                                                                            aTCAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATAtcc                  <  1:1375808/151‑1 (MQ=255)
                                                                                                                              cAACCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATAtcctc                >  2:209071/1‑151 (MQ=255)
                                                                                                                               aaCCGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATAtcctcc               >  2:155660/1‑151 (MQ=255)
                                                                                                                                  cGGCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATATCCTCCTGc            >  1:1289652/1‑151 (MQ=255)
                                                                                                                                  cGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATATCCTCCTGc            >  2:993898/1‑151 (MQ=255)
                                                                                                                                    gCGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATATCCTCCTGCGt          <  2:789948/151‑1 (MQ=255)
                                                                                                                                     cGTTAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATATCCTCCTGCGTc         <  1:1269705/151‑1 (MQ=255)
                                                                                                                                        tAGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATATCCTCCTGCGTCATg      >  2:1236474/1‑151 (MQ=255)
                                                                                                                                         aGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATATCCTCCTGCGTCATgg     <  1:507737/151‑1 (MQ=255)
                                                                                                                                         aGCTGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATATCCTCCTGCGTCATg      <  2:683644/150‑1 (MQ=255)
                                                                                                                                            tGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATATCCTCCTGCGTCATGGTGc  >  1:587667/1‑151 (MQ=255)
                                                                                                                                            tGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATATCCTCCTGCGTCATGGTGc  >  2:282735/1‑151 (MQ=255)
                                                                                                                                            tGGGATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATATCCTCCTGCGTCATGGTGc  >  2:1400801/1‑151 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
AGTAGAGGTCAGACGCGTCTGAGCAATAATAGAGGCGTTAAGCTGTTGGCCTTTCACCGGCGGCGTACCACCGAGCTGACCCGCCGCAACCTGGGCGTTCTGCGCTTTGATGGCGGTAATGACATCAACCGGCGTTAGCTGGAATTTGTTCAGCTCATTCGGGTTCATCCAGATACGCATCGCGTACTGTGAACCGAACAACTGAACATCACCCACGCCCGACGTACGGCTGATGGCATCTTTCATATTCGCCGCCACGTAGTCGGAGATATCCTCCTGCGTCATGGTGC  >  NC_000913/483675‑483964

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