Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,857,867	C→T	9.0%	L169L (TTG→TTA)	USA300HOU_1728 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,857,867	0	C	T	9.0%	161.9 / 2.5	69	L169L (TTG→TTA)	USA300HOU_1728	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (40/21);  new base T (2/5);  total (42/27)
Fisher's exact test for biased strand distribution p-value = 9.69e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGCTTCTTGAATCTCTTGTG  >  CP000730/1857727‑1858006                                                                                                                                             |                                                                                                                                            tcttcttTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTc                                                                                                                                             <  1:308050/141‑1 (MQ=255)     cttTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCCCTTGTGCCTCTTGTGGTGTATCACTCTAat                                                                                                                                         >  2:485372/1‑141 (MQ=255)     cttTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTGCTTTTTTAGCTTCTTGAATTTCTTATGCCTCTTGTGATGTTTACCTCAAat                                                                                                                                         >  1:409702/1‑141 (MQ=17)      ttttATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAatt                                                                                                                                        <  1:638434/141‑1 (MQ=255)        ttATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGATTCTAAAGGATCTCCTTGAATCTCTTGTTCCTCTTGTGATGTATCACTCAAattat                                                                                                                                      >  1:549648/1‑141 (MQ=17)         tATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTGTTGAATCTCTTGTGCGTCTTGTGATGTATCACTCAAattatt                                                                                                                                     >  1:248448/1‑141 (MQ=255)           tCGCTGCTTGTTGTGCTTTCAATGCCACTGATTTTGGTGCGACATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTg                                                                                                                                   >  1:363995/1‑141 (MQ=255)             gctgctTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCTTCTTGTGATGTATCACTCAAATTATTTGCa                                                                                                                                 >  2:353893/1‑141 (MQ=255)                  ttgttgTGCTTTCAATGTTACTGCTTTTTCTTCTTCATTTGTTACAGTTGCATTTTTATCTTTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGc                                                                                                                            <  2:333135/141‑1 (MQ=17)                   tgttttGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGAAACCGCTGAATTAGCGGCGTTTACTGGTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGct                                                                                                                           <  1:517846/141‑1 (MQ=17)                    gttgtGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAATTTATTTGCACTTGctt                                                                                                                          >  1:51137/1‑141 (MQ=255)                    gttgtGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGctt                                                                                                                          >  1:394858/1‑141 (MQ=255)                           tttCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCCTCTTGGAACGCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCttt                                                                                                                   >  1:359751/1‑141 (MQ=255)                              cAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTAt                                                                                                                >  1:151155/1‑141 (MQ=255)                               aaTGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTCAATTATTTGCACTTGCTTCTTCTTTTTTc                                                                                                               >  2:48086/1‑141 (MQ=255)                                   ctaCTGTTTTTGGTTTTTCATTTGTTGCTGCTGCACTTTTGTCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctg                                                                                                           <  2:667298/139‑1 (MQ=17)                                    cACTGCTTTTGGTTCTTCCTTCGGGACTGTGGCACTTTTATCTGTTTCTGCATGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCTTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCgctgc                                                                                                          >  2:539448/1‑141 (MQ=14)                                     aCTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGCTGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTAGCgctgct                                                                                                         >  2:400639/1‑141 (MQ=255)                                      cTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAAACTCTTGTGTTTTTGATCATTTATCACTTAATTTATTTGCACTTGCTTCTTCTTTTATCGCTGCtt                                                                                                        >  2:556446/1‑141 (MQ=14)                                        gCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCttgt                                                                                                      >  1:247041/1‑141 (MQ=255)                                               gTTCTTCATTTGTTACAGCTGCACCTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCtt                                                                                               <  1:194709/141‑1 (MQ=255)                                               gTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATATGCACTTGCTTCTCCTTTTATCGCTGCTTGTTGTGCtt                                                                                               <  1:486830/141‑1 (MQ=255)                                                  cttcATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATGGCTGCTGGTTGTGCTTTCa                                                                                            >  1:753843/1‑141 (MQ=255)                                                   ttcATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTcta                                                                                           >  2:343657/1‑139 (MQ=255)                                                   ttcATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCaa                                                                                           <  1:522477/141‑1 (MQ=255)                                                            tACAGCTGCACTTTTATCTGATTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCGGCTTGTTGTGCTTTCAATGCCACTGc                                                                                  >  1:458210/1‑141 (MQ=255)                                                              cAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCtt                                                                                >  1:650615/1‑141 (MQ=255)                                                              cAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCtt                                                                                >  1:684468/1‑141 (MQ=255)                                                                gCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCtttt                                                                              >  1:575163/1‑141 (MQ=255)                                                                   ccttctgttcctttttCTGCTTGTGCTTCTTTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGAt                                                                           <  2:367204/128‑1 (MQ=11)                                                                      cTTTTATCTGTTTGTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTGATTGCGACTGCTTCTGAtttt                                                                        >  2:638692/1‑139 (MQ=17)                                                                        ttcttCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTTTATATTTCGATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCttt                                                                      <  1:174770/137‑1 (MQ=11)                                                                           ttCTTTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTAtt                                                                   <  2:238180/140‑1 (MQ=255)                                                                            tCTTTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGCCTTTCTTGTGGTTATTATCCTTTTTCACTGAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTAttt                                                                  <  2:731468/141‑1 (MQ=14)                                                                                   cTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAg                                                           <  1:712453/141‑1 (MQ=255)                                                                                        tgtgCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTTTTTATTTGATACAGCcaca                                                      >  2:687367/1‑141 (MQ=255)                                                                                         gtgCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCccac                                                     >  1:137491/1‑141 (MQ=255)                                                                                         gtgCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCacac                                                     >  1:272527/1‑141 (MQ=255)                                                                                            cttcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACttt                                                  <  2:388298/141‑1 (MQ=255)                                                                                            cttcttTTTTAGCTTCTTGAACCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACttt                                                  <  2:430988/141‑1 (MQ=255)                                                                                             ttcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACtttt                                                 >  1:493290/1‑141 (MQ=255)                                                                                              tcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCCATTTACTCTTCTTTTATCGCTGCGTGTTGTACTTTCAATTCCGCTGGTTGTGAGTCTTTTTTTGCTACatcgcacatctca                                                >  2:24131/1‑129 (MQ=9)                                                                                              tcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCCATGACCCCGCTACGCCTGCCATTTTTGATACAACCACACTTTTa                                                >  1:100362/1‑141 (MQ=14)                                                                                              tcttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTa                                                <  1:228853/141‑1 (MQ=255)                                                                                               cttTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACCTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTAt                                               >  2:73449/1‑141 (MQ=255)                                                                                                  ttttAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTAATTTGATACAGCAACACTTTTATCCg                                            >  1:341649/1‑141 (MQ=255)                                                                                                     tAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTTATACAGCCACACTTTTATCCGCtt                                         >  2:640002/1‑141 (MQ=255)                                                                                                     tAGCTTCTTGAACCTCTTGTGCCTCTTGTGATATATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCtt                                         <  2:263968/141‑1 (MQ=255)                                                                                                      aGCTTCTTGAATCTCTTGTGCTTCTGGTTATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTTAATGCCGCTGCTTTTGATTCTTTATTTGATACAGCTACATTTTTATCCGCttt                                        >  2:735871/1‑140 (MQ=17)                                                                                                      aGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACtt                                                                                                                              >  1:680584/1‑55 (MQ=255)                                                                                                      aGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCtt                                                                                               <  1:113475/86‑1 (MQ=255)                                                                                                      aGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTc                                        >  1:383275/1‑141 (MQ=255)                                                                                                        cttcttGAATCTCTTGTGCCTCTTGTGATTTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTg                                      <  1:707121/141‑1 (MQ=255)                                                                                                          tcttGACTCTCTTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCt                                    <  2:385411/141‑1 (MQ=255)                                                                                                                acctctTGTGCCTCTTGTGATGTATCACTTAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCt                              <  1:415654/139‑1 (MQ=255)                                                                                                                 tctcTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCtt                             >  1:514552/1‑141 (MQ=255)                                                                                                                  ctctTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTAACGCACCTTTTTGTTCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTtc                            <  1:175150/141‑1 (MQ=17)                                                                                                                     tttttccTCTTGTTATGTATCACTCAAATTATTCTCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTctct                         <  1:549893/136‑1 (MQ=255)                                                                                                                     ttGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTCATCCGCTTCTGCTTGTGCTTctct                         <  1:485372/141‑1 (MQ=255)                                                                                                                       gtgCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCttt                       >  2:707314/1‑141 (MQ=255)                                                                                                                       gtgCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCttt                       >  2:707714/1‑141 (MQ=255)                                                                                                                           ctctTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGc                   <  2:272527/141‑1 (MQ=255)                                                                                                                              tggtgATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTTTTACCTTGATACAGACCGGCTGTTGTCCGCTTTTGCTTGTGCTTCTGATTTAGcttc                <  1:556446/139‑1 (MQ=14)                                                                                                                                    tGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCGCTGCTTTTTATTCTTTATTTTTTACAGCTACTCTTTTATCTTCTTCTGCTTATGCTTCTCTTTTTTCTTCTTGCAt          >  2:474167/1‑141 (MQ=14)                                                                                                                                     gTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGCTTCTTGAAtc         >  1:535356/1‑141 (MQ=255)                                                                                                                                      tATCCCTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTTTTTGGTACAGCCACACTTTTTTCCGCTTGTGCTTTTGGTTCTTTTTTCGcttctttatttt        >  1:16035/1‑135 (MQ=14)                                                                                                                                       aTCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCGCTACTTTTGATTCTTTTTTTGTTACAGCCACACTTTTATCCGCTTCTGTTTGTGCTTCTCTTTTTGCTTCTTGagtgtt       >  2:133846/1‑136 (MQ=17)                                                                                                                                        tCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGc                                                                                  <  1:269422/65‑1 (MQ=35)                                                                                                                                           cTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTCTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGCTTCTTGAATCTCTtgt   <  1:511333/141‑1 (MQ=255)                                                                                                                                            tCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCGCTGCTTTTGTTTCTTTATTTGTTACAGCCCCACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGCTTCTTGCATCTCTtgtt  >  2:588560/1‑140 (MQ=14)                                                                                                                                             |                                                                                                                                            TCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGGTTCTTCATTTGATACAGCTGCACTTTTATCTGTTTCTGCTTGTGCTTCTTTTTTAGCTTCTTGAATCTCTTGTGCCTCTTGTGATGTATCACTCAAATTATTTGCACTTGCTTCTTCTTTTATCGCTGCTTGTTGTGCTTTCAATGCCACTGCTTTTGATTCTTTATTTGATACAGCCACACTTTTATCCGCTTCTGCTTGTGCTTCTCTTTTAGCTTCTTGAATCTCTTGTG  >  CP000730/1857727‑1858006

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