Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,157,271	0	T	G	42.4%	92.2 / 15.9	60	G165G (GGT→GGG)	ycfH	putative DNase
Reads supporting (aligned to +/- strand):  ref base T (4/30);  new base G (25/0);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 7.84e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCG  >  NC_000913/1157122‑1157420                                                                                                                                                       |                                                                                                                                                           catATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTa                                                                                                                                                           <  1:908103/151‑1 (MQ=255)  atatCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTaa                                                                                                                                                          <  1:514177/151‑1 (MQ=255)    atCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAAt                                                                                                                                                        <  2:1076276/151‑1 (MQ=255)        aGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑Ct                                                                                                                                                   <  1:495261/151‑1 (MQ=255)        aGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGGAAtt                                                                                                                                                       >  2:1123045/1‑148 (MQ=255)          aTTGGCCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGTCGATACGATGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTgg                                                                                                                                                 <  1:496396/151‑1 (MQ=255)              gTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATA‑CTGGGtc                                                                                                                                              >  1:629612/1‑150 (MQ=255)               tCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTA‑TTgggttt                                                                                                                                             >  1:66149/1‑146 (MQ=255)                cGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTGCAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGAtctc                                                                                                                                            <  1:684333/150‑1 (MQ=255)                cGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGAtctc                                                                                                                                            <  1:74267/150‑1 (MQ=255)                   gaactgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGGAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGAt                                                                                                                                        <  1:398850/151‑1 (MQ=255)                     actgaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTTCATTATCTGGaatgggagtt                                                                                                                                       >  2:64368/1‑142 (MQ=255)                        gaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGAtttt                                                                                                                                     <  1:18600/149‑1 (MQ=255)                        gaacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGAtttt                                                                                                                                     <  2:947776/149‑1 (MQ=255)                         aacAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATAGCGGTGGCGTCCTACTCTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACa                                                                                                                                  <  1:92485/151‑1 (MQ=255)                          acAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGGAAATA‑ATGGGTCTCGGGtttt                                                                                                                                     >  1:1037135/1‑147 (MQ=255)                          acAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGCAACTTG‑CTTGATCTTGGTtttt                                                                                                                                     >  2:552940/1‑147 (MQ=255)                           cAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACAtc                                                                                                                                <  2:942654/151‑1 (MQ=255)                            aaGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATA‑ATGGGTGTTGGAatttt                                                                                                                                    >  1:965905/1‑145 (MQ=255)                                cGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑tt                                                                                                                          <  2:341726/151‑1 (MQ=255)                                 ggTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTA‑CTGGGTCTCGGGTTTTTCATATCC‑ctt                                                                                                                         >  1:114024/1‑151 (MQ=255)                                   ttATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTCAATTA‑ATAGGCCTCTGGTTTTTAATCTCC‑tttt                                                                                                                        >  2:650054/1‑150 (MQ=255)                                     aTCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCg                                                                                                                     <  2:973406/151‑1 (MQ=255)                                     aTCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGTAAATA‑CTTGATTTCGGATTTTTAATTTTC‑tttt                                                                                                                        >  1:708881/1‑148 (MQ=255)                                       cGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGc                                                                                                                   <  1:277931/151‑1 (MQ=255)                                        gTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCa                                                                                                                  <  2:308769/151‑1 (MQ=255)                                        gTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATA‑CTGGGTCTCGGGTTTTATATCTCC‑tttttccgc                                                                                                                   >  1:1208813/1‑145 (MQ=255)                                          ccATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCAt                                                                                                                 <  2:1045578/150‑1 (MQ=255)                                           cATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTA‑ATGGGTTTCGGGTTTTTAATTTTC‑tttttccgcct                                                                                                                 >  2:613460/1‑142 (MQ=255)                                                    gACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGGAAATA‑ATGGATCTCGGGTTTTACATCTCC‑TTTTTCGGCATT‑GTG‑ACCTTcc                                                                                                    >  1:17436/1‑151 (MQ=255)                                                       gCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGCGGGGTAAATTA‑‑TGGATCTCAGGTTTTTTAACTCCTTTTTCCGCCATC‑TTG‑GCCTTc                                                                                                     >  2:30937/1‑148 (MQ=255)                                                        cccGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTaa                                                                                                <  2:531661/151‑1 (MQ=255)                                                             gCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATA‑CTTGATCTCGGATTTTTCATCTCC‑CTTTTCCGCATT‑TTG‑ACCTTCCGGCATGCgg                                                                                           >  1:1242844/1‑151 (MQ=255)                                                               cGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAg                                                                                         <  1:1191190/151‑1 (MQ=255)                                                                 aTACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTA‑CTTGATCTCGGGTTTTTTATCTCC‑TTTTTCCGCATT‑TTG‑ACCTTCCCTCATGCgg                                                                                           >  1:627499/1‑147 (MQ=255)                                                                    cGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGGAAATA‑TCGGAGCTCGGATTTTACAAATCC‑TTTTTCGGCAAT‑GTG‑ATCTTCCGTAATgaggag                                                                                         >  2:1233177/1‑146 (MQ=255)                                                                              tCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATA‑CTTGGTCTTCGGTTTTAAATCTCC‑TTTTTCCGCACT‑TTG‑GCCCTTCGTAAAGCGGAGCAAATGCGCGAt                                                                             >  2:963277/1‑148 (MQ=255)                                                                               cTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGGAATTA‑CTGGGTCTCGGGTTTTTTATTCCC‑TTTTTCCGCCCT‑GTG‑GGCCTTCGTAATGCGGGGCAAATCcgcggt                                                                             >  2:208515/1‑145 (MQ=255)                                                                               cTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTA‑CTGGATCTCGGGGTTTACATCTCC‑TTTTTCCGCCTT‑GTG‑ACCTTCCCGAATGCGGGGCAACTGCGCGAtgctg                                                                         >  1:630352/1‑151 (MQ=255)                                                                                  cgcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTTGGGCTCGGTTTTTATATCTCC‑CTTTTCGGCATT‑GTGTACCTT‑CGTGATGTGGAGCAACTGcgcggttgtggg                                                                       >  1:911889/1‑142 (MQ=255)                                                                                  cgcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTA‑CTGGATCTCGGGTTTTTTATCTCC‑TTTTTCCGCCTT‑TTG‑ACCTTCCGTAATTCGGGGCAACTGCGCGAtg                                                                            >  1:521446/1‑145 (MQ=255)                                                                                   gcgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAAATA‑ATGGATCTCCGGTTTTACATCTCC‑TTTTCCGGCATTGGTG‑A‑CTTCCCTAATGCGGAggcacagc                                                                                  >  2:627861/1‑131 (MQ=255)                                                                                    cgAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAGATTA‑ATGGATCTCTGATTTTTCATCTCCTTTTTCCGGCATG‑TTG‑GCCTTCC‑TAAAGCGGGGCAACTGCGGGATGTTGcg‑cg                                                                    >  1:315995/1‑150 (MQ=255)                                                                                       agaaAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTAt                                                                <  1:136613/151‑1 (MQ=255)                                                                                            aaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTccc                                                           <  2:1259972/151‑1 (MQ=255)                                                                                            aaGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTA‑CTGGGTCTCGGGTTTTTCATCTCC‑TTTTTCCGCCTT‑TTG‑GCCTTCCGTAATGCGGGGCAACTGCGCGGTGCTGCG‑CGGTTTTTccc                                                           >  1:165259/1‑151 (MQ=255)                                                                                             aGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGGGAAATA‑CTGGGTCTCGGGTTTTTTATCTCC‑CTTTTCCGCATT‑GTG‑GCCCTTCGTAATTCGGAGGAAATGCGCGATGGTTGG‑CGttttttctcc                                                          >  2:819909/1‑143 (MQ=255)                                                                                                gACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCt                                                        <  1:283861/150‑1 (MQ=255)                                                                                                    gATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATc                                                   <  2:627559/151‑1 (MQ=255)                                                                                                             ggCGTACTACACTGTTTTACAGAGGACAGTGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACtg                                           <  1:736723/150‑1 (MQ=255)                                                                                                                      cacTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTA‑CCGGACCCCCGGTTTTAAATCTCC‑TTTTTCGGCGTT‑GTG‑AACCTCCGCAAATCGGGGGAACTGCGCGGGGCTGCG‑CGTTTTGTCTCCCCGGGTCGGGTACTGGTGGGAAATg                                 >  2:702276/1‑151 (MQ=255)                                                                                                                       acTGTTTTACAGAGGACAGAGAAACGGC‑GGGGAAATTA‑CTGGGTCTCGGGTTTTATATCTCC‑TTTTTCCGCATT‑GTG‑ACCCTCCGTAATGCGGGGGAACTGCGCGGTGCTGCG‑CGGTATTTCCCCCCCGGGCGGGTACTGGGGGgaa                                    >  1:380403/1‑147 (MQ=255)                                                                                                                            tttACAGAGGACAGAGAAACGGC‑GGGGAAAATT‑CTTGATTTCGGGGTTTTCATCCCCTTTTTCCGGCGTT‑‑TG‑ACCCTTCCTAAAGCGGAGGAAATTCGCGCAGCCGGGCCGTTTGTTCCCCCGGGACGGGTTA‑TGGGGGGAACTGACt                              >  1:74448/1‑148 (MQ=255)                                                                                                                              tACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCacc                          <  1:619641/150‑1 (MQ=255)                                                                                                                                    ggACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCaccttacctt                   <  2:1251066/151‑1 (MQ=255)                                                                                                                                         gagaAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGTCTCACCTTACCTTGCGcc              <  1:1013052/151‑1 (MQ=255)                                                                                                                                         gagaAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGcc              <  1:667149/151‑1 (MQ=255)                                                                                                                                            aaaCGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGt           <  1:208521/151‑1 (MQ=255)                                                                                                                                                ggc‑ggGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCg       <  1:1036938/151‑1 (MQ=255)                                                                                                                                                ggc‑ggGGAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCg       >  2:470085/1‑151 (MQ=255)                                                                                                                                                       tAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCg  >  2:63352/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                           CATATCCAGATTGGTCGTGAACTGAACAAGCCGGTTATCGTCCATACCCGTGACGCCCGTGCCGATACGCTGGCAATTCTGCGCGAAGAAAAAGTGACGGATTGCGGTGGCGTACTACACTGTTTTACAGAGGACAGAGAAACGGC‑GGGTAAATTA‑CTGGATCTCGGATTTTACATCTCC‑TTTTCCGGCATT‑GTG‑ACCTTCCGTAATGCGGAGCAACTGCGCGATGCTGCG‑CGTTATGTCCCCCTGGATCGGTTACTGGTGGAAACTGACTCACCTTACCTTGCGCCGGTACCGCATCG  >  NC_000913/1157122‑1157420

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