Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	787,902	G→A	5.5%	R23H (CGC→CAC)	pabC →	possible aminodeoxychorismate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	787,902	0	G	A	5.5%	220.1 / 2.2	72	R23H (CGC→CAC)	pabC	possible aminodeoxychorismate lyase
Reads supporting (aligned to +/- strand):  ref base G (33/35);  new base A (2/2);  total (35/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

TTACAATTGATTCTAAGCCAAAAGATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTATTAAATGCTGAAAGCAAATTTGAA  >  CP000730/787762‑788041                                                                                                                                             |                                                                                                                                            ttACAATTGATTCTAAGCCAAAAGATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCg                                                                                                                                             <  2:629435/141‑1 (MQ=255)       ttGATTCTAAGCCAAAAGATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCACATaaa                                                                                                                                       >  2:628033/1‑141 (MQ=255)       ttGATTCTAAGCCAAAAGATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCACATaaa                                                                                                                                       >  2:627948/1‑141 (MQ=255)           ttCTAAGCCAAAAGATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGt                                                                                                                                   <  1:59937/141‑1 (MQ=255)           ttCTAAGCCAAAAGATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGt                                                                                                                                   <  1:677981/141‑1 (MQ=255)            tCTAAGCCAAAAGATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGtt                                                                                                                                  <  1:606462/141‑1 (MQ=255)            tCTAAGCCAAAAGATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGtt                                                                                                                                  >  2:471888/1‑141 (MQ=255)             cTAAGCCAAAAGATGAAGTGAACGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCACTTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGttt                                                                                                                                 <  2:161286/141‑2 (MQ=255)               aaGCCAAAAGATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGttctt                                                                                                                               >  1:266336/1‑141 (MQ=255)                aGCCAAAAGATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCGATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGttcttt                                                                                                                              >  2:14759/1‑140 (MQ=255)                aGCCAAAAGATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGttcttc                                                                                                                              >  2:210800/1‑141 (MQ=255)                      aaGATGAAGTGAATGAATTTTACGCAAAACCCCAGATTTTGGAGATGTTATAATGCAATGAGTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTACTCGCATCAAATGTTCTACTGAGCg                                                                                                                        >  2:288163/1‑141 (MQ=255)                       aGATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGa                                                                                                                       >  1:5843/1‑141 (MQ=255)                        gATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGAt                                                                                                                      <  1:708134/141‑1 (MQ=255)                        aaTGAAGTTAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGAt                                                                                                                      <  2:304431/140‑1 (MQ=255)                          tGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTa                                                                                                                    >  2:879328/1‑141 (MQ=255)                          tGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTa                                                                                                                    >  2:697838/1‑141 (MQ=255)                          tGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTa                                                                                                                    >  2:159461/1‑141 (MQ=255)                           gAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTaa                                                                                                                   >  2:288126/1‑141 (MQ=255)                              gTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACt                                                                                                                <  2:649441/141‑1 (MQ=255)                              gTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACt                                                                                                                <  1:911235/141‑1 (MQ=255)                                    gaatTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAAt                                                                                                          <  2:757969/141‑1 (MQ=255)                                     aagtttACGCAAAAACAAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAAGTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAAtt                                                                                                         <  1:186341/138‑1 (MQ=255)                                       ttttACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTtg                                                                                                       >  2:279705/1‑141 (MQ=255)                                           aCGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTtgatga                                                                                                   <  1:461589/141‑1 (MQ=255)                                             gccaaaaACAAGATTTTGGAGATGATATAATGCTATTACTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAAc                                                                                                 >  2:893389/4‑141 (MQ=255)                                                aaaaCCAAGAATTTGGAGATGTTATAATGCAATTATTTGAAACAATAAATAATGTTAATGGATATATTACTAGGCTTTCTTTGGATACTAATCGCATAAAATGATCTTTTGAGCGATTAAACTTTAGATTTGATGAAcata                                                                                              >  2:411988/1‑140 (MQ=255)                                                 aaaCCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAAcatgc                                                                                             >  1:868324/1‑141 (MQ=255)                                                  aaCCAAGCTTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAAcatgca                                                                                            <  2:73395/141‑1 (MQ=255)                                                  aaCCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAAcatgca                                                                                            <  1:242288/141‑1 (MQ=255)                                                  aaCCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAAcatgca                                                                                            <  2:781432/141‑1 (MQ=255)                                                      aaGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGc                                                                                        <  2:631500/141‑1 (MQ=255)                                                        gATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGa                                                                                      >  1:439760/1‑141 (MQ=255)                                                           tttGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAAATTTAAATTTGATAAAcatgcctgtttataa                                                                                   >  2:888685/1‑134 (MQ=255)                                                            ttGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAAGCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATg                                                                                  <  2:5843/141‑1 (MQ=255)                                                              ggAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGaa                                                                                >  2:815052/1‑141 (MQ=255)                                                               gagaTGTTATAATGCAATTAGTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAAt                                                                               <  1:846035/141‑1 (MQ=255)                                                                  aTGTAATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTaa                                                                            <  2:254798/141‑1 (MQ=255)                                                                    gTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGCCATACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGACTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAAc                                                                          <  2:518075/141‑1 (MQ=255)                                                                         aaTGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGt                                                                     <  1:581085/141‑1 (MQ=255)                                                                          aTGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTa                                                                    <  2:96274/141‑1 (MQ=255)                                                                              aaTTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTaacaa                                                                <  2:223699/141‑1 (MQ=255)                                                                              aaTTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAAGTAAACGATGTaacaa                                                                <  1:699485/141‑1 (MQ=255)                                                                              aaTTATTTGAAACAATGAAAATTGATAATGGACATAGCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTaacaa                                                                <  2:831891/141‑1 (MQ=255)                                                                               attattTGAAACAATGAAAATTGATAATGGCCATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTaacaac                                                               <  1:510537/141‑1 (MQ=255)                                                                                 tattTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAAcaacaa                                                             <  1:664631/141‑1 (MQ=255)                                                                                 tattTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAAcaacaa                                                             <  2:537432/141‑1 (MQ=255)                                                                                  attTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACaaa                                                            >  2:607761/1‑141 (MQ=255)                                                                                    ttGAAACAGTGAAAATAGATGATGGACTTATCCCTAGACTTACTTATCATACTAATAGAATACAATGTACTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCCTGGCGAAATGAATTAAACGATGTAACAACAAAgg                                                          >  2:32843/1‑140 (MQ=255)                                                                                    ttGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGACGAAATGAATTAAACGATGTAACAACAAAga                                                          >  2:574858/1‑140 (MQ=255)                                                                                           aaTGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAATAACAAAGTATCACAg                                                   <  1:697838/141‑1 (MQ=255)                                                                                           aaTGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAg                                                   <  1:879328/141‑1 (MQ=255)                                                                                            aggAAAATCCATAATGGAAATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATCTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGt                                                  <  1:411988/139‑1 (MQ=255)                                                                                                aaaTTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTc                                              >  1:830957/1‑141 (MQ=255)                                                                                                 aaTTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCa                                             <  1:223178/141‑1 (MQ=255)                                                                                                 aaTTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCa                                             <  1:939710/141‑1 (MQ=255)                                                                                                      aTAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCatttta                                        >  2:444816/1‑136 (MQ=255)                                                                                                          tgaatttatCCCTAGACGTACTTATCAGACTACTCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGACGACCATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACt                                    <  1:388130/135‑1 (MQ=255)                                                                                                              cATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTaaa                                >  2:535639/1‑141 (MQ=255)                                                                                                                      tAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTAtt                        <  1:311016/141‑1 (MQ=255)                                                                                                                         acttacttATCATACTAATCACATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTATTaaa                     <  1:628033/141‑1 (MQ=255)                                                                                                                         acttacttATCATACTAATCACATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTATTaaa                     <  1:627948/141‑1 (MQ=255)                                                                                                                          cttacttaTCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGACATGTGATGAACGATGTAACAACAAAGTATCACAGTGGACAATATAGACTTAAAATCGAATTAAAt                    >  2:334479/1‑141 (MQ=255)                                                                                                                            tacttaTCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTATTAAATGc                  >  1:499908/1‑141 (MQ=255)                                                                                                                                 aTCATACTAATCGCATAACATGTTCTGCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTATTAAATGCTGaaa             >  1:457573/1‑141 (MQ=255)                                                                                                                                 aTCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTATTAAATGCTGaaa             >  2:793858/1‑141 (MQ=255)                                                                                                                                 aTCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTATTAAATGCTGaaa             >  1:531355/1‑141 (MQ=255)                                                                                                                                     tACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTATTAAATGCTGAAAGCaa         <  1:471888/141‑1 (MQ=255)                                                                                                                                        tAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTATTAAATGCTGAAAGCAAAtt      >  1:848237/1‑141 (MQ=255)                                                                                                                                          aTCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTATTAAACGCTGAAAGCAAATTTg    >  1:595389/1‑141 (MQ=255)                                                                                                                                           tCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTATTAAATGCTGAAAGCAAATTTGa   >  1:91390/1‑141 (MQ=255)                                                                                                                                            cGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGCCTTAGACTCGTATAAAATGCTGAAAGGAAATTTGaa  >  1:357121/1‑141 (MQ=255)                                                                                                                                            cGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTATTAAATGCTGAAAGCAAATTTGaa  >  1:105223/1‑141 (MQ=255)                                                                                                                                             |                                                                                                                                            TTACAATTGATTCTAAGCCAAAAGATGAAGTGAATGAATTTTACGCAAAAACCAAGATTTTGGAGATGTTATAATGCAATTATTTGAAACAATGAAAATTGATAATGGACATATCCCTAGACTTACTTATCATACTAATCGCATAAAATGTTCTTCTGAGCGATTAAACTTTAAATTTGATGAACATGCATGGCGAAATGAATTAAACGATGTAACAACAAAGTATCACAGTGGTCAATATAGACTTAAAATCGTATTAAATGCTGAAAGCAAATTTGAA  >  CP000730/787762‑788041

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