Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	112,082	T→C	7.0%	D173G (GAT→GGT)	zapD ←	FtsZ stabilizer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	112,082	0	T	C	7.0%	170.5 / 3.4	57	D173G (GAT→GGT)	zapD	FtsZ stabilizer
Reads supporting (aligned to +/- strand):  ref base T (24/29);  new base C (2/2);  total (26/31)
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.97e-01

AAAACGGCTCTTATGTCCGGAAATTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGC  >  NC_000913/111934‑112223                                                                                                                                                     |                                                                                                                                              aaaaCGGCTCTTATGTCCGGAAATTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAAt                                                                                                                                               >  5:268008/1‑149 (MQ=255) aaaaCGGCTCTTATGTCCGGAAATTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAAt                                                                                                                                               <  7:26406/149‑1 (MQ=255) aaaaCGGCTCTTATGTCCGGAAATTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAAt                                                                                                                                               >  3:302898/1‑149 (MQ=255)  aaaCGGCTCTTATGTCCGGAAATTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATc                                                                                                                                              >  6:56478/1‑149 (MQ=255)      ggCTCTTATGTCCGGAAATTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGTGCCGACTGGCGAATTAAATCCAGc                                                                                                                                          <  4:68935/149‑1 (MQ=255)           ttATGTCCGGAAATTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCACTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCAt                                                                                                                                     >  8:236326/1‑149 (MQ=255)             aTGTCCGGAAATTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATg                                                                                                                                    <  2:47939/148‑1 (MQ=255)              tGTCCGGAAATTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCAGCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATgg                                                                                                                                   <  2:252523/148‑1 (MQ=255)                   ggAAATTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTg                                                                                                                             >  2:254351/1‑149 (MQ=255)                   ggAAATTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAACCCAGCACCATGGTAAGTg                                                                                                                             >  2:325149/1‑149 (MQ=255)                   ggAAATTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAACCCAGCACCATGGTAAGTg                                                                                                                             >  2:329800/1‑149 (MQ=255)                       aTTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTg                                                                                                                         <  1:319530/149‑1 (MQ=255)                        tttGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATTGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTgg                                                                                                                        <  5:73740/149‑1 (MQ=255)                        tttGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTgg                                                                                                                        <  2:201956/149‑1 (MQ=255)                        tttGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTgg                                                                                                                        <  4:61818/149‑1 (MQ=255)                        tttGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTgg                                                                                                                        <  3:310113/149‑1 (MQ=255)                           gctgATAAAGCTGGGAATCGATCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCCTTCAGGCGTGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGt                                                                                                                      <  2:60458/145‑1 (MQ=255)                             ggATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAg                                                                                                                   >  4:130560/1‑149 (MQ=255)                                tAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCAACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCgg                                                                                                                <  4:220836/149‑1 (MQ=255)                                   aGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTATGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGtt                                                                                                             >  3:81624/1‑149 (MQ=255)                                       gtgAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAgg                                                                                                         >  3:87898/1‑149 (MQ=255)                                        tgAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAg                                                                                                          <  5:64435/147‑1 (MQ=255)                                           aTCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTgg                                                                                                     <  5:309707/149‑1 (MQ=255)                                            tCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTgg                                                                                                     >  5:96764/1‑148 (MQ=255)                                            tCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGc                                                                                                    <  5:323976/149‑1 (MQ=255)                                                 cGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATcc                                                                                               <  7:153398/149‑1 (MQ=255)                                                 cGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATcc                                                                                               <  4:77414/149‑1 (MQ=255)                                                  gACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCa                                                                                              <  7:168093/149‑1 (MQ=255)                                                  gACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCa                                                                                              <  6:268008/149‑1 (MQ=255)                                                   aCAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAg                                                                                             <  1:254351/149‑1 (MQ=255)                                                      gATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGtt                                                                                          >  3:140077/1‑149 (MQ=255)                                                           aGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTAc                                                                                     <  4:43553/149‑1 (MQ=255)                                                                           cATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTgc                                                                     <  8:272649/149‑1 (MQ=255)                                                                                ccACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTg                                                                >  3:116570/1‑149 (MQ=255)                                                                                 cACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTgg                                                               <  3:29622/149‑1 (MQ=255)                                                                                 cACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTgg                                                               >  6:196949/1‑149 (MQ=255)                                                                                  aCCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTggg                                                              >  5:5174/1‑149 (MQ=255)                                                                                   ccGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTgggg                                                             >  4:25411/1‑149 (MQ=255)                                                                                         tCCTGATAAAAACCACTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATg                                                       <  7:236326/149‑1 (MQ=255)                                                                                           cTGATAAAAACCATTCAGGCTGGGTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCa                                                     >  4:183245/1‑149 (MQ=255)                                                                                              aTAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTCCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGc                                                   >  2:223760/1‑148 (MQ=255)                                                                                                      cATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAAtgtg                                           >  1:247558/1‑148 (MQ=255)                                                                                                            ggCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAACCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTa                                    <  1:325149/149‑1 (MQ=255)                                                                                                            ggCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAACCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTa                                    <  1:329800/149‑1 (MQ=255)                                                                                                             gCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAg                                   >  8:309485/1‑149 (MQ=255)                                                                                                             gCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTGCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAg                                   <  7:171101/149‑1 (MQ=255)                                                                                                               tGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAgg                                  >  7:49610/1‑148 (MQ=255)                                                                                                                           ggAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTg                     >  7:159610/1‑149 (MQ=255)                                                                                                                           ggAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTg                     <  5:56478/149‑1 (MQ=255)                                                                                                                           ggAAGGCGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTCAGCGGCTTCAGGCTGGCAATCCAGGTATCTACCAGCCTGTCGCGCTGCGCCTGGGGTAGATCCAGCCAAAAGTGCAATGTAGGTAAATCAAAGCTg                     >  8:207557/1‑149 (MQ=255)                                                                                                                                 gggCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAgg                                  <  4:90095/130‑1 (MQ=255)                                                                                                                                 gggCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAgg                                  >  3:90095/1‑130 (MQ=255)                                                                                                                                      gACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAgccgcc          >  5:294626/1‑149 (MQ=255)                                                                                                                                       aCTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCt         <  8:49610/149‑1 (MQ=255)                                                                                                                                         tGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTgg       >  8:284731/1‑149 (MQ=255)                                                                                                                                           gCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTggg      <  4:140077/148‑1 (MQ=255)                                                                                                                                           gCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGa     <  4:87898/149‑1 (MQ=255)                                                                                                                                              aaTTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGc  >  4:191768/1‑149 (MQ=255)                                                                                                                                                   aaTCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGc                                                      >  8:139111/1‑92 (MQ=255)                                                                                                                                                   aaTCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGc                                                      <  7:139111/92‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              AAAACGGCTCTTATGTCCGGAAATTTGCGGATAAAGCTGTGAATCGAGCGACAGATTCAGGCGCAGCAAGTCGGCATCGCCACCGTTATCCTGATAAAAACCATTCAGGCTGGTTTGTTTACGGAAGGGGGCCGACTGGCGAATTAAATCCAGCACCATGGTAAGTGCCTGGGTGAGCGGGTTCAGGCTGGCAATCCAGGTTTCTACCTGGCTGTCGCGCTGCGCCTGGGGTAGATGCAGCCAAATGTGCAATGTAGGTAAATCAAAGCTGCAACAGCCGCCTGGGATGC  >  NC_000913/111934‑112223

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