Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	310,259	G→T	6.1%	P138Q (CCG→CAG)	ecpA ←	ECP pilin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	310,259	0	G	T	6.1%	288.3 / 7.7	82	P138Q (CCG→CAG)	ecpA	ECP pilin
Reads supporting (aligned to +/- strand):  ref base G (35/42);  new base T (2/3);  total (37/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 = 6.80e-01

TACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAA  >  NC_000913/310111‑310403                                                                                                                                                     |                                                                                                                                                 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                                                                                            >  4:319871/1‑149 (MQ=255)                                                       tGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGc                                                                                            <  3:337199/149‑1 (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                                                              <  4:201218/149‑1 (MQ=255)                                                                                     gatgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAAt                                                              >  2:110508/1‑149 (MQ=255)                                                                                      atgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATc                                                             <  2:579198/149‑1 (MQ=255)                                                                                      atgGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATc                                                             >  2:278009/1‑149 (MQ=255)                                                                                           gaAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGc                                                        <  2:14362/149‑1 (MQ=255)                                                                                           gaAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGc                                                        >  4:458769/1‑149 (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                                               <  1:465039/149‑1 (MQ=255)                                                                                                    aCCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCa                                               <  4:86611/149‑1 (MQ=255)                                                                                                    aCCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCa                                               <  1:393657/149‑1 (MQ=255)                                                                                                     ccATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAg                                              >  1:458701/1‑149 (MQ=255)                                                                                                      cATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAgt                                             >  3:95276/1‑149 (MQ=255)                                                                                                      cATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAgt                                             <  4:256229/149‑1 (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                                 <  2:575317/149‑1 (MQ=255)                                                                                                                  tggtACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAg                                 <  4:32629/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)                                                                                                                                                     |                                                                                                                                                 TACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGTGCGGTGGTACGATTGCTGGCATTGTAACCGTTAGCCAGCGGGCTAAGGTTGCCGCCCAGTACGCCGTTGGCGGTATCGATCATCACGGTATCGCCAGTTTTTTCGACTGCCGCGCCGTTATAATCCACGCCCACATTCAGTGTGGAACCTGAGGTATCCAACTGGGTTAATGTGTTGGTGATAA  >  NC_000913/310111‑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: ‑