Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,2010TG40.9% 108.4 / 19.1 66T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (7/32);  new base G (27/0);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 2.28e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  NC_000913/1356054‑1356351
                                                                                                                                                      |                                                                                                                                                          
cGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGGGtt                                                                                                                                                         >  1:90474/1‑151 (MQ=255)
    gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTaaa                                                                                                                                                      <  1:924316/150‑1 (MQ=255)
         aGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCa                                                                                                                                                 <  2:19833/150‑1 (MQ=255)
            cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATa                                                                                                                                             <  1:293473/151‑1 (MQ=255)
             ccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAg                                                                                                                                            <  1:282973/151‑1 (MQ=255)
              cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑TGC‑GGGGTTtaaagac                                                                                                                                                  >  2:533783/1‑140 (MQ=255)
              cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTtaaattcaaaaagg                                                                                                                                           >  1:887087/1‑140 (MQ=255)
              cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGCACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑TTC‑GGTGTTTAATACCAGaaa                                                                                                                                             >  1:1357983/1‑147 (MQ=255)
                 gAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGGGGTAAATGCCAGAAAGGCg                                                                                                                                         >  1:856402/1‑150 (MQ=255)
                  aGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGG‑GGGGGTTAATGCCAGATAAGCg                                                                                                                                         >  1:614214/1‑149 (MQ=255)
                     cgcgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑Gtt                                                                                                                                   <  1:352833/151‑1 (MQ=255)
                       cgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTc                                                                                                                                  <  1:1278589/150‑1 (MQ=255)
                           gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATccc                                                                                                                             <  1:475139/151‑1 (MQ=255)
                              aaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGG‑GGGGGTTAATACCCGAGAGGCGAC‑GGTTATCCCAgg                                                                                                                          >  1:437205/1‑151 (MQ=255)
                                   gagGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCa                                                                                                                     <  1:1311413/151‑1 (MQ=255)
                                     ggCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑g                                                                                                                   <  2:846703/150‑1 (MQ=255)
                                           cgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTT‑GGG‑GGGGTTAAA‑GGCAAAAAGGGGAC‑GTTCATTCCCAGAAGCA‑GC‑ggt                                                                                                              >  2:961707/1‑146 (MQ=255)
                                              ttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑TGC‑GGGGGTTAAAGCCCGAAAGGCGAC‑GTTTCTACCCCGACGCC‑GC‑GGTaaa                                                                                                           >  2:597092/1‑148 (MQ=255)
                                               tCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGc                                                                                                        <  2:965943/150‑1 (MQ=255)
                                                cAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCaa                                                                                                      <  1:114117/151‑1 (MQ=255)
                                                 aGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGGGTTAAATTCCAGAAAAGCGAC‑GGTCATCCCCGGACGCA‑GC‑GTTaaa                                                                                                           >  1:152095/1‑145 (MQ=255)
                                                 aGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGGGGTTAAAGGCAGCGAGGCGGC‑GTTCATTCCCGGACGCA‑GC‑GGTaaaaacaa                                                                                                      >  2:515611/1‑145 (MQ=255)
                                                  gCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGG‑GGGGGTAAAAGACAAAAAGGCGGC‑GTTTATCCCCAGGCGCA‑GC‑ggttaa                                                                                                           >  1:889/1‑139 (MQ=255)
                                                   ccGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAgcg                                                                                                   <  1:79163/151‑1 (MQ=255)
                                                    cGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑TGG‑GGGGTTTAATGCCCAATAAGCGAC‑GTTTATCCCCCGACGCA‑GC‑GGTAAACACCAgcg                                                                                                   >  2:526333/1‑150 (MQ=255)
                                                       aCCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑TGG‑GGGGGTTAATGCCCGAAAGGGGGC‑GTTCATCCCCGGAAGCA‑GC‑GGTTAAAACAAGAGCggg                                                                                               >  2:797407/1‑151 (MQ=255)
                                                        ccATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGgt                                                                                              <  2:513170/151‑1 (MQ=255)
                                                         cATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑TGG‑GGGGTTTAATGCCCGAAAAGCGGA‑GGTCTTTCCCCGGCGCA‑GC‑GGTTAAAGCAAACCCggg                                                                                               >  2:1320040/1‑149 (MQ=255)
                                                         cATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGgtc                                                                                             <  1:1666/151‑1 (MQ=255)
                                                              taGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑TGG‑GGGGGTTAATGGCCGAGAAGGGAC‑GGTTCTCCCACGACGCA‑GCGGTTAAAC‑CAAGCGCgggggcgccc                                                                                        >  2:1256913/1‑144 (MQ=255)
                                                               aGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑GTT‑TGC‑GGGGGTTAATTCCAGAAAAGGGAC‑GGTCATCCCCCGACGCC‑GG‑GGTTAACGCAAGCCCGGGGCGGCAc                                                                                        >  2:143156/1‑150 (MQ=255)
                                                                 aTGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTT‑GGG‑GGGGTTAAA‑GCCCAATAAGCGAC‑GGTCATCCCAAGGCGCA‑GC‑GGTTAAAACAAGCGCgggggg                                                                                            >  1:166677/1‑141 (MQ=255)
                                                                  tGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGGGTTAAATATCAGATAAGGGGC‑GGTCATCCCCCGACGCA‑GC‑GGTAAAAGCAAGAGCGGggcg                                                                                            >  1:462107/1‑143 (MQ=255)
                                                                   ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑TGCGGGGGTTAAATG‑CAGAAAAGGGACGGTTCATCCCCGGA‑GCA‑CC‑GGTAAAAGCCAGCGGGGGGCGGCACGc                                                                                      >  2:674654/1‑148 (MQ=255)
                                                                      ttGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGGG‑GGGGTTAATTACAAAAAAGGGGCG‑GTTCTTCCC‑GGACGCA‑GC‑GGTAAACGCAAGGGCGGGGGGTGACGCaaaa                                                                                  >  1:440040/1‑149 (MQ=255)
                                                                       tGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGG‑GGGGTTTAAAGGCACATATGAGAA‑GGTCATCCCAGGACGCA‑AC‑GTTAAAAGCAAGAGCGGGGCGGCCCCCaaaaa                                                                                 >  1:1287508/1‑148 (MQ=255)
                                                                       tGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATg                                                                                <  2:440519/150‑1 (MQ=255)
                                                                       tGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATg                                                                                <  2:10931/150‑1 (MQ=255)
                                                                          aTCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATa                                                                            <  2:1382476/151‑1 (MQ=255)
                                                                             aGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑TGC‑GGTGTTAAATGCCCGATAAGCGGC‑GTTCATCCCAGGGCGCA‑GC‑GGTAAAAGCAAACGCGGGGCGTCACGCCAAAttggta                                                                            >  1:942697/1‑143 (MQ=255)
                                                                                 gcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGG‑GGGGGTTAAAACCAGAGAAGAGAG‑GGTCCTCCCCGGGAGCAGGG‑GTTAACAGCAA‑CGCGGGGCGGCGCGCaaaaaaggaaacggg                                                                       >  2:986424/1‑138 (MQ=255)
                                                                                       tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcag                                                                <  1:854176/150‑1 (MQ=255)
                                                                                        attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcagc‑a                                                             <  1:475378/151‑1 (MQ=255)
                                                                                        attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcagc                                                               <  2:1204069/150‑1 (MQ=255)
                                                                                        attaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGgcagc                                                               <  2:737197/150‑1 (MQ=255)
                                                                                          taATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑TGC‑GGGGGTTAATAACAGATAAGGGAC‑GGTTTTCCCCGGGCGCA‑GC‑GGTAAACCCAAACGCGGGGCGGCCCGCCAAATGGATAGCTGGCTGGCAGc‑agg                                                           >  1:205209/1‑151 (MQ=255)
                                                                                            aTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGc‑aggcc                                                         <  1:797433/151‑1 (MQ=255)
                                                                                              cGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGc‑aggcca                                                        <  1:866756/150‑1 (MQ=255)
                                                                                                ccAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCATCACGCAAAATGGATAGCTGGCTGGCAGc‑aggccaggc                                                     <  2:1150777/151‑1 (MQ=255)
                                                                                                  aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGG‑GGGGGTTAAAGCCCGATAGGCGAC‑GTTCATCCCCCGGCGCA‑GC‑GGTAAAAGCAAACGCGGGGCGGCACGCAAAAAGGATAGGTGGGCGGCAGC‑aggccaggcc                                                    >  1:727894/1‑150 (MQ=255)
                                                                                                   gTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑Aggccaggcca                                                   <  1:815502/150‑1 (MQ=255)
                                                                                                     gCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGGGGTTAAAGCCAGATAAGCGGC‑GTTCATCCCCAGGCGCA‑GC‑GTTTAACACAAACGCGGGTCGTCACGCCAAAAGGATAGGTTGCTGGCAGC‑aggccaggcc                                                    >  1:886858/1‑147 (MQ=255)
                                                                                                          cgcATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCg                                           <  2:1204668/151‑1 (MQ=255)
                                                                                                             aTGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑TGC‑GGGGTTAAAAGCCAGATAGGCGAC‑GTTCATCCCCGGACCCA‑GC‑GGTTAAAGCCAACGCGGGGCGTCACGCAAAATTGGTAACTGGCTGGCAGC‑AGGGCAGGCCCGAACGTCGc                                          >  1:792294/1‑149 (MQ=255)
                                                                                                             aTGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCGCAt                                        <  2:833747/151‑1 (MQ=255)
                                                                                                              tGGCGGACAGCGGGATTATTTAGCGGCGG‑GTT‑GGG‑GGGGGTAAAATCCAGAAAGGCGAG‑GTTCATCCCCAGACGGA‑GC‑GTTTAACGCAAACGCGGGTCGTCACGCCAAATAGGGAACTGGCTGGCGGG‑AGGACCGGGCAGaaagccg                                           >  1:381563/1‑142 (MQ=255)
                                                                                                               ggCGGACAGCGGGATTATTTAGCGGCGGTTTT‑GGG‑GGGGTTAAA‑GGCAGGTAGGCGAC‑GTTTATTCCAAAAGGCA‑GC‑CATAAAAAAAAGCGGGGGGCGTGAAGCAAAAAAGATAAATAGCTGGCCGCAaggccgggcataaacacgaat                                         >  2:1083222/1‑137 (MQ=255)
                                                                                                                  ggACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCGCATTcccc                                   <  2:887058/151‑1 (MQ=255)
                                                                                                                            attattTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGt                         <  1:953810/151‑1 (MQ=255)
                                                                                                                               attTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGttt                       <  1:1259577/150‑1 (MQ=255)
                                                                                                                                tttAGCTGCGG‑TTT‑GGC‑GGTGTTAAAGGCCAGATAGGCGAC‑GTTCAGCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTc                      <  2:468761/150‑1 (MQ=255)
                                                                                                                                     cggcgg‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGa                <  2:1207807/151‑1 (MQ=255)
                                                                                                                                        cgg‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCCGGCCAGAACGGCGCATTCCCCGGTCAGGGGTTTCGACAg                 >  2:219990/1‑147 (MQ=255)
                                                                                                                                        cgg‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACg              >  2:422652/1‑150 (MQ=255)
                                                                                                                                          g‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTcc           >  1:1231218/1‑151 (MQ=255)
                                                                                                                                             tt‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCgg         >  1:1233407/1‑151 (MQ=255)
                                                                                                                                                     gtgtTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgc   >  2:1191789/1‑151 (MQ=255)
                                                                                                                                                      tgtTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCGCATTCCCCGGGCAGGAGTTTCGACAGACGTCCGGTgccgcc  >  1:568572/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                          
CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTT‑GGC‑GGTGTTAAATGCCAGATAGGCGAC‑GTTCATCCCAGGACGCA‑GC‑GTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCTGGCTGGCAGC‑AGGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGCC  >  NC_000913/1356054‑1356351

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