Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,132,960	C→A	6.6%	A14A (GCC→GCA)	tolC →	agglutination protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,132,960	0	C	A	6.6%	225.3 / 3.8	73	A14A (GCC→GCA)	tolC	agglutination protein
Reads supporting (aligned to +/- strand):  ref base C (30/37);  new base A (3/2);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.51e-01

ATGTACTTAGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTCCG  >  NC_002947/5132813‑5133108                                                                                                                                                    |                                                                                                                                                      aTGTACTTAGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTTCTGTTGGccc                                                                                                                                                      >  8:373976/1‑149 (MQ=255) aTGTACTTAGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGccc                                                                                                                                                      <  5:217596/149‑1 (MQ=255)  tGTACTTAGCTGTGTTGTTAAAAAGCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCt                                                                                                                                                     <  2:335463/149‑1 (MQ=255)      cTTAGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGc                                                                                                                                                  <  4:359343/148‑1 (MQ=255)         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                                                                                                                           <  6:339744/147‑1 (MQ=255)                              cAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGt                                                                                                                           <  8:103980/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                                                                                                <  6:224996/148‑1 (MQ=255)                                                       gTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGcc                                                                                                <  7:135164/148‑1 (MQ=255)                                                        tCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGt                                                                                              >  1:79055/1‑149 (MQ=255)                                                          aTTTTTTTGGCTAACCGCAGGATCAATTCCACGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTgg                                                                                            <  4:327122/149‑1 (MQ=255)                                                          aTTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTg                                                                                             >  8:262885/1‑148 (MQ=255)                                                          aTTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTg                                                                                             <  2:79055/148‑1 (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                                                                              <  7:282016/149‑1 (MQ=255)                                                                        cagcagGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACaa                                                                              <  4:37400/149‑1 (MQ=255)                                                                        cagcagGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACaa                                                                              <  1:29744/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                                                            <  7:373976/149‑1 (MQ=255)                                                                                          gATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAAc                                                            <  8:103068/149‑1 (MQ=255)                                                                                              aGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTaa                                                             <  4:201688/144‑1 (MQ=255)                                                                                              aGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTaa                                                             >  3:201688/1‑144 (MQ=255)                                                                                               gggAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCaa                                                       <  3:70944/149‑1 (MQ=255)                                                                                                     tGTACAATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GAACCGTGGACTACCACCAACAAGTCAGCTACAACCGTAACAGCAAGCTGTc                                                 >  1:280093/1‑149 (MQ=255)                                                                                                       tACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCaa                                                                                                       >  7:163439/1‑94 (MQ=255)                                                                                                       tACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCaa                                                                                                       <  8:163439/94‑1 (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                                                                                              >  4:253234/1‑94 (MQ=255)                                                                                                               gcgTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGt                                                                                              <  3:253234/94‑1 (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                                                       >  5:267756/1‑102 (MQ=255)                                                                                                                                              gTTGGCACTCGCCAGGGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCaa                                                       <  6:267756/102‑1 (MQ=255)                                                                                                                                              gTTGGCACTCGCCAGGGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCAGCCGATGAAGACGTGAAGTTCGCGCgtg                    <  2:346362/137‑1 (MQ=255)                                                                                                                                              gTTGGCACTCGCCAGGGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCAGCCGATGAAGACGTGAAGTTCGCGCgtg                    >  1:346362/1‑137 (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)                                                                                                                                                    |                                                                                                                                                      ATGTACTTAGCTGTGTTGTTAAAAACCTGCAACACTATACGCGATAAGTTTATCGTCATTTTTTTGGCGAACAGCAGGATCAATTCCAGGATCAGGGAGATGTACCATGCGCGTTTTGAACCCCATCACCAGTGCTTTACTGTTGGCCCTGGCGAGTGCCAACGTACAGGCGATGTCGATTACCGAGGCTGTCCAAA‑GCGCCGTGGACTACCACCCACAAGTCAGCTCCAACCGTAACAGCAAGCTGTCGGCCGATGAAGACGTGAAGTTCGCGCGTGGTGGCTACTACCCTTCCG  >  NC_002947/5132813‑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: ‑