Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,356,2010TG41.9% 188.3 / 17.0 93T304P (ACC→CCC) sapAantimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (6/48);  new base G (39/0);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 3.32e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGC  >  NC_000913/1356051‑1356350
                                                                                                                                                        |                                                                                                                                                            
ccTCGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑Ggt                                                                                                                                                              <  1:77924/151‑1 (MQ=255)
   cGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATCGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGtt                                                                                                                                                           >  4:411105/1‑151 (MQ=255)
   cGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑Ggtgt                                                                                                                                                            <  3:33330/150‑1 (MQ=255)
   cGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑Ggtg                                                                                                                                                             <  1:1157160/149‑1 (MQ=255)
   cGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGtt                                                                                                                                                           <  3:795936/151‑1 (MQ=255)
         cATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTtaa                                                                                                                                                        >  3:1406434/1‑148 (MQ=255)
          aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGG‑Ggggtgtaa                                                                                                                                                        >  1:643924/1‑143 (MQ=255)
            aGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCa                                                                                                                                                  <  3:359546/150‑1 (MQ=255)
            aGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCa                                                                                                                                                  <  2:505564/150‑1 (MQ=255)
            aGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAg                                                                                                                                                 <  3:736745/151‑1 (MQ=255)
             ggCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAg                                                                                                                                                 <  4:330274/150‑1 (MQ=255)
             ggCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGa                                                                                                                                                <  4:276189/151‑1 (MQ=255)
             ggCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑ATTGGC‑GGGGTTAAAA‑Gcc                                                                                                                                                   >  2:515146/1‑148 (MQ=255)
              gCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTTAAT‑gac                                                                                                                                                   >  3:251816/1‑145 (MQ=255)
                ccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTTAAT‑GCCAGa                                                                                                                                                >  3:758557/1‑148 (MQ=255)
                ccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑g                                                                                                                                            <  3:1064117/151‑1 (MQ=255)
                 cACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTAAAT‑GCCAGAAA‑gg                                                                                                                                           >  2:423030/1‑151 (MQ=255)
                     aGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTaaaatgcaaaaaggcgg                                                                                                                                          >  1:12401/1‑136 (MQ=255)
                      ggCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACg                                                                                                                                      <  1:628443/151‑1 (MQ=255)
                       gcggggggTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGTTTGCAGCGTGCGCTGGTTATTAATCGCCTGTGCCAGCGCATGACGGACAGTGGGATTAGTTACTGGCGT‑TTTGCC‑GGTGTAAAAG‑GCTAAATC‑GGCGACg                                                                                                                                      <  1:126852/144‑1 (MQ=255)
                       gcgcgcGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGTTTAAT‑GCCAGAAA‑AGCGAc                                                                                                                                       >  3:1291716/1‑149 (MQ=255)
                       gcgcgcGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGGTAAAA‑GCCAGATA‑gg                                                                                                                                           >  2:104574/1‑145 (MQ=255)
                         gcgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTAAAT‑GCCAGAGA‑Ggcggcggt                                                                                                                                    >  1:279436/1‑148 (MQ=255)
                           gcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTAAAT‑GCCAAATA‑AGCGGCGtt                                                                                                                                    >  3:342757/1‑148 (MQ=255)
                                  aaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGa                                                                                                                          <  4:268280/151‑1 (MQ=255)
                                   aTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTTAAT‑GCCAGATA‑GGCGGCGTTCATCCCCgg                                                                                                                           >  3:287956/1‑149 (MQ=255)
                                    tAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTAAAT‑GCCCGATA‑GGGGCCGTTCATCCCCGGAc                                                                                                                         >  4:223698/1‑150 (MQ=255)
                                     agagGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTTAAA‑TCCAGATA‑AGCGACGTTCATcccccg                                                                                                                           >  3:623009/1‑144 (MQ=255)
                                      gagGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGCCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGc                                                                                                                       >  3:949637/1‑150 (MQ=255)
                                          ccgccgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTAAAT‑TCCAGATA‑AGCGGCGGTCATCCCACGACGccg                                                                                                                     >  1:1138458/1‑146 (MQ=255)
                                              cgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTaaa                                                                                                             <  4:522734/151‑1 (MQ=255)
                                              cgTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGGTTAAT‑GCCAGAAA‑AGGGGCGGTCATCCCCCGACGCAGC‑ggt                                                                                                                >  1:97165/1‑146 (MQ=255)
                                               gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAAc                                                                                                            <  2:142809/151‑1 (MQ=255)
                                               gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGGTTAAT‑GCCAGATA‑AGGGACGTTCCTTCCCAGACGCcgc‑ggttaa                                                                                                             >  3:798664/1‑145 (MQ=255)
                                                 ttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACg                                                                                                           <  1:1284713/150‑1 (MQ=255)
                                                   cAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCaa                                                                                                        <  3:467793/151‑1 (MQ=255)
                                                     gCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGGTAAAA‑GCCCGATA‑AGCGGCGTTTATCCCCCGACGCCGC‑GGTTAACGc                                                                                                          >  1:81618/1‑147 (MQ=255)
                                                        gTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTTAAA‑GCCAGAAA‑AGCGACGCTCATCCCCGGACGCAGC‑GTTTAAAACAAacgc                                                                                                    >  1:256609/1‑150 (MQ=255)
                                                        gTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGGTAAAT‑GCCAGATA‑GGGGAGGTTCATCCCAGGACCCAGCGGTTaaagc                                                                                                           >  4:144211/1‑142 (MQ=255)
                                                            cATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTTAAT‑GCCCGATA‑GGCGGCGTTCATCCCCGGACGCAGC‑GGTAAACGCAAGCGCggg                                                                                                 >  3:1035583/1‑149 (MQ=255)
                                                            cATAATAGATGGATTGCATCAGGCGCGGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGgtc                                                                                               <  1:169689/151‑1 (MQ=255)
                                                             ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGgtcg                                                                                              <  3:80670/151‑1 (MQ=255)
                                                              taataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGG‑GGGGTTTAAA‑GCCCGAAA‑GGGGGCGGTTATCCCCAGACGCCGCGGTTAAACGaaaccaggg                                                                                                  >  4:267186/1‑141 (MQ=255)
                                                               aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGAACGCAGC‑GTTAAACGCAAGCGCGGgtcgt                                                                                             <  1:911335/150‑1 (MQ=255)
                                                                ataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑ACCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCa                                                                                           <  3:1019251/151‑1 (MQ=255)
                                                                ataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑ATTTGC‑GGGGGTTAAA‑GACAGAAA‑GGTGACGTTCATCCCCGGAGCCAGC‑GTTAAAAGCAAACGCggg                                                                                                 >  3:1167890/1‑145 (MQ=255)
                                                                   gATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑c                                                                                      <  1:1405804/151‑1 (MQ=255)
                                                                      ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGG‑GGGGTTAAAA‑GCCAGATA‑AGCGGCGGTCATCCCAGGACGCAGC‑GGTAAAAGCAAACGCGGGTCGGCAC‑G‑Caa                                                                                    >  2:32774/1‑150 (MQ=255)
                                                                      ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑Caaa                                                                                   <  3:268310/151‑1 (MQ=255)
                                                                      ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑Caaa                                                                                   <  2:273224/151‑1 (MQ=255)
                                                                      ggATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑Caa                                                                                    <  1:257768/150‑1 (MQ=255)
                                                                       gATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCTG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑Caaaa                                                                                  <  3:145794/151‑1 (MQ=255)
                                                                       gATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGGTAAAA‑GGCAGATA‑GGCGACGTTTATCCCCGGACGCGGC‑GTTAAAAGCAAGCGCGGGGCGGCAC‑G‑Caa                                                                                    >  2:70307/1‑149 (MQ=255)
                                                                          tGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATg                                                                                <  1:79296/150‑1 (MQ=255)
                                                                          tGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATg                                                                                <  2:307493/150‑1 (MQ=255)
                                                                           gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGa                                                                              <  2:163026/151‑1 (MQ=255)
                                                                            cATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGAt                                                                             <  3:416034/151‑1 (MQ=255)
                                                                            cATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGTTAAAA‑GCCAGGTA‑GGCGACGTTCATCCCCCGGAGCAGC‑GTTAAAACCAAGCGCGGGGCGTCCC‑C‑CAAAATgg                                                                               >  2:277061/1‑149 (MQ=255)
                                                                                    gcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAA‑GCCAGATA‑GGCGACGTTCATCCCCGGACGCAGC‑GTTAAAAGCAAGCGCGGGGCGTCAC‑G‑CAAAAAGGATAgct‑gg                                                                      >  2:541936/1‑149 (MQ=255)
                                                                                      tGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTTAAAA‑ACAAAAAA‑AG‑GACGTTTATCCCAGGGCGCAGC‑GGTAAAACCAAAGGCGGGGCGCCAC‑C‑Caaaaatggtaactg                                                                        >  4:63661/1‑136 (MQ=255)
                                                                                       ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGct‑ggctgg                                                                  <  3:32470/150‑1 (MQ=255)
                                                                                             taATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGG‑GGGGTTAAATAGCCA‑AAA‑GGGGGCGTTTATCCCCGGACGCAAC‑GTTAAAAGCAAGCGCGGGGCGTCAA‑G‑CCAAAAGGATAGgt‑ggttgg                                                                  >  3:297923/1‑144 (MQ=255)
                                                                                             taATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGcagcag                                                            <  1:1270338/150‑1 (MQ=255)
                                                                                             taATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGGTAAAA‑GCCAGCTA‑GGCCACGTTTATCCCCGGGCGCAGG‑GGTAAACGCCAGCGCGGGGCGGCCC‑G‑CAAAAAGGATAGCT‑GGCGGgca                                                                >  4:365463/1‑146 (MQ=255)
                                                                                              aaTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTAAAA‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GGTAAAAGCAAACGCGGGGCGTCAC‑G‑CAAAATGGAAAGct‑gggtgg                                                                  >  2:402608/1‑143 (MQ=255)
                                                                                              aaTCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGcaggc                                                          <  3:234349/151‑1 (MQ=255)
                                                                                                 cGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGcaggcca                                                        <  4:405392/150‑1 (MQ=255)
                                                                                                  gCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTAAAT‑GGCAGAAA‑GGCGGCGTTCATCCCAGGACGCAGC‑GGTAAAAGCAAGCGCGGGGCGGCAC‑G‑CAAAATGGGTAGCT‑GGGTGGGAGcagggc                                                         >  4:560594/1‑146 (MQ=255)
                                                                                                  gCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGcaggccagg                                                      <  3:133709/151‑1 (MQ=255)
                                                                                                  gCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGcaggccagg                                                      <  3:901945/151‑1 (MQ=255)
                                                                                                    cAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGTTTAAA‑GCCAGAAA‑AGCGGAGGTCATCCCCGGACGCAGC‑GGTTAACGCAAGCGCGGGGCGGCAC‑G‑CCAAAAGGATAGCT‑GGGTGGCAGcaggcccgg                                                      >  2:128102/1‑149 (MQ=255)
                                                                                                     aGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAggccaggcca                                                   <  3:479843/151‑1 (MQ=255)
                                                                                                           aGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGG‑GGGGTTAAA‑‑GGGAGAAAGGGGGAGGTTATTCCCCGG‑CGCAGC‑GGTTAAAGCAAAAGCGGGGCGGGAC‑G‑CCAAAAGGATAGCT‑GGGCGGGAGcagggccgggcaaaa                                                >  1:611356/1‑137 (MQ=255)
                                                                                                           aGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTTAAT‑GCCAGATA‑AGCGGCGTTTATCCCAGGACGCCGC‑GTTAAAAGCAAGCGCGGGGCGTCAC‑G‑CCAAATGGAAAGCT‑GGCTGGCAGCcggcccggcc                                                    >  3:1301163/1‑144 (MQ=255)
                                                                                                            gcgcATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGTGTTAAAA‑GGCATACA‑GGGGAACGTCAACCCAAAACCAATG‑CTTAAAACCAAGCGCGGGTCGTTGCGG‑CAAAAAGG‑TAACC‑GGCTGTCCGCaggtccggcccc                                                  >  3:480576/1‑143 (MQ=255)
                                                                                                             cgcATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCTGC‑GTCAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCg                                           <  2:166026/151‑1 (MQ=255)
                                                                                                              gcATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑c                                         <  2:466407/151‑1 (MQ=255)
                                                                                                                aTGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTAAAA‑GCCAGAAA‑GGCGGCGTTCATCCCAAGACGCAGC‑GGTAAAAGCCAGCGCGGGGCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGGCCGAAAGTCg                                           >  3:569778/1‑148 (MQ=255)
                                                                                                                   gCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGGTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAAAGCAAGCGCGGGGCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTcc                                    >  4:48536/1‑151 (MQ=255)
                                                                                                                    cGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTccc                                   <  3:12400/151‑1 (MQ=255)
                                                                                                                    cGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTccc                                   <  2:155725/151‑1 (MQ=255)
                                                                                                                     ggACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTTAAA‑TCCCGAAA‑GGCGACGTTTATTCCCAGACGCAGC‑GGTAAAAGCCAACGCGGGGCCGCAC‑C‑CAAAATTGGTAGCTGGGCCGGC‑GCCAGCCAGGCCCGAACGTCG‑CCTTtccc                                  >  3:241003/1‑151 (MQ=255)
                                                                                                                     ggACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTcccc                                  <  2:193004/151‑1 (MQ=255)
                                                                                                                       aCAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTTAAT‑GCCAGATA‑GGCGGCGGTCATTCCAGGGCGCAGC‑GGTGAAAGCCAGACCGGGGTGTCAC‑G‑CAAAAAGGATAGAC‑GGGCGGGAGCAAGGCCGGCCAGaaagacggctttccc                                   >  1:763625/1‑136 (MQ=255)
                                                                                                                          gCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCa                             <  2:560602/151‑1 (MQ=255)
                                                                                                                             ggATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGa                          <  1:319267/151‑1 (MQ=255)
                                                                                                                             ggATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGa                          <  3:842965/151‑1 (MQ=255)
                                                                                                                                  attTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGAGTTTc                     <  1:1098163/151‑1 (MQ=255)
                                                                                                                                  attTAGCGGCGG‑TTTGGC‑GGGGTTAAAA‑GTCAGAAA‑GGGGACGTTCATCCCCCGGCGCAGC‑GGTAAAAGCAAAGGCGGGGCGGCAC‑GCCAAAAAGGAAAGCT‑‑GCTGGCAGCAAGGCACGCCCGAAAGGCGCCCTTCCCCGGCCGGGGGttt                      >  1:695762/1‑151 (MQ=255)
                                                                                                                                      aGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGAGTTTCGAc                  <  3:100705/150‑1 (MQ=255)
                                                                                                                                      aGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGAGTTTCGACa                 <  2:345411/151‑1 (MQ=255)
                                                                                                                                         ggcgg‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGAGTTTCGACAg                >  4:381476/1‑149 (MQ=255)
                                                                                                                                           cgg‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGAGTTTCGACAGACGt            >  3:957350/1‑151 (MQ=255)
                                                                                                                                            gg‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGAGTTTCGACAGACGTc           <  3:561264/151‑1 (MQ=255)
                                                                                                                                             g‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGAGTTTCGACAGACGTcc          >  3:1225988/1‑151 (MQ=255)
                                                                                                                                             g‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGAGTTTCGACAGACGTcc          <  3:335439/151‑1 (MQ=255)
                                                                                                                                                  ggC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGt       <  4:297868/150‑1 (MQ=255)
                                                                                                                                                       gtgtTAAAT‑GCCAGCTA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CCTTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgc  >  4:301833/1‑151 (MQ=255)
                                                                                                                                                       gtgtTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgc  >  1:58345/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                            
CCTCGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAT‑GCCAGATA‑GGCGACGTTCATCCCAGGACGCAGC‑GTTAAACGCAAGCGCGGGTCGTCAC‑G‑CAAAATGGATAGCT‑GGCTGGCAGCAGGCCAGGCCAGAACGTCG‑CATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGC  >  NC_000913/1356051‑1356350

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: ‑