Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	310,310	T→C	6.1%	D121G (GAT→GGT)	ecpA ←	ECP pilin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	310,310	0	T	C	6.1%	260.4 / 6.9	83	D121G (GAT→GGT)	ecpA	ECP pilin
Reads supporting (aligned to +/- strand):  ref base T (34/43);  new base C (3/2);  total (37/46)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

GGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAGTTTAAAGGCGGTAGCCGTTGAGTCACCCTCG  >  NC_000913/310171‑310447                                                                                                                                            |                                                                                                                                          ggtACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTAt                                                                                                                                           >  2:342394/1‑140 (MQ=255)       aatggtcgTACAGCAGAAGATGGAGAAGGTGAAACCATCCTGTGGGGTGGAACGATTGATTGCATAGTAACCCTAAGACAGAGCGCTAAGGTTGGCGCCCAGTACGACGTTGTCGGTATCGATCATCACGGTATCCCCAg                                                                                                                                     >  2:953456/3‑140 (MQ=255)       aTTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTACCGCCAg                                                                                                                                     >  2:2854/1‑140 (MQ=255)       aTTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTACCGCCAg                                                                                                                                     >  2:3134/1‑140 (MQ=255)       aTTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTACCGCCAg                                                                                                                                     >  2:3422/1‑140 (MQ=255)              gtACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTc                                                                                                                              <  1:1611908/140‑1 (MQ=255)                aCCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGa                                                                                                                            >  2:995761/1‑140 (MQ=255)                aCCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGa                                                                                                                            >  2:1239199/1‑140 (MQ=255)                aCCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGa                                                                                                                            >  2:1463333/1‑140 (MQ=255)                aCCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGa                                                                                                                            >  2:1586727/1‑140 (MQ=255)                aCCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGa                                                                                                                            >  2:72600/1‑140 (MQ=255)                aCCGCTGATGAGGGAGTAGGTGAAACCATCCTGTGCGGTGGTACGCTTGCAGGCCTTGTAACCGTTAGCCAGCGGGCTGAGGGAGCCGCCCCGTACGCCGGTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTGGa                                                                                                                            >  2:811400/1‑140 (MQ=255)                aCCGCAGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTAAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGa                                                                                                                            >  2:377773/1‑140 (MQ=255)                 ccGCTGATGATGGAGTAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTATGGTTGCCGCCCAGTACGCCGTGGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGAc                                                                                                                           >  2:834057/1‑140 (MQ=255)                 ccGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGAc                                                                                                                           >  2:394672/1‑140 (MQ=255)                 ccGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGGTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGAc                                                                                                                           >  2:110279/1‑140 (MQ=255)                    cTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGc                                                                                                                        >  1:1367517/1‑140 (MQ=255)                    cTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGc                                                                                                                        >  1:599355/1‑140 (MQ=255)                      gatgatGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCcg                                                                                                                      >  2:457451/1‑140 (MQ=255)                      gatgatGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCcg                                                                                                                      >  2:457456/1‑140 (MQ=255)                      gatgatGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCcg                                                                                                                      <  1:1107994/140‑1 (MQ=255)                        tgatgGAGATGGTGAAACCATCCTGTGCGGTGGAACGATTGCTGGCATGTAAACCGTTAGCCAGCGGGCTAAGGTTGCCCCCCAGTACGCGGTTGGGGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCcgcg                                                                                                                    <  2:832215/140‑1 (MQ=255)                        tgatgGAGAAGGTGACACCATCCTGTGCGGTGGTACGATAGCTGGCATTGTAACCGTTAGCCAGCGGGATAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCcgcg                                                                                                                    <  2:832362/140‑1 (MQ=255)                        tgatgGAGAAGGTGAAACCATCCTGTGCGGTGGTTCGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCcgcg                                                                                                                    <  2:285557/140‑1 (MQ=255)                        tgatgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCcgcg                                                                                                                    <  1:4203/140‑1 (MQ=255)                        tgatgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCcgcg                                                                                                                    <  1:50305/140‑1 (MQ=255)                        tgatgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCcgcg                                                                                                                    <  1:198671/140‑1 (MQ=255)                        tgatgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCcgcg                                                                                                                    <  2:290366/140‑1 (MQ=255)                        tgatgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGAGTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGTCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCcgcg                                                                                                                    <  1:1594320/140‑1 (MQ=255)                         gatgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTATGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCgcgc                                                                                                                   <  1:1547119/140‑1 (MQ=255)                            ggAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGt                                                                                                                >  1:894887/1‑140 (MQ=255)                            ggAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGt                                                                                                                >  1:1425783/1‑140 (MQ=255)                            ggAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGt                                                                                                                >  1:1377028/1‑140 (MQ=255)                              agaAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTta                                                                                                              >  2:1498861/1‑140 (MQ=255)                              agaAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTta                                                                                                              <  2:1088133/140‑1 (MQ=255)                              agaAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTta                                                                                                              <  2:753578/140‑1 (MQ=255)                                     gAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCa                                                                                                       <  2:1314897/140‑1 (MQ=255)                                       aaCCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACg                                                                                                     <  1:555475/140‑1 (MQ=255)                                            tCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCcac                                                                                                >  1:1538365/1‑140 (MQ=255)                                             ccTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTTTATTCCACGCCcaca                                                                                               <  2:422603/140‑1 (MQ=255)                                               tgtgCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACAtt                                                                                             <  1:134320/140‑1 (MQ=255)                                               tgtgCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACAtt                                                                                             <  1:74672/140‑1 (MQ=255)                                               tgtgCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACAtt                                                                                             <  2:1418203/140‑1 (MQ=255)                                                gtgCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTc                                                                                            <  1:672056/140‑1 (MQ=255)                                                 tgCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCa                                                                                           >  2:210312/1‑140 (MQ=255)                                                 tgCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCAACACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCa                                                                                           <  2:1612397/140‑1 (MQ=255)                                                     gtggtACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAgtgt                                                                                       <  1:774256/140‑1 (MQ=255)                                                     gtggtACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAgtgt                                                                                       <  1:1144686/140‑1 (MQ=255)                                                     gtggtACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAgtgt                                                                                       <  1:1076324/140‑1 (MQ=255)                                                      tggaACGATAGCTGGCATTGTAACCGTAAGACAGCGGGCTAAGGTTGCCGCCCAGTACGCCCATGGCGGTATCGATCCTCACGGTATCGCCAGGTGTTTCGACTGCCGCGCCGGTATAATCCACGCCCACAGACAGagtg                                                                                      >  2:657752/1‑140 (MQ=255)                                                       ggtACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTgg                                                                                     <  2:1060193/140‑1 (MQ=255)                                                           cGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAAcc                                                                                 >  1:1581989/1‑140 (MQ=255)                                                                       ttGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCaa                                                                     <  2:648545/140‑1 (MQ=255)                                                                       ttGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCaa                                                                     <  1:1098039/140‑1 (MQ=255)                                                                       ttGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCaa                                                                     <  1:1149535/140‑1 (MQ=255)                                                                       ttGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCaa                                                                     <  1:446257/140‑1 (MQ=255)                                                                       tagTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTATTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCaa                                                                     <  2:80369/138‑1 (MQ=255)                                                                         gTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACt                                                                   <  1:354542/140‑1 (MQ=255)                                                                                  aGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAAtg                                                          >  1:1619505/1‑140 (MQ=255)                                                                                  aGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGATTTTTCGACTGCCGCGCCGATATAATCCACCCCCACGTTCAGTGTGGAACCTGAGGTATCCAACTGGGATaaag                                                          >  1:834226/1‑138 (MQ=255)                                                                                       gCGTGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTg                                                     >  1:1496553/1‑140 (MQ=255)                                                                                        cGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAGCTGGGTTAATGTGTTgg                                                    >  2:1121130/1‑140 (MQ=255)                                                                                            cTAAGGTTGCCGCCCAGTTAGCCGTTGGCGGTATCGATCATCACGTTATCGCCAGATTTTTCGACTGCCGCGCCGTTATAATCCACGCCCCCATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGAt                                                <  2:615087/140‑1 (MQ=255)                                                                                              aaGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATaa                                              <  1:1463333/140‑1 (MQ=255)                                                                                              aaGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATaa                                              <  1:1586727/140‑1 (MQ=255)                                                                                                     ccgccCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGTTTATAAGACGTGa                                       >  1:817180/1‑140 (MQ=255)                                                                                                     ccgccCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGa                                       >  1:660090/1‑140 (MQ=255)                                                                                                        cccAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGt                                    <  1:1137666/140‑1 (MQ=255)                                                                                                        cccAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGt                                    <  1:961660/140‑1 (MQ=255)                                                                                                        cccAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGt                                    <  1:959362/140‑1 (MQ=255)                                                                                                        cccAGTACGCCGTTGGCGGTATCGATCATCACGGAATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGt                                    <  2:648407/140‑1 (MQ=255)                                                                                                        cacAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATACTCCTCGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGt                                    <  1:377773/138‑1 (MQ=255)                                                                                                           aGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAg                                 >  2:1253990/1‑140 (MQ=255)                                                                                                               cGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAGTTTa                             <  1:602835/140‑1 (MQ=255)                                                                                                                       gCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAGTTTAAAGGCGGt                     <  1:1239199/140‑1 (MQ=255)                                                                                                                       gCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAGTTTAAAGGCGGt                     >  1:1680782/1‑140 (MQ=255)                                                                                                                            aTCGATCATCTCGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAGTTTAAAGGCGGTAGCCg                <  1:420097/140‑1 (MQ=255)                                                                                                                            aTCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAGTTTAAAGGCGGTAGCCg                <  1:232379/140‑1 (MQ=255)                                                                                                                                  catcaCGGTAGCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAGTTTAAAGGCGGTAGCCGTTGAGt          <  1:3422/140‑1 (MQ=255)                                                                                                                                  catcaCGGTACCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAGTTTAAAGGCGGTAGCCGTTGAGt          <  1:3134/140‑1 (MQ=255)                                                                                                                                  catcaCGGTACCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAGTTTAAAGGCGGTAGCCGTTGAGt          <  1:2854/140‑1 (MQ=255)                                                                                                                                     caCGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAGTTTAAAGGCGGTAGCCGTTGAGTCAc       >  2:1102552/1‑140 (MQ=255)                                                                                                                                       cGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAGTTTAAAGGCGGTAGCCGTTGAGTCAccc     >  1:162498/1‑140 (MQ=255)                                                                                                                                          tATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAGTTTAAAGGCGGTAGCCGTTGAGTCACCCTCg  >  2:310324/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                          GGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAAGACGTGAGGTCAGTTTAAAGGCGGTAGCCGTTGAGTCACCCTCG  >  NC_000913/310171‑310447

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