Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913503,8200TC43.5% 93.0 / 24.2 69K293E (AAA→GAA) fsrputative fosmidomycin efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (27/12);  new base C (0/30);  total (27/42)
Fisher's exact test for biased strand distribution p-value = 3.59e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCC  >  NC_000913/503670‑503968
                                                                                                                                                      |                                                                                                                                                     
gcTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATAttt                                                                                                                                                       >  1:474012/1‑151 (MQ=255)
   ccTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTc                                                                                                                                                      <  1:182139/149‑1 (MQ=255)
          gTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAAtt                                                                                                                                             <  1:50242/151‑1 (MQ=255)
                cAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATcc                                                                                                                                       >  1:875588/1‑151 (MQ=255)
                   aaTGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCt                                                                                                                                    <  2:73798/151‑1 (MQ=255)
                   aaTGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCt                                                                                                                                    <  1:352895/151‑1 (MQ=255)
                        cAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAgg                                                                                                                               <  2:846441/151‑1 (MQ=255)
                         agagAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGc                                                                                                                             >  2:114573/1‑151 (MQ=255)
                             aaTGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAACCGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCccgc                                                                                                                         <  2:129864/151‑1 (MQ=255)
                             aaTGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCccgc                                                                                                                         >  1:881568/1‑151 (MQ=255)
                             aaTGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCccgc                                                                                                                         >  1:827995/1‑151 (MQ=255)
                                      gCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCg                                                                                                                >  1:983313/1‑151 (MQ=255)
                                        gAGGATAAATCCAATAATCACCTTTAAACCCCCCTTCCATTCCATCCTGGCTTAGTTTAAAATCAGCTTAAACTCCCCAACCCCGATTATAGACCCCCAAATCACATATTTCCCCCCAATTTTATCCCCTACACCCCCTCCTATCAcctt                                                                                                               <  2:102370/150‑3 (MQ=255)
                                        gAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGt                                                                                                               >  2:424577/1‑150 (MQ=255)
                                        gAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTa                                                                                                              <  1:644566/151‑1 (MQ=255)
                                            aTAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTg                                                                                                          <  1:564777/151‑1 (MQ=255)
                                             tAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGc                                                                                                         >  2:522542/1‑151 (MQ=255)
                                                aTCCAATAATCACCGTTAAACCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGc                                                                                                      <  1:788190/151‑1 (MQ=255)
                                                 tCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGc                                                                                                      >  1:679507/1‑150 (MQ=255)
                                                    aataatCACCGTAAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACg                                                                                                  <  1:370960/151‑1 (MQ=255)
                                                     ataatCACCGTAAAACCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTCCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACgg                                                                                                 <  1:751916/151‑1 (MQ=255)
                                                      taatCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGc                                                                                                >  1:234719/1‑151 (MQ=255)
                                                        atCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCa                                                                                               >  1:45713/1‑150 (MQ=255)
                                                            ccGTAAAACCCCCCGCCCAGGGCGGGCTGGCGTAGGGAAAAATCAGCGTAACCGGCGCAACGCCGAGGATAGCCCCCAAAATCACATATTCCCCCCCAATTTTACCCCCTACAGCCCCGCCGATCACCGTACCTGCCGCAACGGCaaa                                                                                             <  2:329157/148‑1 (MQ=255)
                                                              gTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAgg                                                                                         >  1:971269/1‑150 (MQ=255)
                                                                 aacccccccGTCAAGTGCAGGCTGGCGTAGGGAAAAATCAGCGTAAGCGGCGCAACCCCGAGGATAGGCCCCCAAATCCCATTTTCCCCCCCATTTTTACCCCCTACAGGCCCGCAGATCCCCGTACCTGCCGCAACGGCAAACAGGAAgg                                                                                     <  2:114436/147‑1 (MQ=255)
                                                                  aaccccccGCCCAGGGCGGGCGGGCGTGGGGAAAAACCAGCTAAACCGGCGCACCGCGGAGGATAGACCCCCAAACCACATTTTCCCCCCCAATTTTACCCCCTACGGCCCCGCCGATCCCCGTACCTGCCGCAACGGCAAACAGGAAGGc                                                                                    <  2:129458/148‑1 (MQ=255)
                                                                   aaCCCCCGTCCAGTGCAAGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGc                                                                                    >  2:779150/1‑150 (MQ=255)
                                                                    aCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCa                                                                                   >  1:968244/1‑150 (MQ=255)
                                                                     cccccGTCCAGTGCAGGCTGGCGTGGGGTAAAACCAGCGAAACCGGCGCACCGCCGAGGATAGCGCCCAAAACCACATTTTCCCCCCCAATTTTACCCCCTACGGCCCCGCCGATAACCGTCCCTGCCGCAACGGCAAACAGGAAGGCaa                                                                                  <  2:731315/150‑1 (MQ=255)
                                                                          gTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCAAATATTCCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGAt                                                                            <  1:724901/151‑1 (MQ=255)
                                                                             cAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAGCGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGa                                                                          >  2:377646/1‑150 (MQ=255)
                                                                                  cAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTg                                                                     >  1:687401/1‑150 (MQ=255)
                                                                                    ggCTGGCGAAGGAAAAAATCAGCGAAACCGGCGCACCGCCGAGAATAGACCCCAAAACCACATATTCCCCCCCATTTTTACCCCCAACGGGCCCGCGGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAg                                                                   <  2:62197/150‑1 (MQ=255)
                                                                                       tGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCAtt                                                               >  2:302343/1‑151 (MQ=255)
                                                                                        ggCGTAGGAAAAAACCAGCGAAAACGCCGCACCCCCGAGGATAGACCCCAAAACCACATTTTCCCCCCCAATTTTACCCCCTCCGGCCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCAt                                                                <  2:311503/149‑1 (MQ=255)
                                                                                              gggTAAAATCAGCTTAAACGGCCCACCGCCGAGGATAGAGCCCCAAATCACATATTCCCCCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGga                                                          <  1:889331/149‑1 (MQ=255)
                                                                                              ggatAAAACCAGCGAAAACGGCGAAACGCCGAGAATAGGGCCCAAAATCACATATTCCCCCCAAATTTTACCCCCAAAAGCCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGga                                                          <  2:187726/146‑1 (MQ=255)
                                                                                                  aaaaCCAGCTTAAACGGCGCACCGCCGAGGATAGCGCCCAAAATCACATTTTCCCCCCCAATTTTACCCCCAACAGGCCCGCCAATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGgatagata                                                    <  2:810198/151‑1 (MQ=255)
                                                                                                    aaTCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAAt                                                  >  2:340548/1‑151 (MQ=255)
                                                                                                        aGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCg                                               >  2:391645/1‑150 (MQ=255)
                                                                                                         gCGTAACCGGCGCAACGCCGAGGATAGAGCCCCAAATCCCATATTCCCCCCAAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGaa                                             <  1:404416/151‑1 (MQ=255)
                                                                                                         gCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGa                                              >  1:869796/1‑150 (MQ=255)
                                                                                                          cgaaaaCGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTCCCGCCCAATTTTACCCCCTACAGCCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAt                                            <  1:168648/148‑1 (MQ=255)
                                                                                                             aaaCGGCGCACCGCCGAGGATAGAGCCCCAAATCACATATTCCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAt                                            <  1:338852/148‑1 (MQ=255)
                                                                                                             aaaCGGCGCAACGCCGAGGATAGACCCCAAAATCACATATTCCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAtt                                           <  1:165445/149‑1 (MQ=255)
                                                                                                             aaaCGGCCCACCGCCGAGGATAGACCCCCAAATCACATATTCCCCCCAAATTTTACCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAt                                            <  1:843116/148‑1 (MQ=255)
                                                                                                              aaCGGCGCAACGCGGAGGATAGAGCCCAAAATCACATATTCCCCCCCATTTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAttttt                                        <  1:96426/151‑1 (MQ=255)
                                                                                                                cGGCCCACCGCCGAGGAAAGACCCCAAAACCAAATATTCCCGCCCAATTCAACCCCCTACAGCCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTg                                       <  2:723988/150‑1 (MQ=255)
                                                                                                                 ggCGCACCGCCGAGGATAGAGCCCCAAACCACATATTCCCCCCAAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGc                                      <  1:621207/150‑1 (MQ=255)
                                                                                                                   cgcACCGCCGAGGATAGCCCCCAAAATCACATATTCCCCCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCa                                     <  1:505835/149‑1 (MQ=255)
                                                                                                                    gcAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATc                                   >  1:366371/1‑150 (MQ=255)
                                                                                                                         gCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATa                              >  1:800435/1‑150 (MQ=255)
                                                                                                                          ccGAGGATAGAGCCCCAAACCACATATTCCCCCCCAATTTTATCCCCAACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATaa                             <  1:238255/150‑1 (MQ=255)
                                                                                                                            gagaATAGACCCCAAAACCAATTATTCCCGCCAAATTTTACCCCTTCTGGCCCCGCCTATCCCCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGAt                               <  1:145599/146‑1 (MQ=255)
                                                                                                                              ggaaggaccccaAAATCCCATATTCCCCCCAAATTTTTCCCCCAACGGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATaaa                            <  1:265944/139‑1 (MQ=255)
                                                                                                                              ggATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAgg                         >  1:37726/1‑150 (MQ=255)
                                                                                                                              ggAAAGACCCCAAAATCACATATTCCCCCCCAATTTTATCCCCTACGGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAgg                         <  2:988934/150‑1 (MQ=255)
                                                                                                                                   gagCCCAAAATCCCATATTCCCGCCCAATTTTACCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAAt                   <  1:699586/151‑1 (MQ=255)
                                                                                                                                    agCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAAt                   >  2:30299/1‑150 (MQ=255)
                                                                                                                                     cccccAAACCACATATTCCCGCCCAATTTTACCCCCTACGGCCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTa                     <  2:958400/146‑1 (MQ=255)
                                                                                                                                      ccccAAACCACATTTTTCCGCCCAATTTCACCCCCTACAGCCCCGCCGATCTCCGTCCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAgg                         <  2:14317/142‑1 (MQ=255)
                                                                                                                                      ccccAAACCACATTTTCCCGCCCAATTTTACCCCCAACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGgtg                       <  2:915237/144‑1 (MQ=255)
                                                                                                                                         cAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAgctgc             >  2:409612/1‑151 (MQ=255)
                                                                                                                                          aaaTCACATATTCCCCCCCAATTTTATCCCCTACGGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAgctgct            <  1:289005/151‑1 (MQ=255)
                                                                                                                                          aaaCCAAATATTCCCCCCCAATTTTATCCCCTAAAGCCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAgctg              <  2:117647/149‑1 (MQ=255)
                                                                                                                                              cacaTATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATg        >  1:912982/1‑151 (MQ=255)
                                                                                                                                               acaTATTCCCGCCAAATTTTACCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGAt         <  1:162397/149‑1 (MQ=255)
                                                                                                                                                caTATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGc       >  2:938529/1‑150 (MQ=255)
                                                                                                                                                 atatTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTc     >  1:967160/1‑151 (MQ=255)
                                                                                                                                                    tttCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGcc  >  1:991469/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGGCCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCC  >  NC_000913/503670‑503968

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