Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,223,143	0	A	C	40.0%	96.9 / 12.8	63	G162G (GGT→GGG)	iclR	transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (25/11);  new base C (0/24);  total (26/37)
Fisher's exact test for biased strand distribution p-value = 1.39e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCC  >  NC_000913/4222997‑4223290                                                                                                                                                       |                                                                                                                                                       aaTGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCg                                                                                                                                                    <  2:118295/151‑1 (MQ=255)   tGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGc                                                                                                                                                   >  1:24032/1‑150 (MQ=255)     aaTAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCAcc                                                                                                                                                <  2:666079/151‑1 (MQ=255)            ccGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCg                                                                                                                                        <  2:183900/151‑1 (MQ=255)            ccGTTTGCGCGTTTGGGCGAGATCTTCTTTGAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCg                                                                                                                                        <  2:664349/151‑1 (MQ=255)             cGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGt                                                                                                                                       <  1:585989/151‑1 (MQ=255)              gTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTg                                                                                                                                      >  2:313015/1‑151 (MQ=255)                  gcgcGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCAt                                                                                                                                   <  2:80822/150‑1 (MQ=255)                  gcgcGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAACAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCAt                                                                                                                                   <  1:569996/150‑1 (MQ=255)                   cgcgTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCg                                                                                                                                 >  1:691511/1‑151 (MQ=255)                    gcgTTTGGGCGAGATCTTCTTTTAATTGCACAGGAGACGCCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCgg                                                                                                                                <  2:516095/151‑1 (MQ=255)                    gcgTTTGGGCGAGATCTTCTTTTAAATGCACAGTAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCg                                                                                                                                 <  2:1170617/150‑1 (MQ=255)                         tGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAAtt                                                                                                                           >  2:339171/1‑151 (MQ=255)                         tGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAAt                                                                                                                            <  2:630154/150‑1 (MQ=255)                                  tcttctTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑ccgccg                                                                                                                 >  2:450393/1‑150 (MQ=255)                                     tcttTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGa                                                                                                                <  1:876616/148‑1 (MQ=255)                                      cttTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGAAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAgg                                                                                                            >  1:639519/1‑151 (MQ=255)                                          tAAATGCGCAGGAGACACCAGCGTTGCGGGGGTATAGGCATGTAATCCGTT‑GCGGAGCAACAGCTTCGTCACCTG‑TTCTCAGCTGAGTTAGGCTAAAAAGG‑CTTT‑GCCCGCAACGGAT‑GCGTGCATCTGCAATTT‑A‑CCGCCGATAGGcgcg                                                                                                        >  2:658467/1‑151 (MQ=255)                                           aaaTGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGcgcg                                                                                                        >  2:270726/1‑150 (MQ=255)                                             aTGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGAc                                                                                                     >  2:169362/1‑151 (MQ=255)                                                cacaGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACAt                                                                                                   >  2:368945/1‑150 (MQ=255)                                                        acacCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTA‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATc                                                                                            >  1:198923/1‑149 (MQ=255)                                                         cacCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCa                                                                                              >  1:539430/1‑146 (MQ=255)                                                           ccAGGGTGGGGTGGGTAAGGGCATGAACCCCTTG‑GGGGTCCAGCATTTTCGCCCCTTT‑TTCTTCCCTCATTTGGGCAAAAAAGT‑CTTC‑CCCCGCACGGAGA‑GCGTGAACCGCCATTTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATg                                                                                       <  1:241087/151‑1 (MQ=255)                                                             aGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GTGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGAAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCg                                                                                     >  1:981339/1‑151 (MQ=255)                                                              gCGTTGGGTGGGTATAGGCATGTACCCCTTT‑GCGGTGCAGCAGCTTCGTCCCCTG‑TTCTTCGCTCAGTGGGGCAAAAAAGG‑CTTT‑CCCCGCCCCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGt                                                                                    <  2:813140/151‑1 (MQ=255)                                                               gttttcGTGGGTGTAGGAATTTACCCTTTT‑GCGGGGCACGAACTCCGTCTCATT‑TTTTTCCCTCAGTTGGG‑TAAAAAGGTTTTT‑CCCCCCCCCGAAA‑GGGTGATGCGGAAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTa                                                                                   <  2:769335/146‑1 (MQ=255)                                                               cGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGt                                                                                    >  2:176467/1‑150 (MQ=255)                                                                   gcgggagTATAGGAATTTACCCCTTG‑GCGGTGCAGAGGTTTCGCCCCCTT‑TTCTCCGCCAAGTTGGTCTAAAAAGC‑CTTT‑CCCCCCCCGGAAA‑GCGTGCCCGGCCAATTT‑A‑CCGCCGCTAGGCGCGGACATTCGCATCAGATGCGTa                                                                                   <  1:452259/141‑1 (MQ=255)                                                                   gCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACt                                                                               >  1:860272/1‑151 (MQ=255)                                                                     gTGGTTATGGGCTTGTACCCTTTG‑GCGGTGCAGCAGTTCCTTCCCCTT‑TTTTCCCCTCATTGGGGATAAAAAGG‑CTTC‑CCCCGCCCGGGAG‑GCGTGCAGGGGGATTTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTg                                                                              <  1:739295/150‑1 (MQ=255)                                                                     gTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGt                                                                             >  1:944956/1‑151 (MQ=255)                                                                       gggTATAGGAATGTACCCTTTG‑GGGGTGCGGCGTCTCCGTCCCCTT‑TTCTTCGCTAATTGGGGCAAAAAAGG‑TTTT‑CCCCGCACGGGAG‑GCGTGCAGGGGAATTTT‑G‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTaca                                                                                 <  2:43171/145‑1 (MQ=255)                                                                          tataGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTg                                                                        >  2:625465/1‑151 (MQ=255)                                                                              ggCATGTACCCTTTT‑GCGGTGCAGCAGCTTCGTCCCCTG‑TTCTTCGCTCAGTTGGGCAAAAAAGG‑CTTT‑CCCCGCACCGGAG‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCg                                                                    <  2:960171/151‑1 (MQ=255)                                                                                 aTGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGat                                                                  >  1:965370/1‑150 (MQ=255)                                                                                  tGTAACCTTTT‑GGGGTGCAGCAGTTTCGCCCCCTG‑TTCTTCGCTCAGTTGGGCAAAAAAGG‑CTTT‑CCCCGCCCGGGAA‑GCGTGCATCGGAAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGata                                                                 <  2:483739/150‑1 (MQ=255)                                                                                      ccccTTT‑GCGGTGCAGCGGCTGCAAC‑CCTGTTTCTTCGCTTCGTGGGAAAAAAAAGG‑TTTT‑TCCCGCCCGAGAG‑GCGGCCCCCGGAATTTA‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGataata                                                              <  2:769052/148‑1 (MQ=255)                                                                                      aCCCTTTGGGGGTGCAGCAGCTTCGTCCCCTT‑TTCTCCGCTCAGTGGGGCAAAAAAGG‑TTTC‑CCCCGCCCGGGAG‑GCGTGCTTCGGAATTTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGata                                                                 <  2:1108274/147‑1 (MQ=255)                                                                                       cctttt‑GCGGTGCAGCAGCTTCGTCCCCTT‑TTCTTCCCCCAGTGGGCCAAAAAAGG‑TTTC‑CCCCCCCCCGAAG‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATaataat                                                            <  1:994498/147‑1 (MQ=255)                                                                                       cccTTG‑GGGGTGCAGCAGCTTCGTCCCCTG‑TTCTTCGCTCAGTGGGGCAAAAAGGG‑CTTT‑CCCCGCACCGGAG‑GCGTGCTCGGGCATTTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATaataat                                                            <  2:206584/150‑1 (MQ=255)                                                                                           ttggCGGTGCAGCGGTTTCGTCACCTG‑TTCTTCGCTC‑GTTTGGCCAAAAAGG‑GTTTGCCCCGCACCGGAG‑GCGTGCCCGGGCATTTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGc                                                         <  1:169356/147‑1 (MQ=255)                                                                                           tt‑GCGGCGCCGCAGCTGCGCCCCCTG‑TTCTCCCCTCAGGGGGGAAAAAAAGT‑CTTC‑CCCCGC‑CCGGAAGGCGGGCGCCGGAAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCg                                                          <  1:426272/148‑1 (MQ=255)                                                                                              ggggTGCGGCGGGTTCGTCCCCTG‑TTCTCCCCTCGTTTGGACAAAAAGGG‑CTTT‑CCCCGCACGGGAG‑GCGTGCACCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTc                                                      <  2:527716/148‑1 (MQ=255)                                                                                                ggTGCAGCAGCTTCGTCCCCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑CCCCGCCCCGGAA‑GCGTGCATCGGAAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTg                                                   <  2:18316/151‑1 (MQ=255)                                                                                                 gTGCAGCAGCTTCGTCCCCTT‑TTCTCCGCTCAGTGGGGCAAAAAGGG‑CTTT‑CCCCGCCCGGGAG‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGa                                                  <  2:723797/151‑1 (MQ=255)                                                                                                 gTGCAGCAGCTTCGTCCCCTG‑TTCTTCGCTAAGTTGGGCAAAAAAGG‑CTTT‑CCCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGa                                                  <  2:890709/151‑1 (MQ=255)                                                                                                       cagcTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCtt                                            >  1:935491/1‑151 (MQ=255)                                                                                                           tccgccACCTT‑TTCTCGGCCCGGTTGGGCAAAAAAGG‑CTTA‑GCCCCCACGGGAA‑GGTTGCACCGGAATTTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTg                                          <  1:186429/144‑1 (MQ=255)                                                                                                             cGTCCCTTG‑TTCTGCGCTCAGTTGGCCAAAAAGGG‑TTTT‑CCCCGCACCGGAG‑GCGTGCATCGGCATTTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGAt                                        <  1:279229/149‑1 (MQ=255)                                                                                                             cGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATc                                       >  1:1207549/1‑150 (MQ=255)                                                                                                                cACGTG‑TTCTTGTCTAATTTGGGCAAAAAAGC‑CTTT‑CCCCGCCCCGAAC‑GCGTGCTCCGCCATTTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCaa                                     <  2:1189788/149‑1 (MQ=255)                                                                                                                cACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGc                                   >  1:38748/1‑151 (MQ=255)                                                                                                                    tG‑TTCTTCGCTCGGTGGGGCAAAAAGGG‑CTTT‑CCCCGCCCGGGAG‑GCTTGCATCGGAAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCAcc                                <  2:295804/150‑1 (MQ=255)                                                                                                                       tcctccGCTCTGTTGGTCAAAAAGGG‑CTTC‑CCCCGCCCCGGAG‑GCTTCCTTCGCCATTTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGc                              <  2:1097854/145‑1 (MQ=255)                                                                                                                        tcttcGCTCAGTGGGGCAAAAAAGG‑CTTT‑CCCCGCACGGGAA‑GCGTGCATGGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGcc                             <  2:986856/150‑1 (MQ=255)                                                                                                                             gCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTg                       >  1:962201/1‑151 (MQ=255)                                                                                                                                 aGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCg                   >  2:263789/1‑151 (MQ=255)                                                                                                                                      gggaaaaaaGG‑TTTT‑CCCCCCCACGGAG‑GGGGGCATTGGCAATTT‑ACCCGCCGATAGGCGCGGACATCCGCATCAGATGCTTACCCTGTCCCTGGTCGATAATAATCGCTTCGTGATCGTTTTGATCAAGCACCGCCATATTGACCGTTTc               <  2:1016735/147‑1 (MQ=255)                                                                                                                                          aaaaaGG‑TTTT‑CCCCCCCCGGGGA‑GGGGTCCTCGGAAATTTCC‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTCCCTGGTCGATAATAATCGCTTCGTGATCGTTTTGATCAAGCACCGCCATATTGACCGTTTCGCCa           <  2:648703/151‑1 (MQ=255)                                                                                                                                           aaaaGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGcc            >  2:1009825/1‑148 (MQ=255)                                                                                                                                                  cTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTc   >  2:711201/1‑151 (MQ=255)                                                                                                                                                    tt‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTcc  >  2:126844/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                       AATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTT‑GCGGTGCAGCAGCTTCGTCACCTG‑TTCTTCGCTCAGTTGGGCTAAAAAGG‑CTTT‑ACCCGCACCGGAA‑GCGTGCATCGGCAATTT‑A‑CCGCCGATAGGCGCGGACATTCGCATCAGATGCGTACACTGTACCTGGTCGATAATAATCGCTTCGTGATCGCTTTGATCAAGCACCGCCATATTGACCGTTTCGCCAGACTCTTCC  >  NC_000913/4222997‑4223290

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