Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,132,966	G→C	6.6%	A16A (GCG→GCC)	tolC →	agglutination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,132,966	0	G	C	6.6%	233.5 / 4.4	74	A16A (GCG→GCC)	tolC	agglutination protein
Reads supporting (aligned to +/- strand):  ref base G (32/37);  new base C (3/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

AGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTCCG  >  NC_002947/5132821‑5133108                                                                                                                                                  |                                                                                                                                                aGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGt                                                                                                                                              >  2:273245/1‑149 (MQ=255)  gCTGTGTTGTTAAAAAGCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAAGAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTg                                                                                                                                             <  2:237254/149‑1 (MQ=255)        ttgttAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACg                                                                                                                                       >  7:103980/1‑149 (MQ=255)           ttAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCTATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCGAACGTAc                                                                                                                                    >  6:78789/1‑149 (MQ=255)           ttAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTa                                                                                                                                     >  1:346746/1‑148 (MQ=255)            tAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTAATGTTGGCCCTGGCGAGTGCCAACGTa                                                                                                                                     >  3:327122/1‑147 (MQ=255)             aaaaaCCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAg                                                                                                                                  >  4:278485/1‑149 (MQ=255)             aaaaaCCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAAACCAATCACCACTGCTTTACTGTTGGCCCTCGCGAGTGCCAAAGAACac                                                                                                                                  >  1:86736/1‑148 (MQ=255)               aaaCCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGc                                                                                                                                >  8:282016/1‑149 (MQ=255)                      cAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGt                                                                                                                           <  8:103980/147‑1 (MQ=255)                      cAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGt                                                                                                                           <  6:339744/147‑1 (MQ=255)                      cAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTc                                                                                                                          <  1:25938/148‑1 (MQ=255)                           cTATACGCGATAAGTTTATCGTGATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGTGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTAc                                                                                                                    <  7:130268/149‑1 (MQ=255)                            tataCGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTAcc                                                                                                                   >  7:103068/1‑149 (MQ=255)                            tataCGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTAcc                                                                                                                   >  7:79043/1‑149 (MQ=255)                              taCGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCg                                                                                                                  <  1:273245/148‑1 (MQ=255)                                cgcgATAAGTTTATCGTCATTTTTATGTCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGGTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAgg                                                                                                               <  6:76939/149‑1 (MQ=255)                                     tAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTc                                                                                                          <  3:278485/149‑1 (MQ=255)                                         tttATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCaaa                                                                                                      <  4:308982/149‑1 (MQ=255)                                           tATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑gc                                                                                                   >  4:240826/1‑149 (MQ=255)                                               gTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGcc                                                                                                <  7:135164/148‑1 (MQ=255)                                               gTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGcc                                                                                                <  6:224996/148‑1 (MQ=255)                                                tCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGt                                                                                              >  1:79055/1‑149 (MQ=255)                                                  aTTTTTTTGGCTAACCGCAGGATCAATTCCACGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTgg                                                                                            <  4:327122/149‑1 (MQ=255)                                                  aTTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTg                                                                                             <  2:79055/148‑1 (MQ=255)                                                  aTTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTg                                                                                             >  8:262885/1‑148 (MQ=255)                                                       tttGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTa                                                                                        <  3:240826/148‑1 (MQ=255)                                                         tGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTacca                                                                                     >  4:70944/1‑149 (MQ=255)                                                             gAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCaa                                                                                                                                         >  7:27988/1‑94 (MQ=255)                                                             gAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCaa                                                                                                                                         <  8:27988/94‑1 (MQ=255)                                                             gAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTAccaccca                                                                                 >  3:285676/1‑149 (MQ=255)                                                                cagcagGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACaa                                                                              <  1:29744/149‑1 (MQ=255)                                                                cagcagGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACaa                                                                              <  4:37400/149‑1 (MQ=255)                                                                cagcagGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACaa                                                                              <  7:282016/149‑1 (MQ=255)                                                                      gATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGc                                                                        >  5:201734/1‑149 (MQ=255)                                                                                aGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAATGCGCCGTGGACTACCACCCACAAGTCAGCTCCAAc                                                                  <  1:59849/146‑1 (MQ=255)                                                                                aGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAAc                                                                  >  2:59849/1‑145 (MQ=255)                                                                                  gATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAAc                                                            <  8:103068/149‑1 (MQ=255)                                                                                  gATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAAc                                                            <  7:373976/149‑1 (MQ=255)                                                                                      aGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTaa                                                             >  3:201688/1‑144 (MQ=255)                                                                                      aGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTaa                                                             <  4:201688/144‑1 (MQ=255)                                                                                       gggAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCaa                                                       <  3:70944/149‑1 (MQ=255)                                                                                             tGTACAATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GAACCGTGGACTACCACCAACAAGTCAGCTACAACCGTAACAGCAAGCTGTc                                                 >  1:280093/1‑149 (MQ=255)                                                                                               tACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCaa                                                                                                       <  8:163439/94‑1 (MQ=255)                                                                                               tACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCaa                                                                                                       >  7:163439/1‑94 (MQ=255)                                                                                                aCCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGc                                              >  4:390553/1‑149 (MQ=255)                                                                                                    tGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGAt                                          <  1:138654/149‑1 (MQ=255)                                                                                                      cgcgTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATg                                         <  6:201734/148‑1 (MQ=255)                                                                                                       gcgTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGt                                                                                              <  3:253234/94‑1 (MQ=255)                                                                                                       gcgTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGt                                                                                              >  4:253234/1‑94 (MQ=255)                                                                                                        cgTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAg                                      <  4:342499/149‑1 (MQ=255)                                                                                                                aCCCCATCACCAGTTCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCTAGCTGTCGGCCGATGAAGACGTGaa                               <  2:86736/148‑1 (MQ=255)                                                                                                                  cccATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAg                                      <  1:397669/139‑1 (MQ=255)                                                                                                                  cccATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAg                                      >  2:397669/1‑139 (MQ=255)                                                                                                                  cccATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGtt                            >  4:295503/1‑149 (MQ=255)                                                                                                                   ccATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTc                           >  7:6523/1‑149 (MQ=255)                                                                                                                    cATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTcg                          <  7:262885/149‑1 (MQ=255)                                                                                                                        aCCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTcgcgc                       <  4:285676/148‑1 (MQ=255)                                                                                                                          cAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCgtg                    >  3:248012/1‑149 (MQ=255)                                                                                                                              tcTTTACTGTTGGCCCTTGCGAGTGCCAAAGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCACCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCgtggtg                 <  2:280093/147‑1 (MQ=255)                                                                                                                               cTTTAGTGTTGGCGCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCTCTAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGtggtgg                <  7:270061/148‑1 (MQ=255)                                                                                                                                 ttACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGcta             >  1:297923/1‑149 (MQ=255)                                                                                                                                   tacgTTGGCACTCGCCAGGGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGctac            >  8:156887/4‑148 (MQ=255)                                                                                                                                     tGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGc               <  3:90542/143‑1 (MQ=255)                                                                                                                                     tGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGc               >  4:90542/1‑143 (MQ=255)                                                                                                                                     tGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGATCCAACCGTAATAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCtactac         >  4:59760/1‑149 (MQ=255)                                                                                                                                      gTTGGCACTCGCCAGGGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCaa                                                       <  6:267756/102‑1 (MQ=255)                                                                                                                                      gTTGGCACTCGCCAGGGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCaa                                                       >  5:267756/1‑102 (MQ=255)                                                                                                                                      gTTGGCACTCGCCAGGGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCAGCCGATGAAGACGTGAAGTTCGCGCgtg                    >  1:346362/1‑137 (MQ=255)                                                                                                                                      gTTGGCACTCGCCAGGGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCAGCCGATGAAGACGTGAAGTTCGCGCgtg                    <  2:346362/137‑1 (MQ=255)                                                                                                                                          gCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGc                                                         <  7:153119/96‑1 (MQ=255)                                                                                                                                           cccTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTcc   <  7:74195/149‑1 (MQ=255)                                                                                                                                           cccTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTcc   >  6:154309/1‑149 (MQ=255)                                                                                                                                            ccTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTCCg  >  6:101542/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                AGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTCCG  >  NC_002947/5132821‑5133108

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