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

GCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCC  >  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)
                         agagAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑c                                                                                                                             >  2:114573/1‑151 (MQ=255)
                             aaTGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAACCGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑Cccgc                                                                                                                         <  2:129864/151‑1 (MQ=255)
                             aaTGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑Cccgc                                                                                                                         >  1:881568/1‑151 (MQ=255)
                             aaTGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑Cccgc                                                                                                                         >  1:827995/1‑151 (MQ=255)
                                      gCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCg                                                                                                                >  1:983313/1‑151 (MQ=255)
                                        gAGGATAAATCCAATAATCACCTTTAAACCCCCCTTCCATTCCATCCTGGCTTAGTTTAAAATCAGCTTAAACTCCCCAACCCCGATTATAGACCCCCAAATCACATATTTCCCCCCAATTTTATCCCCTACACC‑CCCTCCTATCAcctt                                                                                                               <  2:102370/150‑3 (MQ=255)
                                        gAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGt                                                                                                               >  2:424577/1‑150 (MQ=255)
                                        gAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTa                                                                                                              <  1:644566/151‑1 (MQ=255)
                                            aTAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTg                                                                                                          <  1:564777/151‑1 (MQ=255)
                                             tAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGc                                                                                                         >  2:522542/1‑151 (MQ=255)
                                                aTCCAATAATCACCGTTAAACCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGc                                                                                                      <  1:788190/151‑1 (MQ=255)
                                                 tCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGc                                                                                                      >  1:679507/1‑150 (MQ=255)
                                                    aataatCACCGTAAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACg                                                                                                  <  1:370960/151‑1 (MQ=255)
                                                     ataatCACCGTAAAACCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTCCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACgg                                                                                                 <  1:751916/151‑1 (MQ=255)
                                                      taatCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGc                                                                                                >  1:234719/1‑151 (MQ=255)
                                                        atCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCa                                                                                               >  1:45713/1‑150 (MQ=255)
                                                            ccGTAAAACCCCCCGCCCAGGGCGGGCTGGCGTAGGGAAAAATCAGCGTAACCGGCGCAACGCCGAGGATAGCCCCCAAAATCACATATTCCCCCCCAATTTTACCCCCTACAGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCaaa                                                                                             <  2:329157/148‑1 (MQ=255)
                                                              gTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAgg                                                                                         >  1:971269/1‑150 (MQ=255)
                                                                 aacccccccGTCAAGTGCAGGCTGGCGTAGGGAAAAATCAGCGTAAGCGGCGCAACCCCGAGGATAGGCCCCCAAATCCCATTTTCCCCCCCATTTTTACCCCCTACAGG‑CCCGCAGATCCCCGTACCTGCCGCAACGGCAAACAGGAAgg                                                                                     <  2:114436/147‑1 (MQ=255)
                                                                  aaccccccGCCCAGGGCGGGCGGGCGTGGGGAAAAACCAGCTAAACCGGCGCACCGCGGAGGATAGACCCCCAAACCACATTTTCCCCCCCAATTTTACCCCCTACGGC‑CCCGCCGATCCCCGTACCTGCCGCAACGGCAAACAGGAAGGc                                                                                    <  2:129458/148‑1 (MQ=255)
                                                                   aaCCCCCGTCCAGTGCAAGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGc                                                                                    >  2:779150/1‑150 (MQ=255)
                                                                    aCCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCa                                                                                   >  1:968244/1‑150 (MQ=255)
                                                                     cccccGTCCAGTGCAGGCTGGCGTGGGGTAAAACCAGCGAAACCGGCGCACCGCCGAGGATAGCGCCCAAAACCACATTTTCCCCCCCAATTTTACCCCCTACGGC‑CCCGCCGATAACCGTCCCTGCCGCAACGGCAAACAGGAAGGCaa                                                                                  <  2:731315/150‑1 (MQ=255)
                                                                          gTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCAAATATTCCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGAt                                                                            <  1:724901/151‑1 (MQ=255)
                                                                             cAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAGCGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGa                                                                          >  2:377646/1‑150 (MQ=255)
                                                                                  cAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTg                                                                     >  1:687401/1‑150 (MQ=255)
                                                                                    ggCTGGCGAAGGAAAAAATCAGCGAAACCGGCGCACCGCCGAGAATAGACCCCAAAACCACATATTCCCCCCCATTTTTACCCCCAACGGG‑CCCGCGGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAg                                                                   <  2:62197/150‑1 (MQ=255)
                                                                                       tGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCAtt                                                               >  2:302343/1‑151 (MQ=255)
                                                                                        ggCGTAGGAAAAAACCAGCGAAAACGCCGCACCCCCGAGGATAGACCCCAAAACCACATTTTCCCCCCCAATTTTACCCCCTCCGGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCAt                                                                <  2:311503/149‑1 (MQ=255)
                                                                                              gggTAAAATCAGCTTAAACGGCCCACCGCCGAGGATAGAGCCCCAAATCACATATTCCCCCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGga                                                          <  1:889331/149‑1 (MQ=255)
                                                                                              ggatAAAACCAGCGAAAACGGCGAAACGCCGAGAATAGGGCCCAAAATCACATATTCCCCCCAAATTTTACCCCCAAAAGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGga                                                          <  2:187726/146‑1 (MQ=255)
                                                                                                  aaaaCCAGCTTAAACGGCGCACCGCCGAGGATAGCGCCCAAAATCACATTTTCCCCCCCAATTTTACCCCCAACAGG‑CCCGCCAATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGgatagata                                                    <  2:810198/151‑1 (MQ=255)
                                                                                                    aaTCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAAt                                                  >  2:340548/1‑151 (MQ=255)
                                                                                                        aGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCg                                               >  2:391645/1‑150 (MQ=255)
                                                                                                         gCGTAACCGGCGCAACGCCGAGGATAGAGCCCCAAATCCCATATTCCCCCCAAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGaa                                             <  1:404416/151‑1 (MQ=255)
                                                                                                         gCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGa                                              >  1:869796/1‑150 (MQ=255)
                                                                                                          cgaaaaCGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTCCCGCCCAATTTTACCCCCTACAGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAt                                            <  1:168648/148‑1 (MQ=255)
                                                                                                             aaaCGGCGCACCGCCGAGGATAGAGCCCCAAATCACATATTCCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAt                                            <  1:338852/148‑1 (MQ=255)
                                                                                                             aaaCGGCGCAACGCCGAGGATAGACCCCAAAATCACATATTCCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAtt                                           <  1:165445/149‑1 (MQ=255)
                                                                                                             aaaCGGCCCACCGCCGAGGATAGACCCCCAAATCACATATTCCCCCCAAATTTTACCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAt                                            <  1:843116/148‑1 (MQ=255)
                                                                                                              aaCGGCGCAACGCGGAGGATAGAGCCCAAAATCACATATTCCCCCCCATTTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAAttttt                                        <  1:96426/151‑1 (MQ=255)
                                                                                                                cGGCCCACCGCCGAGGAAAGACCCCAAAACCAAATATTCCCGCCCAATTCAACCCCCTACAGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTg                                       <  2:723988/150‑1 (MQ=255)
                                                                                                                 ggCGCACCGCCGAGGATAGAGCCCCAAACCACATATTCCCCCCAAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGc                                      <  1:621207/150‑1 (MQ=255)
                                                                                                                   cgcACCGCCGAGGATAGCCCCCAAAATCACATATTCCCCCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCa                                     <  1:505835/149‑1 (MQ=255)
                                                                                                                    gcAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATc                                   >  1:366371/1‑150 (MQ=255)
                                                                                                                         gCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATa                              >  1:800435/1‑150 (MQ=255)
                                                                                                                          ccGAGGATAGAGCCCCAAACCACATATTCCCCCCCAATTTTATCCCCAACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATaa                             <  1:238255/150‑1 (MQ=255)
                                                                                                                            gagaATAGACCCCAAAACCAATTATTCCCGCCAAATTTTA‑CCCCTTCTGGCCCCGCCTATCCCCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGAt                               <  1:145599/146‑1 (MQ=255)
                                                                                                                              ggaaggaccccaAAATCCCATATTCCCCCCAAATTTTTCCCCCAACGGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATaaa                            <  1:265944/139‑1 (MQ=255)
                                                                                                                              ggATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAgg                         >  1:37726/1‑150 (MQ=255)
                                                                                                                              ggAAAGACCCCAAAATCACATATTCCCCCCCAATTTTATCCCCTACGGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAgg                         <  2:988934/150‑1 (MQ=255)
                                                                                                                                   gagCCCAAAATCCCATATTCCCGCCCAATTTTACCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAAt                   <  1:699586/151‑1 (MQ=255)
                                                                                                                                    agCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAAt                   >  2:30299/1‑150 (MQ=255)
                                                                                                                                     cccccAAACCACATATTCCCGCCCAATTTTACCCCCTACGGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTa                     <  2:958400/146‑1 (MQ=255)
                                                                                                                                      ccccAAACCACATTTTTCCGCCCAATTTCACCCCCTACAGC‑CCCGCCGATCTCCGTCCCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAgg                         <  2:14317/142‑1 (MQ=255)
                                                                                                                                      ccccAAACCACATTTTCCCGCCCAATTTTACCCCCAACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGgtg                       <  2:915237/144‑1 (MQ=255)
                                                                                                                                         cAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAgctgc             >  2:409612/1‑151 (MQ=255)
                                                                                                                                          aaaTCACATATTCCCCCCCAATTTTATCCCCTACGGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAgctgct            <  1:289005/151‑1 (MQ=255)
                                                                                                                                          aaaCCAAATATTCCCCCCCAATTTTATCCCCTAAAGC‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAgctg              <  2:117647/149‑1 (MQ=255)
                                                                                                                                              cacaTATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATg        >  1:912982/1‑151 (MQ=255)
                                                                                                                                               acaTATTCCCGCCAAATTTTACCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGAt         <  1:162397/149‑1 (MQ=255)
                                                                                                                                                caTATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGc       >  2:938529/1‑150 (MQ=255)
                                                                                                                                                 atatTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTc     >  1:967160/1‑151 (MQ=255)
                                                                                                                                                    tttCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGcc  >  1:991469/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GCTCCTGAGCGTAGACCAGAATGGCAGAGAATGCCGAAGCGAGGATAAATCCAATAATCACCGTTAAAACCCCCGTCCAGTGCAGGCTGGCGTAGGGTAAAATCAGCGTAAACGGCGCAACGCCGAGGATAGAGCCCCAAATCACATATTTCCGCCCAATTTTATCCCCTACAGG‑CCCGCCGATCACCGTACCTGCCGCAACGGCAAACAGGAAGGCAAACAGATGAAGCTGAGCATTCTGGATAGATAATCCGAATTTTTGCATCAGATAAAAGGTGTAATAGCTGCTGATGCTCGCC  >  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: ‑