Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,638,193	(C)6→7	coding (16/1353 nt)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,638,187	1	.	C	98.7%	273.2 / NA	77	S4P (TCC→CCC)	creD	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  major base C (33/42);  minor base G (0/1);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAAT‑CCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGG  >  NC_000913/4638039‑4638329                                                                                                                                                        |                                                                                                                                                gtttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAAT‑c                                                                                                                                                >  1:747904/1‑150 (MQ=255) gtttAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAAT‑c                                                                                                                                                >  2:78791/1‑150 (MQ=255)              gTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTgg                                                                                                                                  >  1:14556/1‑151 (MQ=255)                cACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGaa                                                                                                                                <  2:367115/151‑1 (MQ=255)                cACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGaa                                                                                                                                <  2:437586/151‑1 (MQ=255)                     tGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTGCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGa                                                                                                                           <  2:785365/151‑1 (MQ=255)                           aCGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGc                                                                                                                      >  1:376106/1‑150 (MQ=255)                           aCGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGc                                                                                                                      >  2:348088/1‑150 (MQ=255)                             gTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGc                                                                                                                      >  1:652939/1‑148 (MQ=255)                             gTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGc                                                                                                                      >  1:996026/1‑148 (MQ=255)                             gTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGc                                                                                                                      >  1:427569/1‑148 (MQ=255)                                   gAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGttt                                                                                                                <  2:791225/148‑1 (MQ=255)                                      ggTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCa                                                                                                          >  2:249516/1‑151 (MQ=255)                                       gTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                                                                           <  2:865867/149‑1 (MQ=255)                                       gTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                                                                           <  2:127732/149‑1 (MQ=255)                                       gTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                                                                           <  1:1188575/149‑1 (MQ=255)                                       gTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                                                                           <  2:7558/149‑1 (MQ=255)                                       gTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                                                                           <  1:920991/149‑1 (MQ=255)                                       gTGGCGTGCTGGCCTCGCTTCGACTGCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGc                                                                                                           <  2:457415/149‑1 (MQ=255)                                           cGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                                                                      <  1:200715/150‑1 (MQ=255)                                           cGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                                                                      <  2:529743/150‑1 (MQ=255)                                           cGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                                                                      <  2:259700/150‑1 (MQ=255)                                           cGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                                                                      <  1:869434/150‑1 (MQ=255)                                           cGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                                                                      <  1:507667/150‑1 (MQ=255)                                           cGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                                                                      <  1:920040/150‑1 (MQ=255)                                           cGTGCTAGCCTCGCTTCGTCTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTAt                                                                                                      <  2:1047696/150‑1 (MQ=255)                                            gTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                                                                    >  2:493524/1‑151 (MQ=255)                                            gTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTg                                                                                                    >  2:968757/1‑151 (MQ=255)                                             tGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                                                                   <  1:209141/151‑1 (MQ=255)                                               cTGGCCTCGCTTCTGCTTCACCGTCACTTCGCATAGCTTCAAATACGTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGGAGACTATGTTGAAATGCCCC‑CCTGGTCTGAAAAATGCCTAGCCTGTTTGGTGCAGTATTGCtg                                                                                                 <  2:406926/151‑1 (MQ=255)                                               cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                                                                 <  1:651391/151‑1 (MQ=255)                                               cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                                                                 <  1:99795/151‑1 (MQ=255)                                               cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                                                                 <  2:871113/151‑1 (MQ=255)                                               cTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCtg                                                                                                 <  1:618511/151‑1 (MQ=255)                                                     tCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                                                                           <  1:19512/151‑1 (MQ=255)                                                     tCGCTTCGACGTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGGCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGttgttg                                                                                           <  2:440099/151‑1 (MQ=255)                                                        cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                                                                         <  1:122002/150‑1 (MQ=255)                                                        cTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGAt                                                                                         <  1:213954/150‑1 (MQ=255)                                                                 cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                                                                               >  1:195434/1‑151 (MQ=255)                                                                 cACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATg                                                                               >  1:40326/1‑151 (MQ=255)                                                                  aCCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGc                                                                              <  1:452434/151‑1 (MQ=255)                                                                       cACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                                                                         <  2:516917/151‑1 (MQ=255)                                                                        aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                                                                         >  1:723220/1‑150 (MQ=255)                                                                        aCTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGAtt                                                                         >  1:283408/1‑150 (MQ=255)                                                                         cTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCg                                                                       >  1:1135889/1‑151 (MQ=255)                                                                             acaTAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAg                                                                   >  2:396920/1‑151 (MQ=255)                                                                                  gCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGAt                                                              >  2:734625/1‑151 (MQ=255)                                                                                     tCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGt                                                           <  1:638468/151‑1 (MQ=255)                                                                                      cAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTc                                                          <  2:530995/151‑1 (MQ=255)                                                                                          ttcttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAAc                                                      <  2:740018/151‑1 (MQ=255)                                                                                            cttcCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGt                                                    >  1:502967/1‑151 (MQ=255)                                                                                               cccACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGGGATTGTCGAACGTg                                                   >  2:409407/1‑149 (MQ=255)                                                                                                 cacaTAGTCTTCGTATCCTGCTGCCATTGCAAAGGAGAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTg                                                   >  1:1072677/1‑149 (MQ=255)                                                                                                     tAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTAc                                           <  1:828059/151‑1 (MQ=255)                                                                                                       gTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCg                                         >  2:294576/1‑151 (MQ=255)                                                                                                       gTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCg                                         >  2:1073774/1‑151 (MQ=255)                                                                                                        tCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGt                                        <  1:274038/151‑1 (MQ=255)                                                                                                        tCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGt                                        <  1:1065744/151‑1 (MQ=255)                                                                                                          ttCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAg                                      <  2:1066830/151‑1 (MQ=255)                                                                                                              tATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGa                                   <  2:747588/150‑1 (MQ=255)                                                                                                              tATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGa                                   <  2:848948/150‑1 (MQ=255)                                                                                                              tATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGa                                   <  1:748827/150‑1 (MQ=255)                                                                                                               aTCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                                                                   <  1:1181495/85‑1 (MQ=255)                                                                                                               aTCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGc                                                                                                   >  2:1181022/1‑85 (MQ=255)                                                                                                               aTCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGAtg                                 >  2:1121882/1‑151 (MQ=255)                                                                                                                     ctgcCATTGCAAAG‑‑GAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAg                           <  1:52795/151‑1 (MQ=255)                                                                                                                     ctgcCATTGCAAAG‑‑GAGAAGACTATGTGGAAAG‑CCCCGCCTGTTCTGGAAAATGACTAGACTGTTTGGGGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTGGAACGTGCTGATTACCGTAGCGATGTGGAAg                           <  2:630229/151‑1 (MQ=255)                                                                                                                                  g‑‑gAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCa              >  1:1061213/1‑151 (MQ=255)                                                                                                                                  g‑‑gAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCa              >  1:223589/1‑151 (MQ=255)                                                                                                                                  g‑‑gAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCa              >  1:731887/1‑151 (MQ=255)                                                                                                                                  g‑‑gAGAAGACTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCa              >  1:731038/1‑151 (MQ=255)                                                                                                                                         agaCTATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTa         >  1:562432/1‑151 (MQ=255)                                                                                                                                             tATGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAg     <  1:1000136/151‑1 (MQ=255)                                                                                                                                              aTGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTCGCGATGTGGAAGATGCGATTCGCCAAAGTACCAg     >  2:807667/1‑150 (MQ=255)                                                                                                                                              aTGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAg     >  2:472397/1‑150 (MQ=255)                                                                                                                                              aTGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAg     >  1:990558/1‑150 (MQ=255)                                                                                                                                              aTGTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAg     >  1:771969/1‑150 (MQ=255)                                                                                                                                                gTTGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCgg  >  1:240177/1‑151 (MQ=255)                                                                                                                                                 ttGAAATCCCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCgg  <  2:1056848/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                GTTTAACGGCGAAGTCACGCTGCGCAACGTGCAGGAAGGTGGCGTGCTGGCCTCGCTTCGACTTCACCGTCACTTCACATAGCTTCAAATTCTTCCCACATAGTCTTCGTATCCTGCTGCCATTGCAAAG‑‑GAGAAGACTATGTTGAAAT‑CCCC‑CCTGTTCTGGAAAATGACTAGCCTGTTTGGTGCAGTATTGCTGTTGTTGATTCCGATAATGCTGATTCGGCAGGTGATTGTCGAACGTGCTGATTACCGTAGCGATGTGGAAGATGCGATTCGCCAAAGTACCAGCGG  >  NC_000913/4638039‑4638329

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