Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	858,664	G→T	7.6%	S133* (TCA→TAA)	ybiU ←	DUF1479 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	858,664	0	G	T	7.6%	228.6 / 3.9	80	S133* (TCA→TAA)	ybiU	DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base G (42/31);  new base T (2/5);  total (44/36)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACC  >  NC_000913/858521‑858806                                                                                                                                                |                                                                                                                                                tCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAgcg                                                                                                                                            <  1:642206/149‑1 (MQ=255)   acaTCGGGGTAAAACCATTGCTTTCCATCACTTTCAATTCCCCCCAGCCGATTGACAACCGCCTGCCCATTCGCAATTTCTTCCCTCGGCCGGCCCGCCATTGGCGCCTGCTACCAGTAGTGGGGGAAAACCTCGGGCCGTGACgcgcg                                                                                                                                          <  2:1017559/149‑2 (MQ=255)     aTCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTgc                                                                                                                                        >  1:893302/1‑148 (MQ=255)     aTCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGa                                                                                                                                        <  2:1249905/149‑1 (MQ=255)      tCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAg                                                                                                                                       >  2:1312110/1‑149 (MQ=255)      tCCGGGTTAAACCATTGCGTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCGGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAg                                                                                                                                       <  2:195113/149‑1 (MQ=255)       ccGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGc                                                                                                                                      <  2:1485837/149‑1 (MQ=255)       ccGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCttcac                                                                                                                                      <  2:729816/149‑5 (MQ=255)           gTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTcc                                                                                                                                  <  2:713819/149‑1 (MQ=255)            ttAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCGTGCATTTGCGCCTGCGACCAGTTGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTccc                                                                                                                                 >  2:295045/1‑149 (MQ=255)            ttAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTccc                                                                                                                                 <  2:551152/149‑1 (MQ=255)                aCCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCgaga                                                                                                                             <  2:918710/149‑3 (MQ=255)                aCCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCgaga                                                                                                                             <  2:1531828/149‑3 (MQ=255)                aCCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCgaga                                                                                                                             <  2:1566739/149‑3 (MQ=255)                  cATTGCTTTCCATCACTTTCAAATGTCCCCAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCgaagaa                                                                                                                           <  2:1588270/149‑1 (MQ=255)                    ttGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAAt                                                                                                                         >  2:1562601/1‑149 (MQ=255)                     tGCTTTCCATCACTTTCAATTGTCCACAGACGATTGAGAAACGAGTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCgagatttac                                                                                                                        <  2:853253/149‑8 (MQ=38)                      gCTTTCCATCACTTTCAACTGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTa                                                                                                                       <  1:249850/149‑1 (MQ=255)                       cTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTAt                                                                                                                      >  2:1313206/1‑149 (MQ=255)                       cTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTAt                                                                                                                      >  2:723597/1‑149 (MQ=255)                        tttCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATc                                                                                                                     >  1:429860/1‑149 (MQ=255)                           ccATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGcc                                                                                                                  >  2:340271/1‑149 (MQ=255)                           ccATCACTTTCAAATGTCAACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTTAAGCGCTGAGCGTCCCGATGAAATTATCGcc                                                                                                                  >  2:345977/1‑149 (MQ=255)                            cATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTGGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGATGCGCTGAGCGTCCCGAAGAAATTATCGCC‑g                                                                                                                >  2:797183/1‑149 (MQ=255)                             aTCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑gg                                                                                                               >  2:484618/1‑149 (MQ=255)                                     cAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTg                                                                                                       >  2:1272885/1‑149 (MQ=255)                                      aaaTGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGt                                                                                                      >  2:771367/1‑149 (MQ=255)                                        aTGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAg                                                                                                    >  1:1516669/1‑149 (MQ=255)                                          gTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGAc                                                                                                  <  2:32299/149‑1 (MQ=255)                                             cacaGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGATGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCGTTGTAGACCTc                                                                                               >  2:801953/1‑149 (MQ=255)                                              acaGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGTAAAATCTCGGTACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCg                                                                                              <  1:1297522/149‑1 (MQ=255)                                               caGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGACCTTTGTAGACCTCGt                                                                                             >  2:522843/1‑149 (MQ=255)                                                aGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTc                                                                                            >  1:860473/1‑149 (MQ=255)                                                           aaaCGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGCCCTCGTCAAAGCGGTTGc                                                                                 >  2:981350/1‑149 (MQ=255)                                                           aaaCGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGc                                                                                 >  1:1259071/1‑149 (MQ=255)                                                            aaCGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCg                                                                                >  1:1391928/1‑149 (MQ=255)                                                               gACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTc                                                                             >  1:1118488/1‑149 (MQ=255)                                                                 cTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGAGGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCa                                                                           <  1:341409/149‑1 (MQ=255)                                                                   gcgcATTCGCCATTTCTTCACTCTGGCGGGCCTTCATTTGCGCCTGCGACCCGTAGAGGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGa                                                                         <  1:715202/149‑1 (MQ=255)                                                                         tCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTATATCTCGGGACGTGAAGCGCTGAGCGTCCCGATGAATTTATCGCC‑GGGGCCTTGGTAGCCCTCGTCAAAGCGGTTGCGGTCCAGATAGTcc                                                                   >  2:210715/1‑149 (MQ=255)                                                                         tCGCCATTTCTTCACTCGGGCGGTCCTGCTTTTGCGCCTGCGACCTGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCAGTCCAGATAGTcc                                                                   <  1:78152/149‑1 (MQ=255)                                                                          cGCCATTTATTCACTCTGGCGGGCCTGCATTTGCGCCTGCTACCATTAGATGTGTAAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCa                                                                  <  1:1562601/149‑1 (MQ=255)                                                                           gCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAACGCGGTTGCGGTCCAGATAGTCCAg                                                                 >  2:968438/1‑149 (MQ=255)                                                                           gCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCAGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCACGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAg                                                                 >  2:149910/1‑149 (MQ=255)                                                                            ccATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTa                                                                                                     <  1:92629/112‑1 (MQ=255)                                                                            ccATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTa                                                                                                     >  2:92629/1‑112 (MQ=255)                                                                               tttCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGACGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATc                                                             <  1:772260/149‑1 (MQ=255)                                                                               tttCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGACGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATc                                                             <  1:601545/149‑1 (MQ=255)                                                                               tttCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATc                                                             >  1:351453/1‑149 (MQ=255)                                                                                  cttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGAc                                                          >  1:1029002/1‑149 (MQ=255)                                                                                  attcCCTCTGGCGGGCCGGCTTTGGCGCCTGCGACCTGTAGAGGGTTTAAATCGCGGGACGTGAAGCGCTGAGCGCCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGAc                                                          <  1:1312110/148‑1 (MQ=255)                                                                                   ttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACt                                                         >  2:1268078/1‑149 (MQ=255)                                                                                   ttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACt                                                         >  1:1378440/1‑149 (MQ=255)                                                                                      aCTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGAt                                                      <  2:848820/149‑1 (MQ=255)                                                                                      aCTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGAt                                                      >  2:531689/1‑149 (MQ=255)                                                                                           ggCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAg                                                 >  1:723420/1‑149 (MQ=255)                                                                                             cGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGcc                                               <  2:463973/149‑1 (MQ=255)                                                                                                gCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAg                                            >  1:773641/1‑149 (MQ=255)                                                                                                 ccTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTTAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgc                                           <  2:1484224/149‑1 (MQ=255)                                                                                                  cTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgcg                                          >  1:1396621/1‑149 (MQ=255)                                                                                                  cTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGGTCCGGTCCAGATAGTCCAGCATCGCCTGATCCCAGCCTAgcg                                          >  1:1072050/1‑149 (MQ=255)                                                                                                   tGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgcgc                                         >  1:623094/1‑149 (MQ=255)                                                                                                   ggCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgcgc                                         <  2:903206/148‑1 (MQ=255)                                                                                                                   acccttAGATGGGGCAAATCTCGTGACGTGAAGCGCTGACCGTCCCGAAGAAATTATCGCC‑GGGGCCCTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGTGGaaa                         <  1:103900/144‑1 (MQ=255)                                                                                                                       ttAGTTGGGGTAACTCTCGGTCCGTGCAGCGCTGAGGGTCCCGAAGAAATTATCGCC‑GGGGCCTGTGTAGACCTCGTCAAAGCGGGTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGc                     <  1:981350/148‑1 (MQ=255)                                                                                                                       gTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGc                     >  2:969963/1‑149 (MQ=255)                                                                                                                         aGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCt                   <  1:771367/149‑1 (MQ=255)                                                                                                                             ggggAAAATCTCGGGACGCGACGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTa               <  1:797183/149‑1 (MQ=255)                                                                                                                               ggTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCAAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATc             >  1:312936/1‑149 (MQ=255)                                                                                                                                 tAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCAc           >  1:396066/1‑149 (MQ=255)                                                                                                                                 tAAAGCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCAc           <  1:918367/149‑1 (MQ=255)                                                                                                                                    tcttcGGGACGCTCAGCGCTTCACGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCAAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGc        >  2:140030/4‑149 (MQ=255)                                                                                                                                    aTCACGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGATGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCTGGGGAAATGCCCTTTTATCACCGc        >  2:679226/1‑149 (MQ=255)                                                                                                                                      ctcGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGcac      <  1:1272885/149‑1 (MQ=255)                                                                                                                                       tcGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCCGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGcaca     >  1:271482/1‑149 (MQ=255)                                                                                                                                       tcGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGcaca     >  2:1356892/1‑149 (MQ=255)                                                                                                                                        cGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACaa    <  2:860473/149‑1 (MQ=255)                                                                                                                                         tggACGTGAAGCGATGAGCGTCCCGAAGAAATTATCGCCGGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCAAGCATCGACTGATCCCAGCCTATCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACaa    >  2:780858/2‑149 (MQ=255)                                                                                                                                         gggACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAAc   >  1:777859/1‑149 (MQ=255)                                                                                                                                          ggACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAAcc  <  2:1391928/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                TCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCC‑GGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACC  >  NC_000913/858521‑858806

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