Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,075,528	C→T	A1119T (GCA→ACA)	putA ←	fused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,075,528	0	C	T	94.7%	103.2 / ‑6.4	41	A1119T (GCA→ACA)	putA	fused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (21/15);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGCGCGTATGGACGCCAAGCGTCAGACCATAACCGGAAGCGTTAATCTGCTCGATCA  >  NC_000913/1075282‑1075773                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      gCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAgc                                                                                                                                                                                                                                                   >  2:130429/1‑251 (MQ=255) gCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAg                                                                                                                                                                                                                                                    <  1:56962/250‑1 (MQ=255)          cGGCAATAGACCACACGCGCACACGCAGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGACAGCAATCGTTGTGTTCCTGCCTGCGCCATCTCGCCATCTTGCGTACATAACGCCGGCAATTCTGGACGATTTGCTGCCCATTCCCGCACTGCATTTCGCGGCTGAGTCAATGCGGCTTTCAACTACGCATCGCCCAGATACTTTGCATACTG‑ACGCGCGCGCGGCACTCACAACGCACTTTcc                                                                                                                                                                                                                                          >  2:56045/1‑251 (MQ=255)          cGGCAATACACAACACGAGCTCACGCGGCAGCAGCGTACAGGTGTTGCGTTCAACCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGACCTTTcc                                                                                                                                                                                                                                          >  1:399217/1‑251 (MQ=255)            gCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCgg                                                                                                                                                                                                                                        <  1:437409/251‑1 (MQ=255)                   acaacaCGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTc                                                                                                                                                                                                                                 <  1:276861/251‑1 (MQ=255)                    caacaCGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCg                                                                                                                                                                                                                                <  1:55030/251‑1 (MQ=255)                      acacGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCCGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑CCGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTcccc                                                                                                                                                                                                                              >  2:436489/1‑248 (MQ=255)                       cacGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGccc                                                                                                                                                                                                                             >  1:480016/1‑250 (MQ=255)                          gcgcTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCcagca                                                                                                                                                                                                                          <  1:489127/251‑1 (MQ=255)                                 ggctgaagcatctcccggttgtTGCTTTCACCCGTCCGTCCGGGCCGCAATCGTTGTGGTACTGCCTGAGCCAGCTCGCCATATTGCGTACATATCGCCTGCAATTCTGGACGATTTGCTGCACATTCCCTCAGTGCATTTAGCGGCTGAGTCAATGCGTCTTTCAACTGCGCCTCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGTAGACGGTa                                                                                                                                                                                                                   <  2:93910/232‑1 (MQ=255)                                        gcagcGTCCAGGTTGTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCTCCATATTGCGTACATAACGCCTGCAATTCTTGACGATTTGCTGCCCATTCCCGCAGTGCATTGAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATTCTTTGCATCCTT‑ACGCGCGAGCGTCACTGTCAGCTCACTTTCCTGGCGATTCGCCAGCTGACGGTAGAGATag                                                                                                                                                                                                            <  2:454020/251‑1 (MQ=255)                                             gTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCCCTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCgg                                                                                                                                                                                                       >  2:466705/1‑251 (MQ=255)                                                           ttCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCg                                                                                                                                                                                         >  1:43110/1‑251 (MQ=255)                                                                 cGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCgg                                                                                                                                                                                   >  1:59693/1‑251 (MQ=255)                                                                    cGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTAc                                                                                                                                                                                >  1:272186/1‑251 (MQ=255)                                                                                              ctgcctgcGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGCACCGGACAACCCTTCGCCGCCGAACGGc                                                                                                                                                      >  1:209953/1‑251 (MQ=255)                                                                                                                   atTGCGTACATAACACCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGAGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGccc                                                                                                                                 >  1:380762/1‑251 (MQ=255)                                                                                                                        gTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGTATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTGCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCat                                                                                                                            <  2:172719/251‑1 (MQ=255)                                                                                                                            aTAACGCCTGCAATACTGTACGATTTGCTGCCCATTCCCGCCGTGCATTTAGCGGCAGAGTCACTGCGGCTTTCACCTGCGCATCGTCCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCCCTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGTCCGCCTGCTTTCGGCCCGGTTCCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCTCTGCGCCCACCATATTa                                                                                                                        <  2:165711/251‑1 (MQ=255)                                                                                                                              aaCGCCTGCAATTCTGGACGAGTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACg                                                                                                                      <  1:163265/251‑1 (MQ=255)                                                                                                                                cgtcTGCAATTATGGACGATTTGTTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGTCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACGGTCAGCGCTCTTTCCGGTCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCGGCACACCAACCACTGCGCCCACCATATTCCGGt                                                                                                                    <  2:292265/248‑1 (MQ=255)                                                                                                                                   cTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGCCCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGCTTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCCCACCAACCCCTGCGCCCACCCTATTACGGTTa                                                                                                                  >  2:272097/1‑250 (MQ=255)                                                                                                                                    tGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTCGCGGCTGAGTCAATGCGGCTTTCACCTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGCCGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCCCACCCACCACTGCGCCCACCATATTACGGTTAAc                                                                                                                >  2:466590/1‑251 (MQ=255)                                                                                                                                               acattttgttgcccgatccccgatcgCATTTAGCGGCTGACTCAATGCGGCTTTCAACTGTGCATCGACCGGATATTTTGC‑TACTGGACGCTGGAGCGTCACGTCAAGCGCACTTTCCGGGCGATGCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCCGCCCGGTACCGCACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCAACATATTACGGTTAACATACAGGTTAc                                                                                                     <  2:490406/226‑1 (MQ=255)                                                                                                                                                       ctgcCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATgg                                                                                             >  2:229319/1‑251 (MQ=255)                                                                                                                                                        ggccaaTTCCCGCAGTGCATTTAGCGGCTGAGTCAAGGCGGCTTTCAACCGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGGCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACGGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATggg                                                                                            <  2:390099/246‑1 (MQ=255)                                                                                                                                                          cccATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGcc                                                                                          >  2:386905/1‑251 (MQ=255)                                                                                                                                                                    cAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCCCTGTCAGAGCACTTTCCGGGCGTTTCGCCAGCCGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGCACGGCTGCACACCAACCACTGCGCCCACCATACTACGGTTAACATACAGGTTACCAACACAGGCCGAGCCAGtc                                                                                 >  2:443835/1‑249 (MQ=255)                                                                                                                                                                       tGCATTTAGCGGCTGCGTCAATGAGGCTTTCAACTGCGCATCGACAGGTTACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCACAGCCCTTTCCGGGCGATTCGCCAGCAGCCGGTAGAGATAGAGCGGACAGCCTGCTTTCGGCCCGGTACCGGCCAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCCGATTACGTTTAACATCCCGGTTAACAACCTGGGCAGAGCCAGTGACCt                                                                             >  2:199892/1‑251 (MQ=255)                                                                                                                                                                         cATTTAGCGGCTGAGCAAATGCGGCTTTCAACTGCGCACCGACCGGATTCTTTGCAACCTG‑ACGCGCGGGCGACATAGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGTGAAAGCGGACCCCCTGCTTCCGGCACGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCCCCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTgg                                                                           <  2:35905/251‑1 (MQ=255)                                                                                                                                                                                      aGTCAATGCGTCTTCCAACTGCGCATCGACCGGAGACTTTGCATCCTG‑TCGCGCGAGCGTCCCTGTCAGCGCACTTTCCGGGCGAGTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTCCCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTtcat                                                              <  1:199892/251‑1 (MQ=255)                                                                                                                                                                                               ggCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGcgcg                                                     <  1:471863/251‑1 (MQ=255)                                                                                                                                                                                                          gcgcAGCGACGGGATCCTTTGCATCCTG‑ACGCGCGAGGGGCACTGACAGCGCCCTTTCCGGGCGATTCGCCAGCATACGGTAGAGATAGAGCGGCCCGCCTGCTTTCGGCCCGGTACCGGACAACCCTGCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGCGCGTATGGACGCCa                                          <  2:152716/251‑1 (MQ=255)                                                                                                                                                                                                           cgcATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATTCGCGTATGGACGCCa                                          >  2:89302/1‑250 (MQ=255)                                                                                                                                                                                                                    cGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGCGCGTATGGACGCCAAGCGTCAGAc                                >  2:118504/1‑251 (MQ=255)                                                                                                                                                                                                                       gattttttgcctccacgcgcgcgAGCGGCACTGTCAGCGCACTTTCCGGTCGTTGCGCCAGAAGACGGTAGAGAGAGAGCGGACCGCCTGCGTCCGGCCGGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCCCCATATTACGGTTAACATACAGGTTTCCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGCGCGTATGGACGCCAAGCGTCATACCa                              <  2:274713/234‑1 (MQ=255)                                                                                                                                                                                                                              gCATCCTG‑ACGCGCGAGCGTCACTGTCAGCGCACTTGCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCTCTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGCGCGTATGGACGCCAAGCGTCAGACCATAACCgg                       <  1:466590/250‑1 (MQ=255)                                                                                                                                                                                                                                        cgcgcgAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGCGCGCATGGACGCCACGCGTCAGACCATAACCGGACGCGTTAAt             >  2:218595/1‑251 (MQ=255)                                                                                                                                                                                                                                            cgAGCGTCACTGTCAGCGCACTTTCCGGGCGATTCGCCAGAAGACGGAAGCGCTAGAGCGGACCGCCTGCTTTCGGCCTGGTACCGGACAAACCTTCGCCGCCGAACGGCTGCACACCAACCACTGACCCCACCATATTACGGTTCACATACAGGTTACCAACATGGGCCGAGCCATTGACCTGGGCGATGGTTTCATCAATGCGCGTATGGACGCCAAGCGTCAGACCATAACCGGAAGCGTTACTCTGc         >  1:239835/1‑251 (MQ=255)                                                                                                                                                                                                                                                   cACTGTCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGCGCCAGTGACCTGGGCGCTGGTTTAATCAATGCGCGTATGGACGCCAAGCGTCAGACCATAACCGGAAGCGTTAATCGGCTCGATcc  >  2:256033/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      GCTCATCATCGGCAATACACAACACGCGCTCACGCGGCAGCAGCGTCCAGGTGTTGCGTTCACCCGTCGGCCCCGGCAGCAATCGTTGTGTTCCTGCCTGCGCCAGCTCGCCATATTGCGTACATAACGCCTGCAATTCTGGACGATTTGCTGCCCATTCCCGCAGTGCATTTAGCGGCTGAGTCAATGCGGCTTTCAACTGCGCATCGACCGGATACTTTGCATCCTG‑ACGCGCGAGCGTCACTGCCAGCGCACTTTCCGGGCGATTCGCCAGCAGACGGTAGAGATAGAGCGGACCGCCTGCTTTCGGCCCGGTACCGGACAACCCTTCGCCGCCGAACGGCTGCACACCAACCACTGCGCCCACCATATTACGGTTAACATACAGGTTACCAACATGGGCCGAGCCAGTGACCTGGGCGATGGTTTCATCAATGCGCGTATGGACGCCAAGCGTCAGACCATAACCGGAAGCGTTAATCTGCTCGATCA  >  NC_000913/1075282‑1075773

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