Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	357,890	0	C	A	45.0%	68.5 / 37.8	60	L267F (TTG→TTT)	cynR	transcriptional activator of cyn operon; autorepressor
Reads supporting (aligned to +/- strand):  ref base C (30/3);  new base A (0/27);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 9.23e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA  >  NC_000913/357740‑358035                                                                                                                                                          |                                                                                                                                                    cTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGcaac                                                                                                                                                     >  2:611926/1‑151 (MQ=255)   gATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGcaaca                                                                                                                                                    >  1:950458/1‑150 (MQ=255)       aGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAcc                                                                                                                                               >  2:997443/1‑151 (MQ=255)         cGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgc                                                                                                                                             >  1:1044177/1‑151 (MQ=255)          gTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgcc                                                                                                                                            >  1:1117270/1‑151 (MQ=255)            aGCGCTTAAGGCACTTCCAGCCGCAGACCTGTGTCAGCGCTTACCGTGATTCATTCCCGCCAACACCCGCGCATTTCTCCAACGCC‑ATGTGAAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAAccgccg                                                                                                                                           <  2:947600/150‑1 (MQ=17)              cgcATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTtc                                                                                                                                        >  1:1004836/1‑151 (MQ=255)                        aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                                >  2:66890/1‑149 (MQ=255)                        aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                                >  2:323426/1‑149 (MQ=255)                        aaTTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCagt                                                                                                                                >  2:148421/1‑149 (MQ=255)                            ccAGCCGCAGCCCTGTGCCAGCGGCTCCCGTGATTCTTTTCCCCCAACACCCGCGCATTTACCCACCGCC‑TTGTGCAAAAATGCTTC‑CGCGGCGG‑TTGTCTGCCAGCTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                            <  1:392063/149‑1 (MQ=255)                            ccAGCCGCAGACTTGTGTCAGCGGCTACCGTGTTTCATTCCCCCCAACACCCGCACATTTATCCAACGCC‑ATGTGAAAAAATGCCTC‑CGCGGCGG‑CTTTCTGCCAACTATTTTTCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTg                                                                                                                            <  1:880399/149‑1 (MQ=255)                             cAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                                                        >  2:1027376/1‑151 (MQ=255)                                 cGCAGACCTGTGTCGGCGGCTCCGGTGTTTCATTCCCCCAAACACCCGCGTTTTTACCCACCGCC‑TTGTAAAAAAATGCCTT‑CGGGGGGG‑TTGTTTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑g                                                                                                                        <  1:703562/147‑1 (MQ=255)                                 cGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                                     >  1:591593/1‑150 (MQ=255)                                      aCCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGaa                                                                                                                  >  1:1335014/1‑148 (MQ=255)                                        cTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggg                                                                                                                     >  1:48741/1‑143 (MQ=255)                                           tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑ggggggggag                                                                                                               >  2:1028218/1‑140 (MQ=255)                                           tgtCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAGGA‑Gaa                                                                                                            >  1:112826/1‑148 (MQ=255)                                           tgccAGCGCCTCCGGTGTTTCATTCCCCCCAACACCCGCGCATTAACCCACCGCC‑ATGTGAAAAAATGCTTC‑CGCGGCGG‑CTGTCTGCCAGTTATTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAaga                                                                                                                <  1:693945/142‑1 (MQ=255)                                             tCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑Gaaaaa                                                                                                         >  2:18403/1‑147 (MQ=255)                                              cAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑Gaaaaa                                                                                                         >  1:807636/1‑146 (MQ=255)                                                ggggctccggtttttCTTTCCCGCAACCACCCGCGATTTTACCCACCGCC‑ATGTGAAAAAATGCCTT‑CGGGGCGG‑CTGTCTGCCGGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGc                                                                                                                    <  1:1328060/123‑1 (MQ=255)                                                 cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAggggc‑ggggccgag                                                                                                                >  2:372986/1‑135 (MQ=255)                                                 cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑gggaggagaaaaa                                                                                                            >  2:1162485/1‑133 (MQ=255)                                                 cGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑GGAAAAGAtt                                                                                                     >  1:153014/1‑149 (MQ=255)                                                  gtctcccGTTATTCATTCCCGCCAACACCCGCGCTTTTACCCACCGCC‑ATGTGAAAAAATCCCTC‑CGCGGCGG‑CTGTCTGCCAGCTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATa                                                                                                         <  1:1324420/139‑1 (MQ=255)                                                  ggcccccGTGATTCATTCCCCCCAACACCCGCGCATTTACCCAACCCC‑GTGTGAAAAAATGCCTC‑CGCGGGGG‑TTTTCTGCCAGCTATTTTCCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTa                                                                                                   <  1:651044/145‑1 (MQ=255)                                                   gtttccctttttCATTTCCCCCCAACACCCCCCCTTTTTCCAACCCCTATTTAAAAAAATCCTTC‑CGCGGCGG‑CTTTTTGCCGGTTTTTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCtt                                                                                                     <  1:555549/138‑1 (MQ=255)                                                       ccgtatttCATTCCCCCCAACACCGGCGCTTTTCCCCACCCCC‑TGGCAAAAAAGTCCTTT‑CGCGCGGT‑CTGTCCGCCACTTTTTTTCCCGCCAAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTaa                                                                                                  <  2:950437/140‑1 (MQ=255)                                                       ccGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTA‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGTAGA‑Gaaaaaatat                                                                                                     >  1:374777/1‑137 (MQ=255)                                                         gTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGA‑GAAAAAACCTTaaagccct                                                                                            >  1:602721/1‑144 (MQ=25)                                                          ttattcttcccccccAACACCCGCCCTTTCACCCACCGCT‑TTGTAAAAAAAGCCCTG‑CG‑GGCGGTCTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATa                                                                                                         <  2:1225140/123‑1 (MQ=255)                                                             ttCTTTCCCCCCAACACCCGCCCTTTTACCCACCCCC‑ATG‑GCAAAAAATCCTTCCGCGGGGG‑CTTTCGCCCATTTATTTTCCCCCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCAtg                                                                                        <  1:678902/150‑1 (MQ=255)                                                                aTTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGAGGAAAAAACTTTAAAccgtcacatt                                                                                       >  2:1418514/1‑141 (MQ=25)                                                                     ccccaccaCCCGCGATTTTACCAACCGCC‑TTGTGAAAAAATCCCTT‑CGCGGCGG‑CTGTCGGCCAGCTATTTTCCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGc                                                                               <  1:388291/149‑1 (MQ=255)                                                                      gCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAGGA‑GGAAAAACATTTAACACCTCATGTTTTGTGGcg                                                                              >  2:175200/1‑150 (MQ=25)                                                                         aacaccCCCCCATTTACCCACCCCC‑ATGGAAAAAAATCCCTC‑CGGGGCGG‑CTGTCTGCCAGTTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATc                                                                           <  2:580864/151‑1 (MQ=255)                                                                               cgcgCTTTTACCCACCGCC‑GTGAAAAAAAAGGCTTC‑GGGGGGGG‑CTGTCTGCCTGTTTTTTTCCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCgctgct                                                                     <  2:233205/151‑1 (MQ=255)                                                                               cgcgCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGCGGGGGCAGGG‑GAAGAAAATTTAAACCC‑TCATGTTGTGTgggagacg                                                                          >  2:701792/1‑140 (MQ=25)                                                                                     tttATCCACCGCT‑ATGTGAAAAAATCCCTG‑CGCGGCGG‑TTGTCTCCGAGTTTTTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTa                                                                 <  1:456604/149‑1 (MQ=255)                                                                                         ccaaccccc‑ATGGGAAAAATTGCCTC‑GGGGGCGG‑CTGTCTGCCAGTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACaa                                                             <  1:874915/142‑1 (MQ=255)                                                                                          ccaccccc‑ATGTGAAAAATTGCTTT‑CGGGGGGG‑TTGTCGGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAAtgt                                                          <  1:1055880/145‑1 (MQ=255)                                                                                           cAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGA‑GAAAAAAAATTTAAACCGCCCTGTTGTGTGGGGATCTCTGCggggtaaaa                                                             >  1:805257/1‑141 (MQ=25)                                                                                            accGCC‑AGGTAAAAAAATGCCTC‑GGCGGCGG‑CTGTCTGCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                        <  2:328025/149‑1 (MQ=255)                                                                                              cGCC‑TTGTGCAAAAATCCCTC‑CGGGGCGG‑CTGTCTCCCAGTTATTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                        <  1:525340/149‑1 (MQ=255)                                                                                                  aatgccaaaaaaTGCTTC‑CGCGGCGG‑CTGTCTGCCATTTTTTTTCCCGCCGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTgg                                                        <  2:33202/139‑1 (MQ=255)                                                                                                   aTGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGCAGG‑GAAAAAACATTTAACCCCACATATTTTGTTGGAAGCCCTGCCTGGAACAAAGTGGgaa                                                     >  2:657726/1‑148 (MQ=19)                                                                                                           aaaTGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGG‑GAAAAAACTTTTAACCCCTCCTGTTTTGTGGGGATCGCTGCTTGGGAAAAAGTGGGAAGGGGAgggggg                                          >  1:1183639/1‑146 (MQ=255)                                                                                                            aagcccgt‑CGGGGGGG‑CTGCCGGCCATATTTTTTCCCGCCGAAAAAAAGCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTgcggc                                         <  2:493698/144‑1 (MQ=255)                                                                                                               gCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAAA‑AAAAAAAATTTAAAACCCCCCTTTTGTGTGGGCAAAGCACCTGGGGACAATGGGGGAAGGGGAgggggg                                          >  1:342018/1‑142 (MQ=255)                                                                                                                      gcggcgG‑TTGTCTCCAATCTTTTTTTCCCCAGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCa                                   <  1:372726/148‑1 (MQ=255)                                                                                                                        ggcgg‑cTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGAAGG‑GGGAAAGCTTTTAGCCAAGCCTGTTTGGTGGGCATCGCTGGTTGGGACAAAGGGGAAAGGGGAGTggggggggacagagc                               >  2:448594/1‑139 (MQ=255)                                                                                                                                  cggCCAGTTTTTTTCCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAAcc                        <  1:18402/146‑1 (MQ=255)                                                                                                                                   tGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑GAAAAAAACTTAAAACCCTCCTTTTGTGTGGGGATAGCCGCTGGGGAAAAAGTGGGAAGGGGAGGGGGGGGGATCaaccccccgaacccc                     >  2:1402850/1‑136 (MQ=255)                                                                                                                                     ctgtctttttttCCGCAGAAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcg                       <  2:57875/139‑1 (MQ=255)                                                                                                                                     ccatttATTTTTCCCCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACcgc                      <  1:1407216/142‑1 (MQ=255)                                                                                                                                     caggtttttttcCCGCCGCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGaa                          <  2:590952/136‑1 (MQ=255)                                                                                                                                       aGCTTTTTTCCCGCCCCAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTaa                 <  2:1068246/150‑1 (MQ=255)                                                                                                                                          ttttttcCCGCCACAAAAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTg              <  2:728842/148‑1 (MQ=255)                                                                                                                                            ttttcCCGCCGCAAAAAATCCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGa             <  2:112824/149‑1 (MQ=255)                                                                                                                                                      gcAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGGGAGA‑GGAAAAAAAATAAAACCCCCATGTTGTGTGGGGATCGCTGCTGGGGAAAAAGAGGGAAGGGAAGGGGGGGGGAACAGCCCCCGCACCCCGCCCAATGAGTTAGCCTc    >  2:672824/1‑148 (MQ=255)                                                                                                                                                        aacaaAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCaa  <  1:148428/148‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                    CTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCAGACCTGTGTCAGCGGCTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCC‑ATGTGCAAAAATGCCTT‑CGCGGCGG‑CTGTCTGCCAGCTATTTTTCCGCCGCAACAAAACCGCCGTTCTCTCCAGTAGTGGC‑GGGGCAAGA‑GAAATAGCTTTAAGCCCGTCATGTTGTGTGGCAATCGCTGCTGGTAACAATGTGGAAAGGGAAGTGCGGCGAATCAGCTCCAGAACCGCGCTAATTGAGTTCGCCTCAA  >  NC_000913/357740‑358035

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