Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,919,585	T→A	5.5%	G692G (GGT→GGA)	yebT →	MCE domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,919,585	0	T	A	5.5%	225.1 / 2.1	72	G692G (GGT→GGA)	yebT	MCE domain protein
Reads supporting (aligned to +/- strand):  ref base T (35/33);  new base A (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATT  >  NC_000913/1919440‑1919733                                                                                                                                                  |                                                                                                                                                     ccGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCggtggt                                                                                                                                                   <  6:240671/149‑1 (MQ=255)   ggAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTAc                                                                                                                                                 <  3:310333/149‑1 (MQ=255)             cGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCa                                                                                                                                       <  4:384793/149‑1 (MQ=255)             cGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCa                                                                                                                                       <  2:143648/149‑1 (MQ=255)                tCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTg                                                                                                                                    <  8:104335/149‑1 (MQ=255)                   ggATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTc                                                                                                                                 <  3:129948/149‑1 (MQ=255)                    gATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTca                                                                                                                                >  6:127575/1‑149 (MQ=255)                        ccAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTAcgc                                                                                                                                               >  8:32604/1‑130 (MQ=255)                        ccAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTAcgc                                                                                                                                               <  7:32604/130‑1 (MQ=255)                          aaTTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCa                                                                                                                                       >  4:25307/1‑136 (MQ=255)                          aaTTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCa                                                                                                                                       <  3:25307/136‑1 (MQ=255)                          aaTTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCg                                                                                                                           >  7:449664/1‑148 (MQ=255)                                cTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGAGGTACGCGCTTCTCAGTGGTCACACCGCAAAttt                                                                                                                    >  7:322905/1‑149 (MQ=255)                                cTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGAGGTACGCGCTTCTCAGTGGTCACACCGCAAAttt                                                                                                                    >  7:331138/1‑149 (MQ=255)                                 tATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCa                                                                                                                                       <  3:431651/129‑1 (MQ=255)                                 tATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCa                                                                                                                                       >  4:431651/1‑129 (MQ=255)                                   tctcGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCgg                                                                                                                 >  5:396830/1‑149 (MQ=255)                                   tctcGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCg                                                                                                                  <  5:127575/148‑1 (MQ=255)                                      cGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAg                                                                                                              >  4:391871/1‑149 (MQ=255)                                           ttGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCggtg                                                                                                                                                     <  4:8945/105‑1 (MQ=255)                                           ttGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCggtg                                                                                                                                                     >  3:8945/1‑105 (MQ=255)                                           ttGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGc                                                                                                         >  1:524967/1‑149 (MQ=255)                                              atttCGGGCAAATCTAGACGCTGGATCTGATTACCGAACGCAATGAAGTACTGGCAAATGCGGTACTTTTCTCGGATTATGTCGAGACCTTTTATTGTTTTGGTTCGCGTTTCTCAGTGTTAACACCGCAAATTTCGGCAGCTTGCGtt                                                                                                      <  4:417743/146‑1 (MQ=255)                                              atatCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGtt                                                                                                      >  3:508589/1‑149 (MQ=255)                                                 tCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAg                                                                                                   >  1:203527/1‑149 (MQ=255)                                                 tCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAg                                                                                                   >  2:320244/1‑149 (MQ=255)                                                  cGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGa                                                                                                    >  7:170907/1‑147 (MQ=255)                                                  cGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGa                                                                                                    <  8:170907/147‑1 (MQ=255)                                                       aaaTCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCtt                                                                                             >  7:346471/1‑149 (MQ=255)                                                        aaTCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCtt                                                                                             >  2:250419/1‑148 (MQ=255)                                                            cAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTTCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATa                                                                                         <  6:187454/148‑1 (MQ=255)                                                             aGAAGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGGACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACt                                                                                       <  2:214102/149‑1 (MQ=255)                                                                 gCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTAtcc                                                                                   <  1:250419/149‑1 (MQ=255)                                                                          gATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCg                                                                           >  8:388194/1‑148 (MQ=255)                                                                           aTTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCa                                                                                                                                       >  5:247898/1‑87 (MQ=255)                                                                           aTTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGc                                                                             <  8:11458/145‑1 (MQ=255)                                                                           aTTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGt                                                                          <  3:391871/148‑1 (MQ=255)                                                                           aTTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGCTACGCGCTTCGCa                                                                                                                                       >  5:239809/1‑87 (MQ=255)                                                                                         tGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAAc                                                           <  6:396830/149‑1 (MQ=255)                                                                                           aaGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCa                                                         <  3:446070/149‑1 (MQ=255)                                                                                           aaGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCa                                                         <  2:203527/149‑1 (MQ=255)                                                                                              tACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTc                                                               <  5:312069/140‑1 (MQ=255)                                                                                              tACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTc                                                               >  6:312069/1‑140 (MQ=255)                                                                                               aCAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGc                                                      <  8:346471/148‑1 (MQ=255)                                                                                                cAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCg                                                    >  8:301910/1‑149 (MQ=255)                                                                                                 aGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGTCCgg                                                   <  8:449664/149‑1 (MQ=255)                                                                                                 aGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCgg                                                   >  3:292720/1‑149 (MQ=255)                                                                                                     aaaGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGAt                                                                                          >  1:417041/1‑106 (MQ=255)                                                                                                     aaaGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGAt                                                                                          <  2:417041/106‑1 (MQ=255)                                                                                                     aaaGGCGGTACTTTA‑CCGGAATATGTCCACACCTTTGCTCGCGGTGGTACGCGCATCACAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACAATCCTACCGCCGTATATCAACGTCGAACCAGGCCCGGGCa                                               >  4:281236/1‑148 (MQ=255)                                                                                                      aaGGCGGTACTTT‑CGCGGAATATGTCCAGACCTTTGCTCGCTGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAAt                                             <  3:281236/149‑1 (MQ=255)                                                                                                             tACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTcgc                                       <  2:524967/149‑1 (MQ=255)                                                                                                               cTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTcgcc                                      <  1:320244/148‑1 (MQ=255)                                                                                                                   aCCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGa                                  <  7:301910/148‑1 (MQ=255)                                                                                                                     ccGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACtt                               >  7:546206/1‑149 (MQ=255)                                                                                                                       ggAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTg                             <  7:388194/149‑1 (MQ=255)                                                                                                                          atatGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGTGGCAATCCTCGCCGCGACTTTGAAt                          <  4:292720/149‑1 (MQ=255)                                                                                                                            atGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTa                        >  7:119405/1‑149 (MQ=255)                                                                                                                              gTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCggg                                                  >  6:235501/1‑121 (MQ=255)                                                                                                                              gTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCggg                                                  <  5:235501/121‑1 (MQ=255)                                                                                                                               tCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTAc                       >  8:508808/1‑147 (MQ=255)                                                                                                                                ccAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCCGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTAc                       >  1:137110/1‑146 (MQ=255)                                                                                                                                  aGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAagag                  >  7:530858/1‑149 (MQ=255)                                                                                                                                   gACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGAc                                 <  6:383922/133‑1 (MQ=255)                                                                                                                                   gACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGAc                                 >  5:383922/1‑133 (MQ=255)                                                                                                                                      cTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACttt                              >  3:231104/1‑133 (MQ=255)                                                                                                                                      cTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACttt                              <  4:231104/133‑1 (MQ=255)                                                                                                                                      cTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGcca              >  4:122765/1‑149 (MQ=255)                                                                                                                                       tttGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCa                                                                                                 <  4:330716/65‑1 (MQ=255)                                                                                                                                       tttGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCa                                                                                                 >  3:330716/1‑65 (MQ=255)                                                                                                                                           cTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCATCGCCGCGACGTTGAATTACAAGAGGCCACCAt          >  6:425932/1‑148 (MQ=255)                                                                                                                                              gcgGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACt      >  1:247197/1‑149 (MQ=255)                                                                                                                                              gcgGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACt      >  2:360983/1‑149 (MQ=255)                                                                                                                                              gcgGAGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTa        <  8:322905/147‑1 (MQ=255)                                                                                                                                              gcgGAGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTa        <  8:331138/147‑1 (MQ=255)                                                                                                                                                 gtggtACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTgta                           <  4:5346/125‑3 (MQ=255)                                                                                                                                                 gtggtACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTgta                           >  3:5346/1‑123 (MQ=255)                                                                                                                                                  tggtACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGAtt  >  5:381102/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CCGGAAAACTGGCGGTCGGGATGCCAATTCGCTATCTCGGTATTGATATCGGGCAAATCCAGACGCTGGATCTGATTACCGCACGCAATGAAGTACAGGCAAAGGCGGTACTTTACCCGGAATATGTCCAGACCTTTGCTCGCGGTGGTACGCGCTTCTCAGTGGTCACACCGCAAATTTCGGCAGCTGGCGTTGAGCATCTTGATACTATCCTCCAGCCGTATATCAACGTCGAACCAGGCCGGGGCAATCCTCGCCGCGACTTTGAATTACAAGAGGCCACCATTACTGATT  >  NC_000913/1919440‑1919733

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