Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,858,886	T→C	T250A (ACC→GCC)	yidE ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,858,886	0	T	C	97.3%	272.4 / ‑1.5	75	T250A (ACC→GCC)	yidE	putative transporter
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (35/38);  total (36/39)
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.15e-01

CACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTT  >  NZ_CP009273/3858754‑3859021                                                                                                                                     |                                                                                                                                        cacGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTgg                                                                                                                                   >  1:24903/1‑139 (MQ=255)    gCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCg                                                                                                                                 >  1:812729/1‑138 (MQ=255)        gACGCGCAAATCAGGGCCTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAAGCACCAGGTGCGCATGATGTAAATCCGCTGGCTGACCCACCAGATGCAGGAAATCGCCCAGGGGGATAATGGCATCGGGCGAAgg                                                                                                                             <  2:615068/138‑1 (MQ=255)         aCGCGCAAATCAGTGCCTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGaa                                                                                                                           >  1:1224117/1‑139 (MQ=255)                        ccTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTttctt                                                                                                            <  2:912437/139‑1 (MQ=255)                           ttCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTcgc                                                                                                         >  2:1137214/1‑139 (MQ=255)                            tCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTcg                                                                                                          <  2:1184443/137‑1 (MQ=255)                            tCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTcg                                                                                                          <  2:401617/137‑1 (MQ=255)                            tCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTcg                                                                                                          <  2:340247/137‑1 (MQ=255)                               tGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTcgccaga                                                                                                     <  1:280199/139‑5 (MQ=255)                                    aGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCa                                                                                                   >  1:599707/1‑136 (MQ=255)                                    aGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCa                                                                                                   <  2:599707/136‑1 (MQ=255)                                      cGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTcg                                                                                                          <  1:214548/127‑1 (MQ=255)                                      cGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTcg                                                                                                          >  2:214548/1‑127 (MQ=255)                                      cGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGAcgc                                                                                              >  2:887605/1‑139 (MQ=255)                                      cGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGAcgc                                                                                              >  2:422046/1‑139 (MQ=255)                                      cGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGAcgc                                                                                              >  1:171088/1‑139 (MQ=255)                                      cGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGAcgc                                                                                              >  2:258320/1‑139 (MQ=255)                                                     ctGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTg                                                                               >  2:777901/1‑139 (MQ=255)                                                     ctGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTg                                                                               >  1:1264497/1‑139 (MQ=255)                                                     ctGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTg                                                                               >  2:923153/1‑139 (MQ=255)                                                     ctGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTg                                                                               >  2:564447/1‑139 (MQ=255)                                                        aCCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCg                                                                            <  2:560717/139‑1 (MQ=255)                                                            aTCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGt                                                                        >  2:456743/1‑139 (MQ=255)                                                             tCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGtt                                                                       >  2:182146/1‑139 (MQ=255)                                                               aCCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTg                                                                      <  2:1140684/138‑1 (MQ=255)                                                                ccAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGtt                                                                       >  1:505499/1‑136 (MQ=255)                                                                ccAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTgag                                                                    <  2:1224117/139‑1 (MQ=255)                                                                   gTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAAt                                                                 <  2:975543/139‑1 (MQ=255)                                                                    ttGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATc                                                                >  1:797651/1‑139 (MQ=255)                                                                       cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCAt                                                                                                                                       >  1:738039/1‑65 (MQ=255)                                                                       cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCAt                                                                                                                                       <  2:738039/65‑1 (MQ=255)                                                                       cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCg                                                                            <  2:77367/124‑1 (MQ=255)                                                                       cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCg                                                                            >  1:77367/1‑124 (MQ=255)                                                                       cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCgg                                                              >  1:1069871/1‑138 (MQ=255)                                                                       cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCgg                                                              >  1:1065120/1‑138 (MQ=255)                                                                       cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCgg                                                              >  2:737005/1‑138 (MQ=255)                                                                       cgcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCgg                                                              >  2:643462/1‑138 (MQ=255)                                                                       cgcAGTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGt                                                             <  2:771943/139‑1 (MQ=255)                                                                        gcATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTa                                                            <  2:583229/139‑1 (MQ=255)                                                                          aTTATGTAAATCCGGTGGCTGACCCACCAGAAGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACa                                                          <  2:1187877/139‑1 (MQ=255)                                                                          aTTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACa                                                          <  1:860642/139‑1 (MQ=255)                                                                                tAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCttt                                                     <  1:777901/138‑1 (MQ=255)                                                                                tAAATACGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCAGCTGGCGAAGGAACTTTTAGGGTTGCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTa                                                    <  2:545782/139‑1 (MQ=255)                                                                                   aTCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCg                                                                            <  2:576668/112‑1 (MQ=255)                                                                                   aTCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCg                                                                            >  1:576668/1‑112 (MQ=255)                                                                                      cGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGc                                               >  1:17384/1‑138 (MQ=255)                                                                                      cGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCt                                              <  2:1069871/139‑1 (MQ=255)                                                                                        ctggctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTa                                             >  2:824956/1‑138 (MQ=255)                                                                                        ctggctgACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTa                                             >  1:168427/1‑138 (MQ=255)                                                                                               aCCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCAt                                       <  2:505499/137‑1 (MQ=255)                                                                                               aCCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCAt                                       <  2:752614/137‑1 (MQ=255)                                                                                               aCCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCAt                                       <  2:1065120/137‑1 (MQ=255)                                                                                               aCCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATg                                      <  1:107515/138‑1 (MQ=255)                                                                                               aCCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATg                                      <  1:211651/138‑1 (MQ=255)                                                                                               aCCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTGTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATg                                      <  2:1063998/138‑1 (MQ=255)                                                                                                     caGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTaaa                                           <  1:240514/127‑1 (MQ=255)                                                                                                     caGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTaaa                                           >  2:240514/1‑127 (MQ=255)                                                                                                     caGATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTa                               <  1:1179181/139‑1 (MQ=255)                                                                                                       gATGCAGCAAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAgg                             >  2:1222665/1‑139 (MQ=255)                                                                                                          gcagcaAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTaa                                                   <  1:1206206/114‑1 (MQ=255)                                                                                                          gcagcaAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTaa                                                   >  2:1206206/1‑114 (MQ=255)                                                                                                              caAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtctc                       >  1:1038377/1‑138 (MQ=255)                                                                                                              caAATCGCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTtctc                       >  2:1108070/1‑138 (MQ=255)                                                                                                               aaaTCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTcg                                                                                                          <  2:241574/54‑1 (MQ=255)                                                                                                               aaaTCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTcg                                                                                                          >  1:241574/1‑54 (MQ=255)                                                                                                                    gCCCAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAAcgcg                 <  1:1108070/138‑1 (MQ=255)                                                                                                                      ccAGTTGGATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAat              <  2:168427/139‑1 (MQ=255)                                                                                                                            ggATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGa         <  2:1130609/138‑1 (MQ=255)                                                                                                                            ggATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGa         <  1:887605/138‑1 (MQ=255)                                                                                                                            ggATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGa         <  1:464100/138‑1 (MQ=255)                                                                                                                            ggATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGa         <  1:1222665/138‑1 (MQ=255)                                                                                                                            ggATAATGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGa         <  1:1043856/138‑1 (MQ=255)                                                                                                                                  tGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCtt  >  1:914170/1‑139 (MQ=255)                                                                                                                                  tGGCATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCtt  >  1:1249255/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                        CACGCTCGACGCGCAAATCAGTGCCTTTCGTGGACAGCGAAGTATCGACCTCCTGACCAATCACCAGTTGCGCATTATGTAAATCCGCTGGCTGACCCACCAGATGCAGCAAATCGCCCAGTTGGATAATGGTATCTGGCGAAGGAACTTTTAGGGTTTCTTCGCGTTTCAGACGCGAGCAGATAATTTTGTCGCCGTTGAGAATCGGTACATCTTTAATGGCTAAATCATGCAGGTTAGGGTTCTCAACGCGAATATTGATAGTCTT  >  NZ_CP009273/3858754‑3859021

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