Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	310,256	A→C	6.2%	L139R (CTG→CGG)	ecpA ←	ECP pilin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	310,256	0	A	C	6.2%	279.6 / 10.8	81	L139R (CTG→CGG)	ecpA	ECP pilin
Reads supporting (aligned to +/- strand):  ref base A (34/42);  new base C (2/3);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-01

CTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAA  >  NC_000913/310108‑310403                                                                                                                                                     |                                                                                                                                                    cTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCa                                                                                                                                                     <  3:294779/149‑1 (MQ=255)    tACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCg                                                                                                                                                  <  3:489997/149‑1 (MQ=255)          aaCGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTaa                                                                                                                                            >  3:256229/1‑149 (MQ=255)             gTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCCGCTGGCTaacgg                                                                                                                                         <  1:168290/149‑4 (MQ=255)              tCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGtt                                                                                                                                        <  4:389555/149‑1 (MQ=255)                 ccgcTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTgcc                                                                                                                                     <  1:118721/149‑1 (MQ=255)                              cTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCCGCTGGCTaacggtta                                                                                                                                      >  3:562331/1‑129 (MQ=255)                              cTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCCGCTGGCTaacggtta                                                                                                                                      <  4:562331/135‑7 (MQ=255)                               ttCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGt                                                                                                                       >  4:537543/1‑149 (MQ=255)                                 ccGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTg                                                                                                                     <  2:390499/149‑1 (MQ=255)                                    gTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCg                                                                                                                  >  1:258104/1‑149 (MQ=255)                                       gagTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTa                                                                                                               <  3:313047/149‑1 (MQ=255)                                            cTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGat                                                                                                          <  4:73420/149‑1 (MQ=255)                                             tGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGc                                                                                                                   >  1:30389/1‑139 (MQ=255)                                             tGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGc                                                                                                                   <  2:30389/139‑1 (MQ=255)                                                aaTCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGatcatc                                                                                                      >  3:262500/1‑149 (MQ=255)                                                    tGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACgg                                                                                                  >  4:409150/1‑149 (MQ=255)                                                       tACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTAt                                                                                               >  3:201218/1‑149 (MQ=255)                                                          tGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGc                                                                                            <  3:337199/149‑1 (MQ=255)                                                          tGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGc                                                                                            >  4:319871/1‑149 (MQ=255)                                                          tGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGc                                                                                            >  1:446097/1‑149 (MQ=255)                                                               tggtACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGtt                                                                                       <  2:316198/149‑1 (MQ=255)                                                                  tACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGttttt                                                                                    >  3:237975/1‑149 (MQ=255)                                                                    ccATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTc                                                                                  <  2:245012/149‑1 (MQ=255)                                                                    ccATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTc                                                                                  >  4:181772/1‑149 (MQ=255)                                                                    ccATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTc                                                                                  <  2:333272/149‑1 (MQ=255)                                                                      aTTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGa                                                                                <  3:444317/149‑1 (MQ=255)                                                                       ttGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGAc                                                                               >  3:126389/1‑149 (MQ=255)                                                                         ggtggtACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTg                                                                             <  2:587755/149‑1 (MQ=255)                                                                             gtACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCcgc                                                                         >  2:465039/1‑149 (MQ=255)                                                                                      atgatgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCGGCGCCGTtata                                                                <  3:342107/149‑1 (MQ=255)                                                                                       tgatgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATaa                                                               <  2:258104/149‑1 (MQ=255)                                                                                        gatgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAAt                                                              >  2:110508/1‑149 (MQ=255)                                                                                        gatgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAAt                                                              <  4:201218/149‑1 (MQ=255)                                                                                         atgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATc                                                             <  2:579198/149‑1 (MQ=255)                                                                                         atgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATc                                                             >  2:278009/1‑149 (MQ=255)                                                                                              gaAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGc                                                        >  4:458769/1‑149 (MQ=255)                                                                                              gaAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGc                                                        <  2:14362/149‑1 (MQ=255)                                                                                              gaAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGc                                                        <  1:213884/149‑1 (MQ=255)                                                                                              gaAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTACGACTGCCGCGCCGTTATAATCCACGc                                                        >  4:378338/1‑149 (MQ=255)                                                                                                    gAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACAt                                                  <  1:278009/149‑1 (MQ=255)                                                                                                       aCCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTTCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCa                                               >  4:261335/1‑149 (MQ=255)                                                                                                       aCCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCa                                               <  4:86611/149‑1 (MQ=255)                                                                                                       aCCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCa                                               <  1:393657/149‑1 (MQ=255)                                                                                                       aCCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCa                                               <  1:465039/149‑1 (MQ=255)                                                                                                        ccATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAg                                              >  1:458701/1‑149 (MQ=255)                                                                                                         cATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAgt                                             <  4:256229/149‑1 (MQ=255)                                                                                                         cATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAgt                                             >  3:95276/1‑149 (MQ=255)                                                                                                                    gtggtACGATTGCTGGCATTGTAAGCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACTCCCACATTCAGTGTGGAACCTGa                                  <  2:157999/149‑1 (MQ=255)                                                                                                                    gtggtACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGa                                  <  3:255263/149‑1 (MQ=255)                                                                                                                    gtggtACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGa                                  <  4:262500/149‑1 (MQ=255)                                                                                                                     tggtACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAg                                 <  3:458769/149‑1 (MQ=255)                                                                                                                     tggtACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAg                                 <  4:32629/149‑1 (MQ=255)                                                                                                                     tggtACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAg                                 <  2:575317/149‑1 (MQ=255)                                                                                                                      ggtACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAgg                                <  4:89023/149‑1 (MQ=255)                                                                                                                       gtACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGt                               >  1:546745/1‑149 (MQ=255)                                                                                                                           gATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATcc                           >  1:498015/1‑149 (MQ=255)                                                                                                                           gATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATcc                           >  2:421613/1‑149 (MQ=255)                                                                                                                             ttGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCaa                         >  4:165526/1‑149 (MQ=255)                                                                                                                              tGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAAc                        <  3:261335/149‑1 (MQ=255)                                                                                                                                 tGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCca                                                     <  1:315903/117‑1 (MQ=255)                                                                                                                                 tGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCca                                                     >  2:315903/1‑117 (MQ=255)                                                                                                                                 tGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTgg                     <  3:379360/149‑1 (MQ=255)                                                                                                                                 tGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTgg                     >  4:484386/1‑149 (MQ=255)                                                                                                                                 tGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTgg                     <  4:237975/149‑1 (MQ=255)                                                                                                                                 tGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTgg                     <  2:458701/149‑1 (MQ=255)                                                                                                                                   gCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGt                   >  2:14391/1‑149 (MQ=255)                                                                                                                                     aTTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTa                 >  2:303625/1‑149 (MQ=255)                                                                                                                                      ttGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTaa                >  1:451435/1‑149 (MQ=255)                                                                                                                                       tGTAACCGTTAGCCATCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTCATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAAt               <  2:173303/149‑1 (MQ=255)                                                                                                                                       tGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAAt               <  1:143399/149‑1 (MQ=255)                                                                                                                                         tAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAAtgt             >  4:285482/1‑149 (MQ=255)                                                                                                                                           aCCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATgtgt           >  2:332906/1‑149 (MQ=255)                                                                                                                                             ccttAGCCCGCTGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGtt          >  4:547377/3‑148 (MQ=255)                                                                                                                                             ccttAGCCCGCTGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGtt          <  3:547377/146‑1 (MQ=255)                                                                                                                                             cGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTg         >  4:530742/1‑149 (MQ=255)                                                                                                                                             cGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACCGGGTTAATGTGTTg         <  2:280140/149‑1 (MQ=255)                                                                                                                                              gTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAg                                 <  1:588273/124‑1 (MQ=255)                                                                                                                                              gTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAg                                 >  2:588273/1‑124 (MQ=255)                                                                                                                                              gTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTgg        <  3:530742/149‑1 (MQ=255)                                                                                                                                               ttAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGt       <  3:30275/149‑1 (MQ=255)                                                                                                                                                    cAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCACCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATaa  >  2:487269/1‑149 (MQ=255)                                                                                                                                                     aGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGAt    <  4:529211/146‑1 (MQ=255)                                                                                                                                                     aGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGAt    >  3:529211/1‑146 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAA  >  NC_000913/310108‑310403

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