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

CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376355‑1376648                                                                                                                                                           |                                                                                                                                                cGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAaca                                                                                                                                                 <  2:1070741/151‑1 (MQ=255)    ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAaca                                                                                                                                                 <  2:811400/148‑1 (MQ=255)    ttACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACAccc                                                                                                                                              >  2:834795/1‑151 (MQ=255)     tACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCt                                                                                                                                             >  1:73419/1‑151 (MQ=255)        cgccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GCGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCg                                                                                                                                           >  2:690495/1‑150 (MQ=255)         gccgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGcc                                                                                                                                         <  2:499771/151‑1 (MQ=255)           cgATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCg                                                                                                                                        >  1:112386/1‑150 (MQ=255)               gACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGcc                                                                                                                                   >  1:321781/1‑151 (MQ=255)                    cgccgcGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCccgcc                                                                                                                               <  2:177814/150‑1 (MQ=255)                        gcgtggcgggggcGGCGACCGAAAAGTGGTGATTGCTCCCCTGCGCGTACGGAAACGGGCCGCCGCGCGTACGGGCCCG‑GTGGTGTA‑CCCCAACCT‑TTAAACCCGTTATC‑AGAATCATATGACCACCCCCCTCGCCGATGCCCGCCGtt                                                                                                                            <  1:84560/137‑1 (MQ=255)                        gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTAc                                                                                                                          >  1:108209/1‑151 (MQ=255)                        gcgTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCAGATGCCCGCCGTTAc                                                                                                                          >  2:490623/1‑151 (MQ=255)                           tAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACAtc                                                                                                                       >  1:371094/1‑151 (MQ=255)                              cgtGGACGGCGACCAAAAATTGTGGATTGATCCCCTGCGCTTACTGAAACGGGCCCCCGCCCGTCCGGGAGCC‑GTGTTTTTTCTCGACCTT‑TTAAC‑CCGCTAGC‑AGAATAATATGATCACCCCCCCCGCCGATGCCCGCCGTTACAtcg                                                                                                                      <  1:953437/147‑1 (MQ=255)                              ccgggACGGCGCCCGAAAAATGGTGAGTGATCCCCGGCGCGTACTGGAACAGGTCGCCGCGCGTCCCGAACCC‑GTGGTGTATCTCGACCCG‑TAAAC‑CCGCTATC‑AGAATAATATGATAAACCCCCTCGCCGATGCCCGCCGTTACAtcgtc                                                                                                                    <  2:91797/148‑1 (MQ=255)                                  gACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGaa                                                                                                                 >  2:527551/1‑150 (MQ=255)                                     gccgcccaaaaaaGGGTGGGTATTTCCCGGCGCTTCCGGGAACGGGCCGCCGCGCGTACCGAAACG‑GTGGTTTATC‑CCAACCC‑TTAAACCCGTTACA‑AGAATAAAATGACAACCCCCCCCGCGGAGGCCCGCCGTTACATCGTCGAAAAc                                                                                                              <  1:638694/141‑1 (MQ=255)                                      gCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACg                                                                                                             >  2:951702/1‑150 (MQ=255)                                          ccaaaaaaGGGTGAGTTTTCCCCGGCGGGTACGGGAAGAGCCCCCCGCGCGTCCCAGAACT‑GGGGTTTATCTGAAAGAG‑TAAAA‑CCGCTATC‑AGAATCATATGACCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGa                                                                                                            <  2:407593/143‑1 (MQ=255)                                          ccGAAAAAGGGTGGGTGATCCCCTGCGCGTACGGGAACGGGCCGCCGCGCGTACCGGACCG‑GTGGTGTATCCCAACCCG‑TTAAC‑CCGCTATC‑AGGATAATAGGACAACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGa                                                                                                            <  1:1082458/147‑1 (MQ=255)                                             cAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTaa                                                                                                      >  2:1016869/1‑150 (MQ=255)                                                aaTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACAt                                                                                                  >  1:290833/1‑151 (MQ=255)                                                 aTGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATg                                                                                                 >  2:815217/1‑151 (MQ=255)                                                  tGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGt                                                                                                >  2:206770/1‑151 (MQ=255)                                                    gggagtgaTCCCCTGCGCGTACTGAAACAGGCCGCCGCGCGTACCGGACCC‑GTGGTGTACCTCGACCCT‑TAAAA‑CCGTTATA‑AGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTAc                                                                                              <  1:1274732/149‑1 (MQ=255)                                                         ttatccccGGCGCGTACTGGAACAGGTCCCCGCGCGTCCCGGACCC‑GTGGTTTATCCCAACCCG‑TAAAC‑CCGCTATC‑AGGATAATATAATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAt                                                                                          <  1:204699/145‑1 (MQ=255)                                                         tgaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                         >  2:710092/1‑151 (MQ=255)                                                          gaTTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                         >  1:215949/1‑150 (MQ=255)                                                            ttCCCTGCGCGTACTGGACCAGGCCGCCGCGCGTCCCGGACCC‑GGGGTGTATCTCGACCCG‑TTAAA‑CCCCTATC‑AGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                         <  2:460271/148‑1 (MQ=255)                                                            ttCCCTGCGCGTACTGGAACGGGTCGCCGCGCGTCCCGGACCG‑GTGGTGTATCTCGACCCT‑TTAAC‑CCGCTACC‑AGGATCATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTAt                                                                                       <  2:606885/150‑1 (MQ=255)                                                            tccccGGCGCGTACGGGAACAGGCCGCCGCGCGTCCGGGACCG‑GTGGTTTATCTCGACCCG‑TAAAC‑CCGCTATC‑AGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGAtt                                                                                         <  2:1168672/146‑1 (MQ=255)                                                              cccGGCGCGTACGGGAACAGGTCGCCGCGCGTCCGGAAACG‑GTGGTGTATCTCGACCCT‑TTAAC‑CCGCTATC‑AGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATc                                                                                      <  2:977305/149‑1 (MQ=255)                                                               cttgggcttacggaaaCAGGCCGCCGCGCGCA‑CGGAACCGGTGGTGTACCTC‑AACCC‑TTAAA‑CCCCCAACAAGGATAAAATAACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAg                                                                                            <  1:840805/128‑1 (MQ=255)                                                               ccggcgcGTACTGAAACGGGTCGCCGCGCGTCCCGGACCG‑GTGGTGTATCCCGACCCT‑TTAAC‑CCGCTATC‑AGGATAATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                   <  1:1096036/148‑1 (MQ=255)                                                               ccTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGc                                                                                   >  1:165832/1‑151 (MQ=255)                                                                          tGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         >  1:1177705/1‑150 (MQ=255)                                                                             aaCAGGTCGCCGCGCGTCCCGAACCG‑GGGGTGTATCTC‑AACCGTTTAAC‑CCGCTATA‑AGGATAATATGACAACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATc                                                                         <  2:1220236/147‑1 (MQ=255)                                                                             aaCAGGCCCCCGCGCGTCCCGAACCG‑GTGTTTTTTCTCGACCCT‑TAAAC‑CCGCTATC‑AGGATAATATGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                      <  2:544396/150‑1 (MQ=255)                                                                                 ggcccccgcGCGAACCGGACCC‑GTGGTTTATCTCGACCCG‑TTAAC‑CCCTTATC‑AGGATAATATGCCCAACCCCCTCGCCGATCCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATcaca                                                                      <  1:1342255/141‑1 (MQ=255)                                                                                 ggTCGCCGCGCGTCCCGGACCG‑GTGGTGTATCTCGACCCT‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGa                                                                   <  2:1304724/149‑1 (MQ=255)                                                                                    cgccgcGCGTACCGGACCC‑GTGGTGTATCTCGAACCT‑TAAAC‑CCGCTATC‑AGGATAATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAatat                                                               <  2:496906/150‑1 (MQ=255)                                                                                    ccccgcGCG‑ACCCGAACCTGTGGTTTTCCTCAACCCG‑TTAAC‑CCGTTTTG‑AGGATCATTTGATCAACCCCCTCCCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAata                                                                <  1:1155472/147‑1 (MQ=255)                                                                                      ccgcGCGTCCCGGAACG‑GTGGTGTCTCTCAACCCT‑TTAAC‑CCGCTACA‑AGGATAATATGATAACCCCCCTCGCCATTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATc                                                              <  1:421640/149‑1 (MQ=255)                                                                                       cgcgcgTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCga                                                            >  1:868373/1‑150 (MQ=255)                                                                                         cgcgTCCGGGACCG‑GTGGTGTATCTCG‑ACCT‑TTAAACCCGCTATC‑AGGATCATATGATCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaag                                                          <  2:793153/150‑1 (MQ=255)                                                                                             tACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         >  1:298915/1‑147 (MQ=255)                                                                                             tACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCgaaga                                                         >  1:703442/1‑147 (MQ=255)                                                                                               ccGAAACC‑GTGTTGTATCTCGACCCT‑TTAAC‑CCGCTATC‑AGGATCATATAATAAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATa                                                    <  1:970437/150‑1 (MQ=255)                                                                                                  ggaccggtggtgTATCTCGACCCG‑TTAAC‑CCGCTATC‑AGTATCATATGATCAACCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAcc                                                 <  2:1064783/145‑1 (MQ=255)                                                                                                  ggaactgggtggtATCCCAACCCG‑TTAAC‑CCGCTATC‑AGGATTATATGTCCAACACCCTCGCGGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAAc                                                  <  2:460136/139‑1 (MQ=255)                                                                                                  gAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCt                                                >  2:541758/1‑151 (MQ=255)                                                                                                        gtggtgTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGgcgc                                           >  1:835344/1‑151 (MQ=255)                                                                                                            tgtATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                       >  1:210493/1‑151 (MQ=255)                                                                                                            tgtATCTCGAACCG‑TAAAC‑CCGCTATC‑AGGATCATAGGACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAgg                                        <  2:104831/150‑1 (MQ=255)                                                                                                             gtATCTCGACCCG‑TAAAC‑CCCCTATC‑AGGATAATATGATCAACCCCCCCGCCGTTGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGgc                                       <  2:148134/150‑1 (MQ=255)                                                                                                                     aaCTG‑TTAAC‑CCATTATC‑AGTATAATATAACCACCCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATg                               <  2:195384/150‑1 (MQ=255)                                                                                                                      acccttaaaa‑CCCCTATC‑AGGATCATATGACCAACCCCCTCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCa                                 <  2:215952/141‑1 (MQ=255)                                                                                                                      aCCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGAc                             >  2:264053/1‑151 (MQ=255)                                                                                                                              aacccGCTATC‑AGGATCATATAACCAACCCCCGCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAAcc                         <  1:613028/146‑1 (MQ=255)                                                                                                                                      aTC‑AGAATAATATAACCACCCCCCCCGCCGAGGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                 <  1:831547/149‑1 (MQ=255)                                                                                                                                        c‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg             >  2:814619/1‑151 (MQ=255)                                                                                                                                        c‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg             >  2:983871/1‑151 (MQ=255)                                                                                                                                           ggATCATATGACCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTg                 <  2:1197447/145‑1 (MQ=255)                                                                                                                                           ggATAATATAACCACCCCCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCg             <  2:211264/149‑1 (MQ=255)                                                                                                                                              tCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAt        >  1:1003199/1‑151 (MQ=255)                                                                                                                                                 tatGATCAACCCCCCCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTc              <  1:206767/142‑1 (MQ=255)                                                                                                                                                    gATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat  >  2:250281/1‑151 (MQ=255)                                                                                                                                                     aTCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCat  >  1:830528/1‑150 (MQ=255)                                                                                                                                                      tCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCAtg       >  2:355463/1‑144 (MQ=255)                                                                                                                                                           |                                                                                                                                                CGCTTACCGCCGATGACCTCGCCGCGTAGCCTGGACGGCGACCGCAAAATGGTGAGTGATTCCCTGCGCGTACTGGAACAGGTCGCCGCGCGTACCGGAACC‑GTGGTGTATCTCGAACCG‑TTAAA‑CCGCTATC‑AGGATCATATGATCAACACCCTCGCCGATGCCCGCCGTTACATCGTCGAAAACGATCTTAAACATGTACAGATTATCGGCGATTTCTATCACATGAATATCGAAGAAGATAACCTGGCGCAGGCGCTGCATGACAACCGCGACCTGCTCGGTCATGTGCAT  >  NC_000913/1376355‑1376648

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