Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,188,861	C→T	V306I (GTA→ATA)	phoQ ←	sensor histidine kinase PhoQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,188,861	0	C	T	95.0%	110.6 / NA	40	V306I (GTA→ATA)	phoQ	sensor histidine kinase PhoQ
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (25/13);  total (25/15)
Fisher's exact test for biased strand distribution p-value = 1.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCACGTTGCCCATCACCTCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTACCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGTTATGTTCTTCCAGTTCG  >  NC_000913/1188617‑1189090                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       aGCACGTGGCCCAACACCTCGACAAAATCGTTATGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGCGAGATATTGACCCCTTTGCGTTTATACACTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCGGAGCAATGTCCCGCCGCGCATACGGGCACGATGCAGG‑TAGTATCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCCGCCTTATCGGCTc                                                                                                                                                                                                                                  >  2:346438/1‑251 (MQ=255)    aCGTTGCCCATCACCTCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGCGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGc                                                                                                                                                                                                                               >  1:540180/1‑251 (MQ=255)          cccagctccTCTACAACATCGGTCTTCTCACCGACAAAGCTGATCGCTGGCGAAATATCGAGAGAGATATTGACCTCTTTGCGTTGATACACTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACTATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATctct                                                                                                                                                                                                                         <  2:431507/244‑4 (MQ=255)          cccATCACCTCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑t                                                                                                                                                                                                                        >  2:216992/1‑251 (MQ=255)               cACCTCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGAGTTGATACACTTTGTTCA‑GAGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGCTGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGc                                                                                                                                                                                                                   >  1:62720/1‑251 (MQ=255)                cccTCGACAAAATCGTTCCGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCt                                                                                                                                                                                                                  <  1:113618/250‑1 (MQ=255)                       cAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCttt                                                                                                                                                                                                           >  1:429747/1‑251 (MQ=255)                                        ccGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAAGGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCTAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACg                                                                                                                                                                                          <  1:222846/251‑1 (MQ=255)                                                        ctGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTGTTCACTACGCAGAGAACGCAGCGTACTTTgca                                                                                                                                                                            <  1:491380/249‑1 (MQ=255)                                                          ggCGCAATATCGTGAGAGATATTGACCCCTTTGCGTTCATACACCTTGTTCA‑GCGCTGAGGTGAGCTTGTCCAGCCGTGGAGCGCCCGGCTGC‑AGCTCGCGGCTGCGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGCAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATAAC‑TGACGCTCATCTTTTCACTCCGCCGAGAACGCAGCGTACTTTGAAGCAc                                                                                                                                                                        >  2:577675/1‑251 (MQ=255)                                                           gCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCA‑GCGCTGAGGTGAGAGTGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGTCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCAcc                                                                                                                                                                       <  2:578019/251‑1 (MQ=255)                                                                      ggagtgatatgcaCCCTTTTGCGTTGATACCCTTTGTCCA‑CCGCTGGGGTGAAGTTCTCCAGCATGGGGGCGACCGGATGC‑AGCTCGCGGCTGACCAAGTTCCCGCCCCTCATTCTGCCACGAGGCAGG‑TATTAGCCAATTTG‑CGGTAAAA‑GGCGCCTGATTTGCTCCAGCATTTTCGGCCCAGCATCCC‑GGCCGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGt                                                                                                                                                            <  2:92859/239‑1 (MQ=255)                                                                            atatTGACCCCTTTGCGTTGATACCCTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGCAGGGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCCCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAg                                                                                                                                                      <  1:216992/251‑1 (MQ=255)                                                                                        ttGCGTTGATACACTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATCCTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTc                                                                                                                                            <  1:392872/249‑1 (MQ=255)                                                                                         tggttggATACACTTTGT‑CACGCGCTGGGGTGGGATTGTCCAGCGGTGGGGCGCAGGGAGCC‑AGCCCGCGGCGGACCAATTCCCCGCCGCGCATACGGGCACGATGCAGG‑GAGTAGCCAATTTG‑CTGTAAAA‑GGCGGCTGATTTGCTCAAGCATTATCGGCTCAGCATCCC‑TGACGCTCATCTTTTCCCTACGCAAAGAACGCAGCGTACTGTGCAGAACCGCCAGTGGCGTTTTCAGACTATGGGTCAg                                                                                                                                          <  2:101591/244‑1 (MQ=255)                                                                                          gctttGATACACTTTGTCAA‑GCGGTGAGGTGAGATTGCCCAGCAGGGGGGCGGCCG‑ATGGGAGCTCGCGGCGGAGCAGGGTCCCGCCGGGCTGACTGGCAGGATG‑AGGATAGTA‑TCAATTTGGCTGTGAAA‑GGCGGGTGATTTACTCCAGCATTATCGGCTCAGCACCAC‑TGACGCAGATGTTTGCACTCCGCAGAGAGCGCAACGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGt                                                                                                                                        <  1:357712/248‑1 (MQ=255)                                                                                             ttGATACACTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCgg                                                                                                                                     >  2:88402/1‑251 (MQ=255)                                                                                                 tACACTTTGTTCA‑GCGATGAGGAGAGATTGTCCAGCAGGAGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCGGCCAAGCATACTGGCACGATGCAGA‑TAGGAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCGAGCATTATCGGCTCAGCA‑AACATGACGCTCATCTTTTAACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGACAGCGGCGTTGTCAGACACTGGGGCAGGTCgggg                                                                                                                                   >  1:544793/1‑247 (MQ=255)                                                                                                  acacTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGAAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCCGGTCGGTGAGc                                                                                                                                >  1:224223/1‑251 (MQ=255)                                                                                                                              tGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATGGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGACCGCGGCGTCCTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTc                                                                                                      <  1:214885/250‑1 (MQ=255)                                                                                                                                tCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTcgcg                                                                                                   >  1:288592/1‑251 (MQ=255)                                                                                                                                 ccAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGCGAGCGTCGTACGGTATTTGTCGTACCGTTCGCGt                                                                                                  >  1:18451/1‑251 (MQ=255)                                                                                                                                   agcagTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTc                                                                                                >  1:340553/1‑251 (MQ=255)                                                                                                                                           ggCGCCCGGAGGC‑AGCCCGCGGCGGAGCAAGGTCCCGCCGCGAATACTGGCACGAGGCAGT‑TAGTAGCAAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCCCCAGCATTACCGGCTCAGCATCCC‑TGACGTTCATTTTTTCACTACGCAGAGAACCCGGCGTCCTTGGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTa                                                                                        <  1:551119/251‑1 (MQ=255)                                                                                                                                              gACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGACGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAC‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCCACGCAAGTGGCGTTTTAAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAAc                                                                                      >  2:234527/1‑250 (MQ=255)                                                                                                                                                            ccggggCTGAGCAATTCCCCGCCGGTTATACTGGCACGATGGAGG‑TAGTAGCCATTTTG‑CT‑TGAAAATGCGTCTGATTTGCTCCGTCATTACCGGCTAAGCATAAC‑GGACGCTCATTTTTTCCATACGAAAAGAACGCGGCGTACTGTGCGGCCCCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGGCGTAACGTTCGCGTTCACTGTTTAACAATCGGTTCAGGtt                                                                        <  2:546695/247‑1 (MQ=255)                                                                                                                                                                           gTCCCGCCGAGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGt                                                         >  1:191076/1‑251 (MQ=255)                                                                                                                                                                                  cgcgCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTtct                                                  >  2:459255/1‑251 (MQ=255)                                                                                                                                                                                    cgcATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTcg                                                >  1:412324/1‑251 (MQ=255)                                                                                                                                                                                            ggCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGAATCAA‑TGACGCTCATATTTTCACTACGAAGAGAACGCAGCGTACTTTGCAGCACCGCCAGGGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTAAGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTgg                                        >  1:339906/1‑251 (MQ=255)                                                                                                                                                                                                 gATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTgg                                    >  1:94846/1‑250 (MQ=255)                                                                                                                                                                                                 gATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATGTGCTACAGCATTATCGCCTCCGCATCAC‑TGACGCTCCTCTTTTCACTACGCCGAGAACGCAGCGTACTTTGCCGCACCGCCAGTGGAGTTTTCAGACTATGGGTCAGGTCGGTGCGCGTCGTACGGAATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTggg                                   >  2:567517/1‑250 (MQ=255)                                                                                                                                                                                                              agCCAATTTG‑CTGTGAAA‑TGCGGCTGTTTTGCTCCAGAATTATCGGCTCAGCATCCC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGGGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTcg                       <  1:535003/251‑1 (MQ=255)                                                                                                                                                                                                               gCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTcgc                      >  2:425317/1‑251 (MQ=255)                                                                                                                                                                                                                 cAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCgg                    >  1:265953/1‑251 (MQ=255)                                                                                                                                                                                                                  aaTTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGt                   >  2:259771/1‑251 (MQ=255)                                                                                                                                                                                                                              aaagtGCGGCTGTTTTGCCCCAGATTTATCGGCTCAGCAACAC‑GGACGCTCATCTTTCCACTGCGCGGAGAACGCCGGGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGTTATGttcttc         <  2:84048/247‑1 (MQ=255)                                                                                                                                                                                                                               aa‑TGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGATCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGTTATGTTCTTcc        >  2:428685/1‑250 (MQ=255)                                                                                                                                                                                                                                   gCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGGGCAATTCGCGGTTATGTTCTTCCAGtt    >  2:192021/1‑251 (MQ=255)                                                                                                                                                                                                                                     ggCTGATTTGCTCCAGCATTCTCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCCAAGAACGCAGAGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGATCGCGTTCCCTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGTTTATGTTCTTCCAGTTcg  >  2:64814/1‑251 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       AGCACGTTGCCCATCACCTCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCA‑GCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGC‑AGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGG‑TAGTAGCCAATTTG‑CTGTGAAA‑TGCGGCTGATTTGCTCCAGCATTACCGGCTCAGCATCAC‑TGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGTTATGTTCTTCCAGTTCG  >  NC_000913/1188617‑1189090

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