Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,575,620	G→A	5.9%	A667V (GCA→GTA)	yddB ←	putative TonB‑dependent outer membrane receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,575,620	0	G	A	5.9%	226.4 / 2.3	68	A667V (GCA→GTA)	yddB	putative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base A (2/2);  total (35/33)
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.00e+00

GCTGATGTATTGCGCATTTGTCTTACCAAGGATAATGCGAGCTTTACGCGCTTCTTGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGATATAGCTAAAGACAATTTGTGGGTTAATATCTACCT  >  NC_000913/1575475‑1575767                                                                                                                                                  |                                                                                                                                                    gCTGATGTATTGCGCATTTGTCTTACCAAGGATAATGCGAGCTTTACGCGCTTCTTGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGcc                                                                                                                                                   <  6:290305/148‑1 (MQ=255)      tGTATTGCGCATTTGTCTTACCAAGGATAATGAGAGCTTTACGCGCTTCTTGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTg                                                                                                                                               >  2:82802/1‑147 (MQ=255)        ttttGCGCATTTGTCTTTCCAAGTATAATGCGAGCTTTACGCGCTTCTTGCCAGGCGAGCGTATTCGCCCATACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGATATATTTAATGGGTTATTAAATTCTGCCACTGGa                                                                                                                                             <  4:257864/145‑1 (MQ=255)           tGCGCATTTGTCTTACCAAGGATAATGCGAGCTTTACGCTCATCTTGCAAGGCCAGAGTATTCGCCCACACTAACCCGCTCGCTAGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTACCACTGCAACGc                                                                                                                                         >  1:159980/1‑148 (MQ=255)           tGCGCATTTGTCTTACCAAGGATAATGCGAGCTTTACGCGCTTCTTGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTACCACTGGAACGc                                                                                                                                         >  7:242090/1‑148 (MQ=255)             cgcATTTGTCTTACCAAGGATAATGCGAGCTTTACGCGCTTCTTGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTa                                                                                                                                       >  3:269340/1‑148 (MQ=255)                  ttGTCTTACCAAGGATAATGCGAGCTTTACGCGCTTCTTGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGt                                                                                                                                   >  3:84219/1‑147 (MQ=255)                  ttGTCTTACCAAGGATAATGCGAGCTTTACGCGCTTCTTGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGt                                                                                                                                   >  5:82083/1‑147 (MQ=255)                   tGTCTTACCAAGGATAATGCGAGCTTTACGCGCTTCTTGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGt                                                                                                                                   >  1:196129/1‑146 (MQ=255)                                         gCTTTACGCGCTTCTTGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATaaa                                                                                                           <  4:269340/148‑1 (MQ=255)                                               cgcgCTTCTTGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGTGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCt                                                                                                     <  8:77517/148‑1 (MQ=255)                                                   cttcttGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGtt                                                                                                 <  2:196129/148‑1 (MQ=255)                                                      gttGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTACCACTCGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTa                                                                                                <  2:159980/145‑1 (MQ=255)                                                      cttGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTACCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTa                                                                                                <  8:242090/146‑1 (MQ=255)                                                         gccaggaaagCGTATTCTACAACACTAACAAGCTCTGTTGATGCGTGAAATGCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTTGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATc                                                                                              <  1:262253/145‑1 (MQ=255)                                                            aggccagCGTATTCGCCCACACTACCCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGt                                                                                        <  1:82802/148‑1 (MQ=255)                                                            aggccagCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGt                                                                                        <  4:84219/148‑1 (MQ=255)                                                                  gCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGct                                                                                  >  2:9300/1‑148 (MQ=255)                                                                   cGTATTCGCCCACACTAACCCGCTCGGTTGATGCGGGAAATCCATGTTTAAGTAGATCTTTAATGGGGTATTAAAATCTGCCACTTGTACGCTATCGTTAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGct                                                                                  <  6:230582/147‑1 (MQ=255)                                                                    gTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGctct                                                                                >  3:299375/1‑148 (MQ=255)                                                                      aTTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTctcc                                                                              >  6:290572/1‑148 (MQ=255)                                                                       ttCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTctcct                                                                             >  3:120973/1‑148 (MQ=255)                                                                       ttCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTctcc                                                                              <  4:299375/147‑1 (MQ=255)                                                                            ccACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATaa                                                                        >  4:315988/1‑148 (MQ=255)                                                                              acacTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACATGGTTATCACCGGTATTGCTCTCCTCATAAc                                                                       >  1:382843/1‑147 (MQ=255)                                                                               cacTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGtt                                                                                                                <  3:213268/105‑1 (MQ=255)                                                                               cacTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGtt                                                                                                                >  4:213268/1‑105 (MQ=255)                                                                                acTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAagag                                                                                                                              >  2:212403/1‑90 (MQ=255)                                                                                acTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAagag                                                                                                                              <  1:212403/90‑1 (MQ=255)                                                                                acTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTAtt                                                                                     >  2:165458/1‑131 (MQ=255)                                                                                acTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTAtt                                                                                     <  1:165458/131‑1 (MQ=255)                                                                                acTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCAACTGGTTATAACCGGTATTGCTATCCTCATAACCAt                                                                    >  2:108150/1‑148 (MQ=255)                                                                                  tAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTg                                                                  >  2:197884/1‑148 (MQ=255)                                                                                  tAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTg                                                                  >  4:57356/1‑148 (MQ=255)                                                                                   aaaccGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACAAGATTACCGTTATAAAACACCTGGTTATAACCGGTATTGCTCTCCTCATAACCATTga                                                                 >  6:254842/4‑147 (MQ=255)                                                                                   aaCCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTa                                                                                       >  7:334527/1‑126 (MQ=255)                                                                                   aaCCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTa                                                                                       <  2:166286/126‑1 (MQ=255)                                                                                   aaCCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTa                                                                                       <  8:334527/126‑1 (MQ=255)                                                                                   aaCCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTa                                                                                       >  1:166286/1‑126 (MQ=255)                                                                                   aaCCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATACCCATTGt                                                                 >  2:70921/1‑148 (MQ=255)                                                                                        gCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAATCATTGTTTAAc                                                            >  1:11268/1‑148 (MQ=255)                                                                                            ggTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGa                                                          <  3:57356/146‑1 (MQ=255)                                                                                            ggTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGa                                                          <  1:197884/146‑1 (MQ=255)                                                                                              ttGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGc                                                                                   <  4:26668/119‑1 (MQ=255)                                                                                              ttGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGc                                                                                   >  3:26668/1‑119 (MQ=255)                                                                                              ttGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAg                                                      <  6:82083/148‑1 (MQ=255)                                                                                              ttGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAg                                                      >  4:112046/1‑148 (MQ=255)                                                                                              ttGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAg                                                      <  2:382843/148‑1 (MQ=255)                                                                                                    cGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGccc                                                <  7:348982/148‑1 (MQ=255)                                                                                                       gAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCTCCGGTATTGCTCTCGTCATAACCATTGTTTACCGACAAGTTGCCCTTg                                             >  6:232487/1‑148 (MQ=255)                                                                                                               tGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGat                                     <  1:9300/148‑1 (MQ=255)                                                                                                               tGTTTAAGGAGATATTTAATGTGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGATACCAGATTACCGTAATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGat                                     <  5:254842/148‑1 (MQ=255)                                                                                                                gTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGACACTGGAACGCTATCGTAAGAGACCAGATTACCTTTATAAACCACCTGGTTATCACCGGTATTGATCTCCTCATAACCATTGTTTAACGACAAGTTGACCTTGCTCTTGata                                    >  1:271531/1‑148 (MQ=255)                                                                                                                 tttAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGata                                    >  2:283122/1‑147 (MQ=255)                                                                                                                 tttAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGata                                    >  7:369428/1‑147 (MQ=255)                                                                                                                       gagaTCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGATATAGCTaa                             <  1:70921/148‑1 (MQ=255)                                                                                                                           tCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTAAAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTTCCCTTGCTCTTCATATAGATAAAGa                          >  6:79300/1‑147 (MQ=255)                                                                                                                            cTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGATATAGCTAAAGAc                         >  8:198918/1‑147 (MQ=255)                                                                                                                              ttAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGa                                      <  8:325671/132‑1 (MQ=255)                                                                                                                              ttAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGa                                      >  7:325671/1‑132 (MQ=255)                                                                                                                               tAATGTGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGAACAGATTAACGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCGCTTGCTCTTGATATAGCTAAAGACa                        <  5:79300/145‑1 (MQ=255)                                                                                                                               tAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGATATAGCTAAAGACaa                       <  7:198918/146‑1 (MQ=255)                                                                                                                               tAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGATATAGCTAAAGACa                        <  1:283122/145‑1 (MQ=255)                                                                                                                                    ggTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGAAAAGATGCCCTTGCTCTTGATATAGCTAAAGACAATTAGTgg                >  4:211590/1‑148 (MQ=255)                                                                                                                                      ttattaAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGATATAGCTAAAGACAATTTGTGGGt              <  4:120973/148‑1 (MQ=255)                                                                                                                                       tattaAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGATATAGCTAAAGACAATTTGTGGGt              >  1:394838/1‑147 (MQ=255)                                                                                                                                            aaaTCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGATATAGCTAAAGACAATTTGTGGGTTAatat        <  1:108150/148‑1 (MQ=255)                                                                                                                                            aaaTCAGACACTGGAACGCTATCGTAAGAGACCATATTACCGTTATAAACCACCTGGTTATCACCGGTATTTCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGATATAGCTAAAGACAATTTGTGGGTTAata         <  5:232487/147‑1 (MQ=255)                                                                                                                                                  gCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGATATAGCTAAAGACAATTTGTGGGTTAATATCTACCt  >  1:255207/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GCTGATGTATTGCGCATTTGTCTTACCAAGGATAATGCGAGCTTTACGCGCTTCTTGCCAGGCCAGCGTATTCGCCCACACTAACCCGCTCGGTTGATGCGTGAAATCCATGTTTAAGGAGATCTTTAATGGGTTATTAAAATCTGCCACTGGAACGCTATCGTAAGAGACCAGATTACCGTTATAAACCACCTGGTTATCACCGGTATTGCTCTCCTCATAACCATTGTTTAACGACAAGTTGCCCTTGCTCTTGATATAGCTAAAGACAATTTGTGGGTTAATATCTACCT  >  NC_000913/1575475‑1575767

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: ‑