Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	756,598	T→A	5.8%	intergenic (‑86/‑127)	bacA ← / → USA300HOU_0707	undecaprenol kinase/ABC superfamily ATP binding cassette transporter, ABC protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	756,598	0	T	A	5.8%	194.7 / 4.8	69	intergenic (‑86/‑127)	bacA/USA300HOU_0707	undecaprenol kinase/ABC superfamily ATP binding cassette transporter, ABC protein
Reads supporting (aligned to +/- strand):  ref base T (37/28);  new base A (2/2);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-01

TAATCCTTCGACGACCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAAC  >  CP000730/756459‑756738                                                                                                                                            |                                                                                                                                             tagtCCTTCGACGACCCCTAAGATAATACCTTTAATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTa                                                                                                                                             <  2:427943/138‑1 (MQ=255) tAATCCTTCGACGACCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATCCTa                                                                                                                                             >  2:583347/1‑141 (MQ=255)  aaTCCTTCGACGACCCCTAGGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAg                                                                                                                                            >  1:497573/1‑141 (MQ=255)  aaTCCTTCGACGACCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAg                                                                                                                                            >  1:400104/1‑141 (MQ=255)  aaTCCTTCGACGACCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGAAATTCTat                                                                                                                                            >  2:634356/1‑140 (MQ=255)    tCCTTCGACGACCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGtt                                                                                                                                          <  2:117681/141‑1 (MQ=255)    tCCTTCGACGACCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGtt                                                                                                                                          <  2:172276/141‑1 (MQ=255)        tCGACGACCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTGATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATACTAGTTataa                                                                                                                                      >  1:46528/1‑141 (MQ=255)        tCGACGACCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTataa                                                                                                                                      >  1:225166/1‑141 (MQ=255)           acggcCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATGATAc                                                                                                                                   <  2:7966/141‑1 (MQ=255)           acgacCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATAc                                                                                                                                   <  1:577948/141‑1 (MQ=255)            cgacCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACt                                                                                                                                  >  1:534259/1‑141 (MQ=255)             gacCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCACTACTCAATTTAATTTAGTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTa                                                                                                                                 <  2:236976/141‑1 (MQ=255)             cacCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTa                                                                                                                                 <  1:309825/140‑1 (MQ=255)               ccccTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCAAGTTATAATACTAGa                                                                                                                               >  1:259637/1‑141 (MQ=255)               ccccTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCAAGTTATAATACTAGa                                                                                                                               >  1:516059/1‑141 (MQ=255)                 ccTAATATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACt                                                                                                                             <  2:125876/141‑1 (MQ=255)                         aataCCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTATAGATGt                                                                                                                     <  2:367418/141‑1 (MQ=255)                           taCCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGtat                                                                                                                   >  1:221286/1‑141 (MQ=255)                           taCCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGtat                                                                                                                   <  1:13340/141‑1 (MQ=255)                               ttttATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTatatat                                                                                                               <  2:331092/141‑1 (MQ=255)                               ttttATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGAAATGCTAGTTATAATACTAGACTAAAGATGTatatat                                                                                                               <  1:460881/141‑1 (MQ=255)                                            aTGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTa                                                                                                  <  1:292983/141‑1 (MQ=255)                                                 aaaCATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCCTTTGGTGAGGTGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGt                                                                                             <  1:624770/141‑1 (MQ=255)                                                    cATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTAACATATCATGATAGCATATTCATTTAAAGACATGCTAGTTGTAATACTAGACTAAAGATGTATATATTAATTTTCTTTTACTTGTaaa                                                                                          >  1:235609/1‑141 (MQ=255)                                                      tAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAAc                                                                                        >  1:247747/1‑141 (MQ=255)                                                      tAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAAc                                                                                        >  1:237876/1‑141 (MQ=255)                                                      tAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGAATAAAGATGTATATATTCATTTTCTTTTACATGTaaaat                                                                                        >  2:609232/1‑140 (MQ=255)                                                                   ttCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATCTTTTAt                                                                           <  2:367915/141‑1 (MQ=255)                                                                       ttACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAg                                                                       >  1:115040/1‑141 (MQ=255)                                                                       ttACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAg                                                                       >  1:538697/1‑141 (MQ=255)                                                                          cTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTa                                                                    >  2:175853/1‑141 (MQ=255)                                                                           tCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTAt                                                                   >  1:96824/1‑141 (MQ=255)                                                                                 taatttaTTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCAAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAAtt                                                             <  2:259637/141‑1 (MQ=255)                                                                                 taatttaTTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCAAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAAtt                                                             <  2:516059/141‑1 (MQ=255)                                                                                  aatttaTTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAAttt                                                            <  2:117922/141‑1 (MQ=255)                                                                                  aatttaTTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAAttt                                                            <  1:377054/141‑1 (MQ=255)                                                                                  aatttaTTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAAttt                                                            <  2:596465/141‑1 (MQ=255)                                                                                   atttatttAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGATCTTTTTAAtttt                                                           >  2:577532/1‑140 (MQ=255)                                                                                   atttatttAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTg                                                           <  1:175295/141‑1 (MQ=255)                                                                                    tttatttaAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGa                                                          >  1:309667/1‑141 (MQ=255)                                                                                     ttatttaAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGAt                                                         >  2:537512/1‑141 (MQ=255)                                                                                          taAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTaa                                                    >  2:224164/1‑141 (MQ=255)                                                                                           aaaTATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAg                                                   >  1:543838/1‑141 (MQ=255)                                                                                           aaaTATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAg                                                   >  1:539228/1‑141 (MQ=255)                                                                                             atatCAAAATTACCATATAATGATAGCATATTCATTTAAAGACATGCTCGTTATAATACTAGACTAAAGATGTTTATATTCATTTTCTTTTATATGTAATAATACAATATTTTTTTGAGCTATTTAATTTGATTTTAAGGa                                                 >  2:231776/1‑141 (MQ=255)                                                                                              tatCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGaa                                                >  1:570697/1‑141 (MQ=255)                                                                                              tatCAAAATTACCATATCATGATAGCATATTCATTTAAAGAAATGCTAGTTATATTACTAGACTAAAGTTGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTTTTTAATTTGTTTTTTAGGaa                                                >  2:265638/1‑141 (MQ=255)                                                                                                   aaaTTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCt                                           <  2:409406/141‑1 (MQ=255)                                                                                                      ttACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCtttt                                        >  1:249244/1‑141 (MQ=255)                                                                                                        aCCTTATCACAATAGTATATTCATTTAATCACATGCTAGTTATAATACTAGACTAAATATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGACAACCTTTtat                                      <  2:214809/141‑1 (MQ=255)                                                                                                        aCCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTtat                                      >  2:320351/1‑141 (MQ=255)                                                                                                            tatCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTTAGGAAAACCTTTTataata                                  <  1:399072/141‑1 (MQ=255)                                                                                                            tatCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTataata                                  <  2:73706/141‑1 (MQ=255)                                                                                                            tatCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTataata                                  >  1:412416/1‑141 (MQ=255)                                                                                                             ttCATGATAGCATAGTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAg                                 <  2:620815/140‑1 (MQ=255)                                                                                                             ataaTGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTTTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAg                                 <  2:363/138‑1 (MQ=255)                                                                                                                  gATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTa                            >  1:133331/1‑141 (MQ=255)                                                                                                                  gATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGAAAAACCTTTTATCAATGGTTTa                            >  2:327380/1‑141 (MQ=255)                                                                                                                      gCATATTCAGTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCAGTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTg                        <  2:69727/141‑1 (MQ=255)                                                                                                                           ttCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGAttta                   >  2:659556/1‑141 (MQ=255)                                                                                                                             cATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAAt                 >  2:579332/1‑141 (MQ=255)                                                                                                                             cATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAAt                 >  1:329093/1‑141 (MQ=255)                                                                                                                             cATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAAt                 >  1:410071/1‑141 (MQ=255)                                                                                                                                gtAAATACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACAGGTAAAACTACAATATTTTATTGAGCTATTTAAATTGATTTTAAGTTATACGTTTTATAATAGGTTTAGGTGATATAATtgt              <  2:405660/140‑1 (MQ=255)                                                                                                                                    aGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTTAGCTATTTAATTTGATTTGAAGGAAAACCTTTTATAATAGGTTTAGGTGATATACTTGTGaaa          >  2:82335/1‑141 (MQ=255)                                                                                                                                    aGACATGCTAGTTACAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACATTTTATAATAGGTTTAGGTGATATAATTGTGaaa          >  1:347217/1‑141 (MQ=255)                                                                                                                                        aTGCTAGTTATAATACTAGACTAAAGCTTTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTGTATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAAt      <  2:269647/141‑1 (MQ=255)                                                                                                                                          gCTAGTTATAATACTAGACTAAAGATTACTATATTCATTTTCTTTTACATGTAAAACTACAATAAATTATTGCGCTATTTAATTTGATTTTAAGGAAAATCTTTTATAATAGGTTTCGCTAATATAATTGTGAAAAAAtca    >  2:220381/1‑139 (MQ=255)                                                                                                                                           ttAGTTATTATTATAGACTTCCGATGTATTTTTTTCTTTTTTTTTACATGTAAAAGTACAATTGTTTATTGAGTTTGTTTATTTGATTTTAAGGAAGACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTaa   <  2:222522/140‑1 (MQ=255)                                                                                                                                            tAGTTATAATACTAGACTAAAGATGTATATATTGATTTTCCTGTACATGTAAAACGACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTaac  >  1:319593/1‑141 (MQ=255)                                                                                                                                            |                                                                                                                                             TAATCCTTCGACGACCCCTAAGATAATACCTTTTATTAATTCAATGATAAACATAATGTACCCACTTTCATTACTCAATTTAATTTATTTAAATATCAAAATTACCATATCATGATAGCATATTCATTTAAAGACATGCTAGTTATAATACTAGACTAAAGATGTATATATTCATTTTCTTTTACATGTAAAACTACAATATTTTATTGAGCTATTTAATTTGATTTTAAGGAAAACCTTTTATAATAGGTTTAGGTGATATAATTGTGAAAAAATTAAC  >  CP000730/756459‑756738

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