Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,100,032	C→T	C320C (TGC→TGT)	nupG →	nucleoside permease NupG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,100,032	0	C	T	94.8%	260.6 / 2.6	77	C320C (TGC→TGT)	nupG	nucleoside permease NupG
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (49/24);  total (51/26)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

TACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGCGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCGGCTGTATCCTCGGCGGCATCGT  >  NZ_CP009273/3099904‑3100163                                                                                                                                 |                                                                                                                                    tACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGAct                                                                                                                           >  1:796166/1‑139 (MQ=255) tACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGCGCATTCGAct                                                                                                                           >  1:815406/1‑139 (MQ=255)           aGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACAt                                                                                                                 >  1:284130/1‑139 (MQ=255)             aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGCGCATTCGACTTCTTCAACAtct                                                                                                               <  1:874577/139‑1 (MQ=255)               cGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATctct                                                                                                             >  2:381338/1‑139 (MQ=255)                gTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTg                                                                                                            >  2:984855/1‑139 (MQ=255)                gTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTg                                                                                                            >  2:527591/1‑139 (MQ=255)                gTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTg                                                                                                            >  1:989262/1‑139 (MQ=255)                gTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTg                                                                                                            >  2:1083620/1‑139 (MQ=255)                  aaTGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGt                                                                                                          >  1:1132135/1‑139 (MQ=255)                          tCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGtt                                                                                                  >  1:55853/1‑139 (MQ=255)                                   tgGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGGACTGTACTGCTGGTACTGTCGATGATCGGTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGaaa                                                                                          >  2:1131946/1‑138 (MQ=255)                                        tGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTAGTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAaagaag                                                                                    >  1:212942/1‑139 (MQ=255)                                        tGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGAc                                                                                                                            >  1:962595/1‑99 (MQ=255)                                        tGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGAc                                                                                                                            <  2:962595/99‑1 (MQ=255)                                        tGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAaagaag                                                                                    >  1:325614/1‑139 (MQ=255)                                        tGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAaagaag                                                                                    >  1:281958/1‑139 (MQ=255)                                               tGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGccc                                                                             >  1:1091650/1‑139 (MQ=255)                                               tGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGccc                                                                             >  1:670526/1‑139 (MQ=255)                                               tGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGccc                                                                             >  2:430686/1‑139 (MQ=255)                                               tGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGccc                                                                             >  1:522779/1‑139 (MQ=255)                                                 cGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCgg                                                                           <  2:522779/139‑1 (MQ=255)                                                      tGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAAt                                                                       <  1:430686/138‑1 (MQ=255)                                                        cgcTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTcg                                                                    <  1:984855/139‑1 (MQ=255)                                                             tttGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCgcgc                                                                 >  2:61501/1‑137 (MQ=255)                                                               tGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAg                                                              <  2:281958/138‑1 (MQ=255)                                                                    aCGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGca                                                          <  2:55853/137‑1 (MQ=255)                                                                     cGGCGACCCGACTCCGTTCGGTAGTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGcaca                                                        <  2:212942/138‑1 (MQ=255)                                                                     cGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGcaca                                                        <  1:61501/138‑1 (MQ=255)                                                                     cGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGcac                                                         >  1:665534/1‑137 (MQ=255)                                                                       gCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAgg                                                     >  2:454600/1‑139 (MQ=255)                                                                       gCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAgg                                                     >  2:266273/1‑139 (MQ=255)                                                                       gCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAgg                                                     >  1:377984/1‑139 (MQ=255)                                                                         gACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGa                                                   <  2:377984/139‑1 (MQ=255)                                                                         gACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGa                                                   <  1:454600/139‑1 (MQ=255)                                                                         gACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGa                                                   <  2:1091650/139‑1 (MQ=255)                                                                          aCCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGAt                                                  >  1:803554/1‑139 (MQ=255)                                                                           cccGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCa                                                               <  1:1245092/125‑1 (MQ=255)                                                                           cccGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCa                                                               >  2:1245092/1‑125 (MQ=255)                                                                           cccGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATg                                                 <  1:681417/139‑1 (MQ=255)                                                                                     ttCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACAtct                                                                                                               >  2:66235/1‑67 (MQ=255)                                                                                     ttCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACAtct                                                                                                               <  1:66235/67‑1 (MQ=255)                                                                                     ttCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCtgatga                                       >  1:793134/1‑139 (MQ=255)                                                                                     ttCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCtgatga                                       >  1:456696/1‑139 (MQ=255)                                                                                     ttCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCtgatga                                       >  2:1246606/1‑139 (MQ=255)                                                                                            ctgtactgCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTa                                 <  2:284130/138‑1 (MQ=255)                                                                                               tactgCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACgg                             >  1:293053/1‑139 (MQ=255)                                                                                               tactgCTGGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACgg                             >  2:188994/1‑139 (MQ=255)                                                                                                 ctgctgGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGc                            <  1:187745/138‑1 (MQ=255)                                                                                                 ctgctgGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGc                            <  2:486300/138‑1 (MQ=255)                                                                                                    ctgGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCg                        >  2:668141/1‑139 (MQ=255)                                                                                                    ctgGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCg                        >  2:894663/1‑139 (MQ=255)                                                                                                    ctgGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCg                        >  1:933648/1‑139 (MQ=255)                                                                                                    ctgGTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCg                        >  2:751957/1‑139 (MQ=255)                                                                                                      ggTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCGGc                      >  2:1111073/1‑139 (MQ=255)                                                                                                      ggTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCGGc                      >  2:528909/1‑139 (MQ=255)                                                                                                       gTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTa                                 <  2:988067/127‑1 (MQ=255)                                                                                                       gTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTa                                 >  1:988067/1‑127 (MQ=255)                                                                                                       gTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCtt                          >  1:952828/1‑134 (MQ=255)                                                                                                       gTACTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCtt                          <  2:952828/134‑1 (MQ=255)                                                                                                         aCTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCgcgc                                                                 <  2:326829/93‑1 (MQ=255)                                                                                                         aCTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCgcgc                                                                 >  1:326829/1‑93 (MQ=255)                                                                                                         aCTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGcaca                                                        <  2:650566/102‑1 (MQ=255)                                                                                                         aCTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGcaca                                                        >  1:650566/1‑102 (MQ=255)                                                                                                         aCTGTCGATGATCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCGGCTg                    >  2:150061/1‑138 (MQ=255)                                                                                                               gatgatCGTTTACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCGGCTGTATCCt              <  1:150061/138‑1 (MQ=255)                                                                                                                  gatCGTTTACGGTTGCGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGaaaaa                                                                                        >  1:724413/1‑61 (MQ=255)                                                                                                                  gatCGTTTACGGTTGCGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGaaaaa                                                                                        <  2:724413/61‑1 (MQ=255)                                                                                                                       tttACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCGGCTGTATCCTCGGCGGCAt     >  2:658462/1‑139 (MQ=255)                                                                                                                       tttACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCGGCTGTATCCTCGGCGGCAt     >  1:1059222/1‑139 (MQ=255)                                                                                                                       tttACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCGGCTGTATCCTCGGCGGCAt     >  1:1281083/1‑139 (MQ=255)                                                                                                                       tttACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCGGCTGTATCCTCGGCGGCAt     >  1:667587/1‑139 (MQ=255)                                                                                                                       tttACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCGGCTGTATCCTCGGCGGCAt     >  2:1187884/1‑139 (MQ=255)                                                                                                                       tttACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCGGCTGTATCCTCGGCGGCAt     >  2:465685/1‑139 (MQ=255)                                                                                                                         tACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCa                                                               >  1:650219/1‑79 (MQ=255)                                                                                                                         tACGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCa                                                               <  2:650219/79‑1 (MQ=255)                                                                                                                           cGGTTGTGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCGGCTGTATCCTCGGCGGCATCGt  <  1:674920/138‑1 (MQ=255)                                                                                                                                 |                                                                                                                                    TACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTTTACGGTTGCGCATTCGACTTCTTCAACATCTCTGGTTCGGTGTTTGTCGAAAAAGAAGTTAGCCCGGCAATTCGCGCCAGTGCACAAGGGATGTTCCTGATGATGACTAACGGCTTCGGCTGTATCCTCGGCGGCATCGT  >  NZ_CP009273/3099904‑3100163

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