Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	107,553	(T)9→8	100%	intergenic (+79/‑152)	lpxC → / → secM	UDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	107,545	0	T	.	91.7%	269.1 / 19.8	72	intergenic (+71/‑160)	lpxC/secM	UDP‑3‑O‑acyl N‑acetylglucosamine deacetylase/regulator of secA translation
Reads supporting (aligned to +/- strand):  ref base T (5/1);  new base . (37/29);  total (42/30)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

TCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGA  >  NC_000913/107446‑107639                                                                                                    |                                                                                               tCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGt                                                                                                >  2:2032629/1‑100 (MQ=255)   aaaGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttt                                                                                              >  2:956988/1‑100 (MQ=255)     aGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttt                                                                                            >  2:2150292/1‑100 (MQ=255)     aGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttt                                                                                            >  2:883552/1‑100 (MQ=255)      gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttt                                                                                           >  2:1782400/1‑100 (MQ=255)      gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttt                                                                                           >  2:2177996/1‑100 (MQ=255)      gcgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttt                                                                                           >  2:126516/1‑100 (MQ=255)       cgcCTTCAGCTGTACTTGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGAAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                          <  2:1281587/100‑1 (MQ=255)       cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                          <  2:1455968/100‑1 (MQ=255)       cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                          >  1:570272/1‑100 (MQ=255)       cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                          <  1:2578403/100‑1 (MQ=255)       cgcCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttt                                                                                          <  1:2507943/100‑1 (MQ=255)        nccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGtttttttt                                                                                         >  1:2184495/2‑100 (MQ=255)         ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttt                                                                                        <  2:1009158/100‑1 (MQ=255)         ccTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                        <  1:652935/100‑2 (MQ=255)         ccTTCAGCTGCACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACACTCTCCGGCCAGGTGAACCAGTCGttttttttg                                                                                        <  1:1912140/100‑2 (MQ=255)          cTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGttttttttga                                                                                       >  2:1040916/1‑98 (MQ=255)             cAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAAtt                                                                                   >  2:124315/1‑100 (MQ=255)               gCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAAtttt                                                                                 >  1:505357/1‑100 (MQ=255)                cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTta                                                                                >  1:746339/1‑100 (MQ=255)                cTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTta                                                                                <  1:2748428/100‑1 (MQ=255)                  gTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTtata                                                                              <  1:2813458/100‑1 (MQ=255)                    aCTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAag                                                                            <  2:2651336/100‑1 (MQ=255)                       ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAaga                                                                           >  1:2248213/1‑98 (MQ=255)                       ggCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGc                                                                         >  2:63403/1‑100 (MQ=255)                       gacATAACGACATGTATACTGTCGTATAAAATTCGACTGGCAGATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGc                                                                         >  2:1968575/3‑100 (MQ=255)                        gCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGTTTTTTATAagagtt                                                                        >  1:2135559/1‑98 (MQ=255)                         cATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCta                                                                       >  2:940849/1‑100 (MQ=255)                         cATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCta                                                                       >  1:573691/1‑100 (MQ=255)                          aTAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCtat                                                                      >  2:2205289/1‑100 (MQ=255)                             aCGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTATTTTGGATTTTATAAGAGCTATaaa                                                                   >  1:150185/1‑100 (MQ=255)                                 cagttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                               <  1:171321/97‑1 (MQ=255)                                 cATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAAc                                                               >  2:365868/1‑100 (MQ=255)                                   tttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACgg                                                             >  2:485318/1‑100 (MQ=255)                                    ttATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGt                                                            >  2:2605178/1‑100 (MQ=255)                                     tataCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTg                                                           <  2:2823013/100‑1 (MQ=255)                                       taCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCg                                                         <  1:2547394/100‑1 (MQ=255)                                       taCTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCg                                                         <  2:2501122/100‑1 (MQ=255)                                         cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                       >  1:2217754/1‑100 (MQ=255)                                         cTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGaa                                                       >  1:1231514/1‑100 (MQ=255)                                           gTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACg                                                     >  2:1706319/1‑100 (MQ=255)                                            tCGTATAAAATTCGAGTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGc                                                    <  2:2859658/100‑1 (MQ=255)                                              gTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTT‑TTTTTTTTAAGAGCTATAAAAAACGGTGCGAACGCTg                                                  >  1:1060496/1‑100 (MQ=255)                                               tataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTT‑TTTTTTTTAAGAGCTATAAAAAACGGTGCGAACGCTGt                                                 >  1:538262/1‑100 (MQ=255)                                               tataAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCtat                                                 >  2:1502820/1‑98 (MQ=255)                                                 taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                               <  2:1155475/100‑1 (MQ=255)                                                 taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                               >  1:2459043/1‑100 (MQ=255)                                                 taAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGttt                                               <  2:2459054/100‑1 (MQ=255)                                                   aaaTTCGACTGGCAAATCTGGCACTCTCTCCGGCAAGGTGAAATAGTCG‑GTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGAGCGAACGCTGTTTTc                                             >  1:2591443/1‑100 (MQ=255)                                                    aaTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCt                                            <  1:2741640/100‑1 (MQ=255)                                                     aTTCGACTGGGAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                           <  1:709796/100‑1 (MQ=255)                                                     aTTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCtt                                           <  1:2266376/100‑1 (MQ=255)                                                      ttCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTa                                          <  1:2870822/100‑1 (MQ=255)                                                        cGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAg                                        <  2:2889559/100‑1 (MQ=255)                                                         gACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGc                                       <  1:1103781/100‑1 (MQ=255)                                                          aCTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCa                                      >  1:2364134/1‑100 (MQ=255)                                                           cTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                     >  2:1877087/1‑100 (MQ=255)                                                             ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCAc                                     >  1:1158949/1‑98 (MQ=255)                                                             ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                   >  2:1711082/1‑100 (MQ=255)                                                             ggCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACtt                                   >  2:927494/1‑100 (MQ=255)                                                              gCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACttt                                  >  2:2571940/1‑100 (MQ=255)                                                                  aTCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                              <  1:602811/100‑1 (MQ=255)                                                                  aTCTGGCACTCTCTCCGACCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCg                              <  1:681764/100‑1 (MQ=255)                                                                   tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCATAAGCACTTTTCCGc                             >  2:2868927/1‑100 (MQ=255)                                                                   tCTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTCTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                             <  2:2042094/100‑1 (MQ=255)                                                                    cTGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGca                            <  2:1682551/100‑1 (MQ=255)                                                                    cTGGCACTCTCTCCGGACAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGc                             <  1:2175348/99‑1 (MQ=255)                                                                     tGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                           >  2:722628/1‑100 (MQ=255)                                                                     tGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                           >  1:1769844/1‑100 (MQ=255)                                                                     tGGCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                           <  1:1998251/100‑1 (MQ=255)                                                                      ggCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGcac                           <  1:1592443/100‑1 (MQ=255)                                                                       gCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACaa                         >  2:1759104/1‑100 (MQ=255)                                                                       gCACTCTCTCCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACaa                         >  1:395588/1‑100 (MQ=255)                                                                           tctctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTa                     >  1:203728/1‑100 (MQ=255)                                                                            ctctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTAt                    >  2:1755831/1‑100 (MQ=255)                                                                             tctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATc                   >  2:630675/1‑100 (MQ=255)                                                                              ctcCGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCt                  <  2:1017878/100‑1 (MQ=255)                                                                              ccccGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTa                     >  2:2081761/3‑97 (MQ=255)                                                                                ccGGCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTc                >  2:2410682/1‑100 (MQ=255)                                                                                  ggCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCAt              >  1:1902737/1‑100 (MQ=255)                                                                                  ggCCAGGTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCAt              <  2:1442973/100‑1 (MQ=255)                                                                                     cAGGTGAACCAGTCGTTTTTTTTT‑TTTTTTTTAAAACAGAAAAACGACGGCGAAAACGCGATGTGCTTAAGAACTTTTCCGCACAACTTactatcat              >  1:1229707/1‑90 (MQ=255)                                                                                       ggTGAACCAGTCGTTTTTTTTT‑TTTTTTTTAAGAGCAAAAAAACACGGTTATACCGCTGTTTTTTTAAGCCCTTTTCCTGCCCACTTATCTTCATTCg           >  1:501271/1‑98 (MQ=255)                                                                                       ggTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTg         <  2:989090/100‑1 (MQ=255)                                                                                        gTGAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGc        <  2:963740/100‑1 (MQ=255)                                                                                          gAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtg      <  2:407325/100‑1 (MQ=255)                                                                                          gAACCAGTCG‑TTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCtg      <  2:682144/100‑1 (MQ=255)                                                                                           aaCCAGTCG‑TTTTTTTTGAATTTTATAAAAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTCCtg      <  1:1341480/99‑1 (MQ=255)                                                                                                gTCGTTTTTTTTT‑AATTTTATAAGAGCTATAAAAAACGGTTCGAACGCTTTTTTCGTTAGCACTTTTCCGCACAACTTATCTTCATTCGTTCTGTGGa  >  1:577485/1‑98 (MQ=255)                                                                                                    |                                                                                               TCAAAGCGCCTTCAGCTGTACTGGCATAACGACATTTATACTGTCGTATAAAATTCGACTGGCAAATCTGGCACTCTCTCCGGCCAGGTGAACCAGTCGTTTTTTTTTGAATTTTATAAGAGCTATAAAAAACGGTGCGAACGCTGTTTTCTTAAGCACTTTTCCGCACAACTTATCTTCATTCGTGCTGTGGA  >  NC_000913/107446‑107639

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