Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,376,505	0	A	C	39.4%	107.9 / 19.3	66	T158P (ACC→CCC)	ycjR	putative TIM alpha/beta barrel enzyme
Reads supporting (aligned to +/- strand):  ref base A (33/7);  new base C (0/26);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 5.16e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATT  >  NC_000913/1376360‑1376651                                                                                                                                                       |                                                                                                                                                   aCCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGAGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTc                                                                                                                                               >  1:500265/1‑151 (MQ=255)    gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGcc                                                                                                                                            <  1:1260108/151‑1 (MQ=255)     ccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATGATATGATCAACACCCTCGCCg                                                                                                                                           <  1:1210122/151‑1 (MQ=255)      cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CA‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCAGATGATCAACACCCTCGCCg                                                                                                                                           <  1:922395/150‑1 (MQ=255)       gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                         >  1:1162251/1‑151 (MQ=255)       gATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                         >  2:1228586/1‑151 (MQ=255)        aTGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                        >  2:1073789/1‑151 (MQ=255)         tGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCTTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAC‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGAt                                                                                                                                         <  2:1064352/149‑1 (MQ=255)         tGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATg                                                                                                                                        >  1:243742/1‑150 (MQ=255)               cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                                  <  2:790037/150‑1 (MQ=255)                 ccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                               >  1:683249/1‑151 (MQ=255)                  cgcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGtt                                                                                                                               >  2:146515/1‑150 (MQ=255)                                ggCGACCGCAAAATGGTGAGTGTTTCCCTGCGCGTACTGGAACGGGTCGCCGCGCGTCCGGGAA‑CC‑GTGGTGTATCTCGACCCGTT‑AAC‑CC‑GCTATCAGGATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAAc                                                                                                                 <  2:857048/150‑1 (MQ=255)                                  cgcccGAAAAAGGGTGAGTGTTTCCCTGCGCGTACGGAAACAGGTCGCCGCGCGTCCCGGAC‑CC‑GTGGTGTATCTC‑AACCCTT‑AAACCC‑GCTATCAGGATCATATGACCAACCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGAt                                                                                                              <  2:113455/148‑1 (MQ=255)                                   gACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                             >  2:1131332/1‑151 (MQ=255)                                   gACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATc                                                                                                             >  1:26842/1‑151 (MQ=255)                                    cccGAAAAAGGGTGAGTGTTCCCCTGCGTTTACGGAAACGGGGCCCCGCGCGTACCGGAC‑CG‑GGGGTTTATCTC‑AACCGTT‑AAACCC‑GCTACGGGGATAAAATAACCAACCCCCCCGCCGATCCCCGCCGTTACATCGTCGaaaa                                                                                                                  <  1:1203249/144‑1 (MQ=255)                                    aCCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATTAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGCTACATCGTCGAAAACGATCt                                                                                                            >  2:14235/1‑151 (MQ=255)                                         aaagaggTGAGTTATCCCCTGCGGTTGCGGAAACAGGTCGCGGCGCGAACCGGAC‑CG‑GTGGTGCCTCCCGACCCGTA‑AAC‑CC‑GCTATCGGGATAATATCACCAACCCCGGCGGCGAGCCCCGCCGTTACATCGTCGAAAACGATc                                                                                                             <  1:721841/140‑1 (MQ=255)                                              ggTGAGTGATCCCCTGCGCGTACTGAAACAGGCCGCCGCGCGTCCCGGAC‑CG‑GTGGTGTATCTCGACCCTTT‑AAC‑CC‑GCTATAAGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                    <  2:267720/149‑1 (MQ=255)                                                   gttctccccGGCGCGTACTGAAACAGGCCGCCCCGCCTCCCGGAC‑CC‑GTGGTTTCCCACGACCTTTT‑AAC‑CC‑TCTAACAGGATAATATAAACCCCCCCCCCGCAAATGTCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                             <  1:80979/145‑1 (MQ=255)                                                    tgtttCCCGGCGAGTACTGAAACGGGTCGCGGGGGGCACCGGAA‑CG‑GTGGTGTCTCTCGACCCGTT‑AAA‑CC‑CCTACCGGTCCCATTTCACCCAACCCCCCGCCGCTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAg                                                                                               <  2:21182/145‑1 (MQ=255)                                                    tgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                            >  1:970213/1‑151 (MQ=255)                                                     ttttCCCGGCGGGTACGGGAACAGGTCGCCGCGCGTACCGGAA‑CG‑GTGTTGTATCCCGACCCGTT‑AAC‑CC‑GCTATGAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                             <  2:377211/147‑1 (MQ=255)                                                       ttCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                         >  1:519040/1‑151 (MQ=255)                                                       tccccGGCTCGTACGGAAACGGGTCCCCGCGCGCACCGAAC‑CG‑GTGGTGTATCAGGCACCGTT‑AAA‑CC‑GTTAAAAGGATAATATGACCACCCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACa                                                                                                <  1:1232233/142‑1 (MQ=255)                                                             gcgcGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAt                                                                                   >  1:939366/1‑151 (MQ=255)                                                                 gTACTGAAACGGGTCGCCGCGC‑TACCGGAAGCC‑GTGGT‑TTTCCCCAACCCTT‑AAA‑CCTCCTATCAGGTTAAAATAATAACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGAtt                                                                                  <  1:840544/148‑1 (MQ=255)                                                                     gggACCAGGCCGCCGCGCGAACGGGAA‑CG‑TGGGTGTATCTCAACCCTTT‑AAC‑CC‑TCATTCAGAACAATATGATCAAAACCCTCCCCTGTTCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTa                                                                              <  1:90393/147‑1 (MQ=255)                                                                     gggAACGGGCCCCCGCGCGCCCCGACC‑CG‑GTTGTGTATCTCGAACCGTA‑AAC‑CC‑CATATCAGGATATTATCATCACCCCCCTCCCGGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                            <  1:604502/149‑1 (MQ=255)                                                                        aaCAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACAt                                                                        >  2:175318/1‑151 (MQ=255)                                                                           aGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                      >  1:667348/1‑150 (MQ=255)                                                                            ggcccccgcGCGTACCGGAC‑CG‑GTGATGTATCCCAAACCTTT‑AAC‑CC‑TCTATAGGAATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                         <  1:425396/141‑1 (MQ=255)                                                                            ggTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAat                                                                    >  1:445651/1‑151 (MQ=255)                                                                             ggccgcggggcgcccGGAA‑CCGGTGGTGTTCCTCAAACCTTA‑AAC‑CC‑GTTATAAGGATAATATAAACACCCCCCGCGCCGTGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                      <  1:1030904/136‑1 (MQ=255)                                                                               cgccgcGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGCCGATTTCTATCACATGAATATc                                                                 >  1:562336/1‑151 (MQ=255)                                                                                gccgcGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCg                                                                >  2:1086167/1‑151 (MQ=255)                                                                                  cgcgcgCACCGAAC‑CG‑GTGGTGTACCTCGACCCGTA‑AAA‑CC‑GCAATAAGGATCATATGACAACCCCCCGCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                              <  2:1162316/151‑1 (MQ=255)                                                                                   gcgcgTCCCGGAA‑CG‑GTGGTTTA‑CCCCAACCTTTAAAA‑CC‑TTTATCAGTATAAAATGATAAACCCCCCCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                               <  1:330997/149‑1 (MQ=255)                                                                                        cggccgggaccgggggtaTATCCCGAACGGTA‑AAA‑CC‑CCTATCAGGTTCATATAATCAACCCCCTCGCAGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                            <  2:23759/135‑1 (MQ=255)                                                                                          ccGGAC‑CG‑GTGGTGTATCTCGAACCGTT‑AAC‑CC‑GCTATCAGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                       <  2:238218/150‑1 (MQ=255)                                                                                          ccGGAA‑CT‑GTGTTTTACCGCCAACCGTA‑AAC‑CC‑GTTACAAGGATAATATGACCAACCCCCGCCCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaa                                                              <  2:294394/143‑1 (MQ=255)                                                                                           ccgaaccc‑GTGGTGTATCTCGAACCGTT‑AAC‑CC‑GCTATCAGGATCATATAAAAAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATaa                                                      <  2:1304558/145‑1 (MQ=255)                                                                                                accggtggttTATCTCAAACCGTA‑AAC‑CC‑GCTATCAGGATAATATGATAAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                     <  2:644371/144‑1 (MQ=255)                                                                                                     tggtgTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                              >  1:379000/1‑150 (MQ=255)                                                                                                      ggtgTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                              >  2:1201394/1‑149 (MQ=255)                                                                                                        tgtATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                           >  1:206311/1‑150 (MQ=255)                                                                                                            tcccaaaCCTTA‑AAC‑CC‑GTTAACAGGATAAAATAATCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                          <  1:660972/142‑1 (MQ=255)                                                                                                                gAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCAt                                   >  1:860727/1‑150 (MQ=255)                                                                                                                 aaCCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                                  >  2:1059465/1‑150 (MQ=255)                                                                                                                   ccGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                                >  1:575468/1‑150 (MQ=255)                                                                                                                    cGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                               >  2:177519/1‑150 (MQ=255)                                                                                                                      tt‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAc                             >  2:843231/1‑150 (MQ=255)                                                                                                                      cggtaaaccccatATCAGGATCATATGACAAACCCCCTCGCAGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACa                               <  2:528758/139‑1 (MQ=255)                                                                                                                       t‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACcg                           >  1:350508/1‑151 (MQ=255)                                                                                                                                gCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCt                     >  2:313285/1‑151 (MQ=255)                                                                                                                                     cAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg                >  1:1271045/1‑151 (MQ=255)                                                                                                                                       gtctaatatTATCAGCCCCCGCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg               <  1:1173046/145‑1 (MQ=255)                                                                                                                                       ggATCATATGACCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCgg               <  2:1239661/150‑1 (MQ=255)                                                                                                                                       gaaTCATAATATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGt              <  2:1277961/149‑1 (MQ=255)                                                                                                                                        aataataaaaccAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCt                  <  1:775067/135‑1 (MQ=255)                                                                                                                                          tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt           >  2:1300911/1‑151 (MQ=255)                                                                                                                                            atatGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg          >  1:290212/1‑150 (MQ=255)                                                                                                                                              aggaccAACACCCGCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtgtg        <  1:500675/145‑1 (MQ=255)                                                                                                                                                gATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat     >  2:345939/1‑151 (MQ=255)                                                                                                                                                   cAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATAtt  >  1:1000310/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                   ACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAA‑CC‑GTGGTGTATCTCGAACCGTT‑AAA‑CC‑GCTATCAGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCATATT  >  NC_000913/1376360‑1376651

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