Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,356,201	0	T	G	42.6%	119.8 / 12.1	68	T304P (ACC→CCC)	sapA	antimicrobial peptide transport ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base T (8/31);  new base G (29/0);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 4.34e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCT‑GGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGC  >  NC_000913/1356054‑1356350                                                                                                                                                      |                                                                                                                                                          cGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑Ggtgt                                                                                                                                                         <  1:269488/150‑1 (MQ=255)  gTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTa                                                                                                                                                       <  1:646746/151‑1 (MQ=255)     gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAAt                                                                                                                                                    <  2:919468/151‑1 (MQ=255)     gTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGTTaaa                                                                                                                                                     >  2:122854/1‑150 (MQ=255)     gTCATAGGCCCACGAGGCGCGCGGTAAAATAAAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTaaa                                                                                                                                                     <  1:887105/150‑1 (MQ=255)        aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGGTTTTGGG‑GGGGTTaaagac                                                                                                                                                  >  1:932692/1‑148 (MQ=255)        aTAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑GTTTGC‑GGGGGTAAATggc                                                                                                                                                 >  1:230413/1‑149 (MQ=255)          aGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTTAATGcc                                                                                                                                                 >  2:634306/1‑149 (MQ=255)           ggCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑a                                                                                                                                             <  2:944257/151‑1 (MQ=255)            gCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGG‑GGGGTTAAAAGCCa                                                                                                                                                >  1:1137717/1‑148 (MQ=255)             cccACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGG‑GGGGTTAAATGCCAG‑a                                                                                                                                             >  1:45030/1‑149 (MQ=255)                   aGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGTGGTAAATGCCAG‑ATa                                                                                                                                           >  2:342486/1‑145 (MQ=255)                   aGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGAc                                                                                                                                     <  2:526814/151‑1 (MQ=255)                       gcgcgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGG‑GGGGGTAAATGCCAG‑aaaaggggaggt                                                                                                                                  >  1:1155482/1‑139 (MQ=255)                          cgGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGCGGGGGTTAAAGGCC‑G‑AAAAGCGGCGTTCa                                                                                                                                >  1:695594/1‑149 (MQ=255)                           ggTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATca                                                                                                                             <  1:986891/151‑2 (MQ=255)                            gTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATccc                                                                                                                            <  2:582284/151‑1 (MQ=255)                              aaaaTAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGAGTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCa                                                                                                                           <  1:249758/150‑1 (MQ=255)                                    gagGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCa                                                                                                                   <  1:397460/151‑1 (MQ=255)                                           gccgttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTaa                                                                                                            <  2:76651/147‑1 (MQ=255)                                             gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAAc                                                                                                          <  2:624533/151‑1 (MQ=255)                                             gTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAAc                                                                                                          <  2:1155511/151‑1 (MQ=255)                                               ttCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACg                                                                                                         <  2:445127/150‑1 (MQ=255)                                                tCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGc                                                                                                        <  1:991582/150‑1 (MQ=255)                                                tCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGGTAAATGACAG‑ATAGGGGGAGTGCATCCCAGGG‑CGCAGCGTTaaa                                                                                                           >  2:165559/1‑147 (MQ=255)                                                    ccGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGGTTAAAGTCCG‑ATAAGGGGAGTTCATCCCAGGA‑CGCAGCGTTTAACGCaa                                                                                                      >  1:499304/1‑148 (MQ=255)                                                     cGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTAAATGCCAG‑ATAGGCGGCGTTTATCCCCGGG‑CGCAGCGGTAAAAGCaa                                                                                                      >  1:21193/1‑147 (MQ=255)                                                     cGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGTTAAATGCCAG‑ATAGGCGGCGTTCATCCCCGGA‑CGCAGCGTTAAAAGc                                                                                                        >  2:612993/1‑145 (MQ=255)                                                      gTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGcgcg                                                                                                 <  2:1031732/151‑1 (MQ=255)                                                      gTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGTTAAATGCCAG‑ATAAGCGGCGGTCATCCCCCGA‑CGCAGCGGTaaa                                                                                                           >  1:165350/1‑141 (MQ=255)                                                       aaCCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGG‑GGGGTTAAAAGCCAG‑AAAAGCGACGTTCATCCCAGGA‑CGCAGCGTTAAAAGCAAAGGCgg                                                                                                >  2:968788/2‑151 (MQ=255)                                                          cATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGTTAAAAGCCCA‑ATAGGGGAAGTTCATTCCAGGG‑CGCAGCGTTTAACGCaaaccccgg                                                                                               >  1:628425/1‑142 (MQ=255)                                                           ataataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGgtcg                                                                                            <  2:185022/151‑1 (MQ=255)                                                             aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTTAATTACAGAAAAGGCG‑CGTTTATTCCCAGA‑CCCAGCGGTAAAAACAAACGCggg                                                                                               >  1:647828/1‑146 (MQ=255)                                                             aataGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGgtcgt                                                                                           <  1:93807/150‑1 (MQ=255)                                                                 gATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAAGCAGCGCTCATCCCCGGT‑AGCAACGTTAAACCCAAATGCGGGTCGTCACGc                                                                                      >  2:171796/1‑151 (MQ=255)                                                                  aTGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGG‑GGGGGTAAAAGCCAG‑ATAAGCGGCGTTTATCCCCGGA‑CGCAGCGGTAAAAGCAAAAGCgggggggca                                                                                         >  1:2451/1‑141 (MQ=255)                                                                  aTGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCa                                                                                     <  1:240773/151‑1 (MQ=255)                                                                       ttGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGGTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGGTAAACGCAAGCGCGGGGCGGCACGCaaaa                                                                                  >  2:130252/1‑149 (MQ=255)                                                                        tGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATg                                                                                <  2:1146276/150‑1 (MQ=255)                                                                         gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATgg                                                                               <  2:1147162/150‑1 (MQ=255)                                                                         gCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATgg                                                                               <  1:574007/150‑1 (MQ=255)                                                                             cAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGGGTTAAATGCCAG‑ATAGGCGACGGTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAAAGGATa                                                                            >  2:1064254/1‑149 (MQ=255)                                                                                 cgcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑ATTCGC‑GGTGTTGAAAACCAG‑GCAGGAGACCTTCATCCCAGGAGCCCGGCGTT‑AACGGCAGCGGGGGGCGGCGCGGCAAATGGATAgc                                                                          >  2:248112/1‑147 (MQ=255)                                                                                  gcTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAgct‑ggct                                                                    <  2:1144893/151‑1 (MQ=255)                                                                                     ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGct‑ggctgg                                                                  <  1:855449/150‑1 (MQ=255)                                                                                     ggTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑GTTGGC‑GGGGGTAAATGGCAA‑ATAAGGGGAGTTCATCCCAGGG‑CGCAACGGTTAAAGCCAGCCCGGGGCGTCACGCAAAAAGGAAAAct‑ggctgg                                                                  >  1:867834/1‑150 (MQ=255)                                                                                     ggTGATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGct‑ggctgg                                                                  <  1:113690/150‑1 (MQ=255)                                                                                        tattaATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGTTTAATGCCAC‑CTAAGCGCCGCTCATCCCCCGA‑CGCCCCGGTTAACGCCAACGCGGGTCGTCACCCAAAAAAGGTAGgt‑gggtgg                                                                  >  1:973399/1‑147 (MQ=255)                                                                                               cGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCT‑GGCTGGCAGca‑ggcca                                                       <  2:350214/150‑1 (MQ=255)                                                                                               cGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑CTTGGC‑GGGGGTAAAAGGCCG‑AGAGGGGGGGGTTATCCCCGCA‑CGCCCCGGTTAAAGCCAACGCGGGGCGTCACGCAAAAAAGATAGAG‑GGCTGGCAGga‑ggccc                                                       >  1:672347/1‑149 (MQ=255)                                                                                                  cAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGG‑GGTGGTAAATGGCAG‑ATAAGGGACGGTCATCCCAGGA‑CGCAGCGGTAAAAGCAAGCGCGGGGCGTCAACCAAAATGGGGAGCT‑TGGTGGGAGca‑gggc                                                        >  2:37732/1‑144 (MQ=255)                                                                                                  cAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTTAATGCTAG‑ATAGGCGGCGTTCATCCCCGGA‑CGCAGCGTTTAAAGGAAGCGCGGGGCGTCACGCAAAAAGGATAGCT‑GGCTGGGAGca‑ggcccgggc                                                   >  1:866291/1‑149 (MQ=255)                                                                                                        cAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCT‑GGCTGGCAGCA‑GGCCAGGCCAGAACg                                             <  2:479089/151‑1 (MQ=255)                                                                                                              aTGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCT‑GGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCAt                                       <  1:634297/151‑1 (MQ=255)                                                                                                                      cAGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTTAATATCCC‑ATAGGCGACGTTCATCCCCGGG‑CGCAGCGGTAAACGCAAGCGCGGGGCGTCCCCCAAAAAAGATAGCT‑GGCTGGCCGCA‑GGCCAGGCCCGAAACGCGCATTtccc                                  >  1:368959/1‑148 (MQ=255)                                                                                                                       aGCGGGATTATTTAGCGGCGG‑TTTTGC‑GGGGGTAAATGCCAG‑AAAGGGGGCCTTTATCCCCGGG‑CGCAGCGGTAAAAGCAAACGCGGGGGGTCACGCAAAAAGGATAGCT‑GGGTGGCAGCA‑GGACCGGGCACAACGGCGCAATcccc                                  >  1:362064/1‑147 (MQ=255)                                                                                                                        gCGGGATTATTTAGCGGCGG‑CTTGGA‑GGGGGTTAATTCCAG‑AGAGGCGAGGGTTTTTCCCGGG‑CGCCGCGGTAAAAACAAACCGGGGGCGTCCCGCCAAAAGGATAGATCGGCT‑GCAGCA‑AGCCCAGCCCGAACAACGCCTTcccc                                  >  1:242136/1‑146 (MQ=255)                                                                                                                          gggATTATTTAGCGGCGG‑TTTGGG‑GGGGGTTAAAGCCAG‑ATTGGCGACGTTTATCCCCGGA‑GGCAGCGGTAAAAACAAAAGCGGGGCGGCACCCAAAAAGGAAAGCC‑GGCTGGGAGCACGGCCCGCCCAAAACG‑CGCATTcccccg                                >  1:813504/1‑144 (MQ=255)                                                                                                                          gggATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCT‑GGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAgg                           <  2:974626/151‑1 (MQ=255)                                                                                                                           ggATTATTTAGCGGCGG‑TTTTGC‑GGGGTTAAATTCCCG‑AAAAGCGACGTTCATCCCAGGA‑CGCAGCGTTTAACGCAAACGCGGGGCGTCACGCCAAATGGATAGCT‑GGCTGGCAGCC‑GGCCCGGCCCGAACCGCGCATTTCCCCGTCCgg                           >  2:865175/1‑150 (MQ=255)                                                                                                                                attTAGCGGCGG‑TTTGGG‑GGGGTTTAATGCCAG‑ATAAGCGACGTTCATCCCAGGG‑CGCCGCGTTAAACGCAAACGCGGGGCGTCACGCAAAATGGATAGCT‑GGCTGGCAGCA‑GGCCAGGCCAGGACGTCGCATTCCCCCGGCAGGAGgttt                     >  1:212504/1‑150 (MQ=255)                                                                                                                                   tAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCT‑GGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGa                   <  1:871234/150‑1 (MQ=255)                                                                                                                                    aGCGGCGG‑TTTTGC‑GGGGTTAAATGCCAG‑ATAGGCGACGTTCATCCCCGGA‑CGCAGCGTTTAACGCAAGCGCGGGGCGTCACGCCAAAAAGATAGCT‑GGCTGGCAGCA‑CGCCAGGCCAGAACGTCGCATTCCCCCGTCAgggggtt                      >  1:688134/1‑141 (MQ=255)                                                                                                                                     gcggcgG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCT‑GGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGGTTCGAc                  >  2:72997/1‑149 (MQ=255)                                                                                                                                        gcgg‑TTTGGC‑GGTGTTAAATGCCAG‑ATAAGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCT‑GGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGa               >  2:1100192/1‑149 (MQ=255)                                                                                                                                          gg‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCT‑GGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTc           <  2:69623/151‑1 (MQ=255)                                                                                                                                                ggC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCT‑GGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGt       <  2:1093630/150‑1 (MQ=255)                                                                                                                                                  c‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCT‑GGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgcc    >  1:602495/1‑151 (MQ=255)                                                                                                                                                     gtgtTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCT‑GGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTgccgc  >  1:567885/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                          CGTTGTCATAGGCCCACGAGGCGCGCGGTAAAATAGAGGCCGCCGTTTCAGCCGTACCATAATAGATGGATTGCATCAGGCGCTGGTTATTAATCGCCAGTGCCAGCGCATGGCGGACAGCGGGATTATTTAGCGGCGG‑TTTGGC‑GGTGTTAAATGCCAG‑ATAGGCGACGTTCATCCCAGGA‑CGCAGCGTTAAACGCAAGCGCGGGTCGTCACGCAAAATGGATAGCT‑GGCTGGCAGCA‑GGCCAGGCCAGAACGTCGCATTCCCCGGTCAGGAGTTTCGACAGACGTCCGGTGCCGC  >  NC_000913/1356054‑1356350

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