Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,785,565	T→C	*238* (TAA→TGA)	ydiV ←	anti‑FlhDC factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,785,565	0	T	C	95.8%	239.9 / 1.3	72	*238* (TAA→TGA)	ydiV	anti‑FlhDC factor
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (37/32);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

GGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTTATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGCCAATATTTCTGCCGTATCAATGCCTCCCGCAA  >  NZ_CP009273/1785431‑1785688                                                                                                                                       |                                                                                                                            ggCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCg                                                                                                                         <  2:834714/139‑1 (MQ=255)   cGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCt                                                                                                                       >  1:605963/1‑139 (MQ=255)   cGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCt                                                                                                                       <  2:520169/139‑1 (MQ=255)     tAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTg                                                                                                                      >  2:1203480/1‑138 (MQ=255)     tAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTg                                                                                                                      >  2:1160095/1‑138 (MQ=255)     tAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTg                                                                                                                      >  1:1036094/1‑138 (MQ=255)       ggggAAAGTGCTGTATA‑GTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAAcc                                                                                                                  <  1:600242/139‑1 (MQ=255)         ggAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCa                                                                                                                 <  1:1203480/139‑1 (MQ=255)             aGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAAcgt                                                                                                             >  1:1245551/1‑139 (MQ=255)              gTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAAcgtc                                                                                                            <  2:1139815/139‑1 (MQ=255)                         gTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGAtt                                                                                                 >  1:903617/1‑139 (MQ=255)                            tCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTAt                                                                                              >  2:1167489/1‑139 (MQ=255)                              tGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATgg                                                                                            >  2:1089598/1‑139 (MQ=255)                              tGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATgg                                                                                            >  2:358811/1‑139 (MQ=255)                              tGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATgg                                                                                            >  2:206364/1‑139 (MQ=255)                                    ttCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACgg                                                                                      <  1:1160095/139‑1 (MQ=255)                                         cgcgCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGc                                                                                 >  2:407240/1‑139 (MQ=255)                                               aGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTTATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           <  1:577930/139‑1 (MQ=255)                                               aGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           <  2:1238379/139‑1 (MQ=255)                                               aGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           <  1:407240/139‑1 (MQ=255)                                               aGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           <  1:482799/139‑1 (MQ=255)                                               aGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           <  2:1036094/139‑1 (MQ=255)                                               aGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           <  1:939595/139‑1 (MQ=255)                                               aGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           <  1:631218/139‑1 (MQ=255)                                                      aaCAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACgg                                                                                      >  2:278546/1‑121 (MQ=255)                                                      aaCAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACgg                                                                                      <  1:278546/121‑1 (MQ=255)                                                              aCAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCa                                                            >  1:213355/1‑139 (MQ=255)                                                                aGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATg                                                          >  1:144754/1‑139 (MQ=255)                                                                     cATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           <  1:401736/117‑1 (MQ=255)                                                                     cATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           >  2:401736/1‑117 (MQ=255)                                                                            cctcctGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAggg                                              >  2:275583/1‑139 (MQ=255)                                                                             ctcctGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAgggg                                             >  1:956603/1‑139 (MQ=255)                                                                             ctcctGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAgggg                                             >  2:522541/1‑139 (MQ=255)                                                                              tcctGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACgg                                                                                      >  2:927127/1‑97 (MQ=255)                                                                              tcctGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACgg                                                                                      <  1:927127/97‑1 (MQ=255)                                                                              tcctGTTTGATAGTGTAGACGGTTA‑‑CTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGga                                                                 >  1:222566/1‑116 (MQ=255)                                                                              tcctGTTTGATAGTGTAGACGGTTA‑‑CTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGga                                                                 <  2:222566/116‑1 (MQ=255)                                                                                ctGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGt                                           >  1:1048119/1‑138 (MQ=255)                                                                                  gTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCcaca                                                                             >  2:342747/1‑102 (MQ=255)                                                                                  gTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCcaca                                                                             <  1:342747/102‑1 (MQ=255)                                                                                  gTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           >  1:717522/1‑104 (MQ=255)                                                                                  gTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           <  2:717522/104‑1 (MQ=255)                                                                                   tttGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAAt                                        <  1:206364/138‑1 (MQ=255)                                                                                   tttGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATc                                       <  2:956603/139‑1 (MQ=255)                                                                                          gtgtAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTTATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGcc                                 >  2:280672/1‑138 (MQ=255)                                                                                             tAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCAt                                                                        >  2:1192525/1‑96 (MQ=255)                                                                                             tAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCAt                                                                        <  1:1192525/96‑1 (MQ=255)                                                                                                   gTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGCCAATATTTCTg                       <  1:1011126/139‑1 (MQ=255)                                                                                                    ttACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTg                                                                                                                      <  1:537024/43‑1 (MQ=255)                                                                                                    ttACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTg                                                                                                                      >  2:537024/1‑43 (MQ=255)                                                                                                     tACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           <  2:965029/85‑1 (MQ=255)                                                                                                     tACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           >  1:965029/1‑85 (MQ=255)                                                                                                     tACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCAt                                                           <  2:814389/101‑1 (MQ=255)                                                                                                     tACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCAt                                                           >  1:814389/1‑101 (MQ=255)                                                                                                     tACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGa                                                  >  2:33119/1‑110 (MQ=255)                                                                                                     tACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGa                                                  <  1:721637/110‑1 (MQ=255)                                                                                                     tACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGa                                                  >  2:721637/1‑110 (MQ=255)                                                                                                     tACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGa                                                  <  1:33119/110‑1 (MQ=255)                                                                                                     tACTCTCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGCCAATATTTCTGcc                     >  1:299583/1‑139 (MQ=255)                                                                                                       ctctCGTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGCCAATATTTCTGCCg                    >  1:472763/1‑138 (MQ=255)                                                                                                            gTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           <  1:1021639/78‑1 (MQ=255)                                                                                                            gTCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAgg                                                                           >  2:1021639/1‑78 (MQ=255)                                                                                                             tCTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCCGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGCCAATATTTCTGCCGTATCAAt             <  2:1245551/139‑1 (MQ=255)                                                                                                              cTTAAACACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGCCAATATTTCTGCCGTATCAATg            >  1:431116/1‑139 (MQ=255)                                                                                                                 aaaCACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGCCAATATTTCTGCCGTATCAATg            <  2:903617/136‑1 (MQ=255)                                                                                                                 aaaCACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGCCAATATTTCTGCCGTATCAATGcc          <  2:431116/138‑1 (MQ=255)                                                                                                                  aaCACCAGCAAACAGAAGGGTTATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGCCAATATTTCTGCCGTATCAATGCCTc        >  1:240203/1‑139 (MQ=255)                                                                                                                  aaCACCAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGCCAATATTTCTGCCGTATCAATGCCTc        >  1:1231834/1‑139 (MQ=255)                                                                                                                       cAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGa                                                  <  1:362796/92‑1 (MQ=255)                                                                                                                       cAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGa                                                  >  2:362796/1‑92 (MQ=255)                                                                                                                       cAGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGCCAATATTTCTGCCGTATCAATGCCTCCCGc    >  2:825324/1‑138 (MQ=255)                                                                                                                        aGCAAACAGAAGGGTCATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGCCAATATTTCTGCCGTATCAATGCCTCCCGCaa  >  2:891153/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                            GGCGTAGGGGAAAGTGCTGTATAGGTTTCTGGCAATTCATCGCGCCAGCGGGTAACAAAAGACAGGGTCATAGATCCTCCTGTTTGATAGTGTAGACGGTTACTCTCGTCTTAAACACCAGCAAACAGAAGGGTTATCGCTGAACCAACGTCGTTATCTGATTTATGGGAACGGCAGGCCACAGGCATCCCTGGAGAGCATGGAAATCGAAGGGGGTAATCTGTGCCAATATTTCTGCCGTATCAATGCCTCCCGCAA  >  NZ_CP009273/1785431‑1785688

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