Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	850,742	A→G	59.3%	D98G (GAC→GGC)	ompX →	outer membrane protein X

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	850,742	0	A	G	59.3%	65.8 / 88.5	88	D98G (GAC→GGC)	ompX	outer membrane protein X
Reads supporting (aligned to +/- strand):  ref base A (20/16);  new base G (22/30);  total (42/46)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTT  >  NC_000913/850647‑850837                                                                                                |                                                                                                cACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTgg                                                                                             <  2:1979402/100‑1 (MQ=255) cACTTACACCGAAAAAAGCCGTACTGCAAGCTCTGGAGACTAAAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGAATTAACGACTgg                                                                                             <  2:1695824/100‑1 (MQ=255)   cTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGc                                                                                           >  1:171122/1‑100 (MQ=255)   cTTACACCGAGAAAAGCCGTACTGCAAACTCTGTTCACCACAACAAAAACCAGTACTACGACACCACTGCTGGACCGGCTTACCGCATTAACGACTg                                                                                              >  2:1684794/1‑97 (MQ=255)    ttACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGc                                                                                           <  1:1388529/99‑1 (MQ=255)       cacCGAGAAAAGCCGTACTGCAAGCTATGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGc                                                                                       >  2:2019991/1‑100 (MQ=255)         ccGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGcn                                                                                      <  1:1191824/99‑2 (MQ=255)         ccGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCAt                                                                                     >  1:1708107/1‑100 (MQ=255)         ccGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCa                                                                                      <  1:725090/99‑1 (MQ=255)          cGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATc                                                                                    <  1:619026/100‑1 (MQ=255)           gagaAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCt                                                                                   >  2:431377/1‑100 (MQ=255)             gaAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTAc                                                                                 >  2:37643/1‑100 (MQ=255)               aaaGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACgg                                                                               <  1:406429/100‑1 (MQ=255)                aaGCCGTGCTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGgt                                                                              >  1:1519737/1‑100 (MQ=255)                 aGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGgtg                                                                             <  2:2277387/100‑1 (MQ=255)                 aGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGgtg                                                                             <  1:1638657/100‑1 (MQ=255)                  gCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGgtgt                                                                            <  1:2266167/100‑1 (MQ=255)                     gTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgt                                                                         <  1:1922992/100‑1 (MQ=255)                         tGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtgggt                                                                     >  2:400295/1‑100 (MQ=255)                         tGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggt                                                                     <  2:2179142/100‑1 (MQ=255)                          gCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtgggtg                                                                    <  1:2340162/100‑1 (MQ=255)                          gCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtg                                                                    >  1:1083221/1‑100 (MQ=255)                           cAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtgggtgt                                                                   >  1:1075779/1‑100 (MQ=255)                           cAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgt                                                                   >  2:164598/1‑100 (MQ=255)                              gCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAgtgggtgtggg                                                                <  1:248899/100‑1 (MQ=255)                               ctctGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAgtgggtgtgggt                                                               >  1:2088736/1‑100 (MQ=255)                                  tGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTAt                                                            <  2:1443783/100‑1 (MQ=255)                                   ggTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATg                                                           >  2:2333594/1‑100 (MQ=255)                                    gTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATgg                                                          <  1:471691/100‑1 (MQ=255)                                     tcaCTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGt                                                         <  2:1447432/98‑1 (MQ=255)                                       aCTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTaa                                                       >  1:2092233/1‑100 (MQ=255)                                        cTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTaaa                                                      >  2:807897/1‑100 (MQ=255)                                        cTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTaaa                                                      >  1:695280/1‑100 (MQ=255)                                         tACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAAt                                                     >  2:25724/1‑100 (MQ=255)                                          acaacaAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAAtt                                                    >  1:2066426/1‑100 (MQ=255)                                           caacaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTaaa                                                      >  1:316843/1‑97 (MQ=255)                                             acaaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCa                                                 >  1:2197122/1‑100 (MQ=255)                                              caaAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAg                                                <  1:1387661/100‑1 (MQ=255)                                               aaaaaCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGa                                               >  2:1636336/1‑100 (MQ=255)                                               aaaaaCCAGTACTACGGCATCACTGCTGGTCAGGCTTACCGAATTAACGACTGGGCAAGCATCTACGGTGTATTGGGTGTGGGTTATGGTAAATTCCAg                                                <  2:20632/99‑1 (MQ=255)                                                 aaaCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGAATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGAcc                                             <  2:1974139/100‑1 (MQ=255)                                                  aaCCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCa                                            >  2:1776832/1‑100 (MQ=255)                                                     cAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTg                                         <  2:2308377/100‑1 (MQ=255)                                                     cAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTg                                         <  2:1926576/100‑1 (MQ=255)                                                     cAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTg                                         <  1:1636356/100‑1 (MQ=255)                                                      aGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGa                                        >  2:1248069/1‑100 (MQ=255)                                                         actacGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATa                                     <  1:1441924/100‑1 (MQ=255)                                                         actacGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATa                                     >  2:1863518/1‑100 (MQ=255)                                                             cGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCg                                 <  1:989631/100‑1 (MQ=255)                                                             cGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCg                                 <  2:1079790/100‑1 (MQ=255)                                                             cGGCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCg                                 <  2:23497/100‑1 (MQ=255)                                                               gCATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGAc                               <  1:1370976/100‑1 (MQ=255)                                                               gCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGAc                               <  2:34265/100‑1 (MQ=255)                                                                cATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCg                                 <  1:376397/97‑1 (MQ=255)                                                                cATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGAcc                              <  2:2377418/100‑1 (MQ=255)                                                                cATCACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGAcc                              <  2:1849533/100‑1 (MQ=255)                                                                 aTCACTGCTGGTCCGTCTTACCGCATTAACGGCTGGGCAATCATCTACGGTGTAGTGGGTGTGGGTTATAGTAAATTCCAGACCACTGAATACCCg                                 <  1:379451/96‑1 (MQ=255)                                                                   cACTGCTGGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTAc                           >  1:1769294/1‑100 (MQ=255)                                                                   cACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCAACATa                            >  1:543811/1‑99 (MQ=255)                                                                     ctgctgGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACaa                         >  2:715490/1‑100 (MQ=255)                                                                      tgctgGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACaaa                        >  2:773058/1‑100 (MQ=255)                                                                        ctgGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAaca                      <  1:507344/100‑1 (MQ=255)                                                                        ctgGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAaca                      <  2:703538/100‑1 (MQ=255)                                                                        ccgatctggCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAaca                      <  2:499761/93‑1 (MQ=255)                                                                         tgGTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAacac                     <  1:157396/100‑1 (MQ=255)                                                                          ggTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACg                    >  1:1244255/1‑100 (MQ=255)                                                                          ggTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACg                    <  1:37644/100‑1 (MQ=255)                                                                          ggTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACg                    >  1:2404290/1‑100 (MQ=255)                                                                          ggTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACg                    >  2:1251731/1‑100 (MQ=255)                                                                           gTCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCAACCTAAAAacacaa                   >  1:744992/1‑98 (MQ=255)                                                                           gTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGa                   >  1:2125733/1‑100 (MQ=255)                                                                            tCCGGCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGac                  >  1:1604307/1‑100 (MQ=255)                                                                            tCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGac                  <  2:2062531/100‑1 (MQ=255)                                                                            tCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGac                  <  2:2308677/100‑1 (MQ=255)                                                                             ccGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGaca                 >  1:516081/1‑100 (MQ=255)                                                                             ccGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGaca                 <  2:2404190/100‑1 (MQ=255)                                                                              cGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGacac                >  2:202308/1‑100 (MQ=255)                                                                               ggCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAacaagatcg                >  1:499764/1‑93 (MQ=255)                                                                                gCTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACAccn              <  1:2201126/100‑2 (MQ=255)                                                                                 cTTACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAg             <  1:448598/100‑1 (MQ=255)                                                                                  ttAACGCATTAACGGCTGTGCAAGCATCTACGGTGTAGTGTGTGTGGGTTATGGTAAATTCCAGACTACTGAATACCCGACCTACAAACACGACACCAGc            <  2:2241754/100‑1 (MQ=255)                                                                                   tACCGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCg           <  1:2333622/100‑1 (MQ=255)                                                                                    aCCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGa          >  1:437750/1‑100 (MQ=255)                                                                                     ccGCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCAGACCTACAAACACGACACCAGCGAc         >  2:172914/1‑100 (MQ=255)                                                                                       gCATTAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTa       >  2:1629443/1‑100 (MQ=255)                                                                                          ttAACGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCATTGAATACCCGACCTACAAACACGACACCAGCGACTACgg    <  2:694116/100‑1 (MQ=255)                                                                                           tAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACCGAATACCCGACCTACAAACACGACACCAGCGACTACgg    >  2:1760718/1‑99 (MQ=255)                                                                                            aaCGGCTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGtt  >  1:1772376/1‑100 (MQ=255)                                                                                                |                                                                                                CACTTACACCGAGAAAAGCCGTACTGCAAGCTCTGGTGACTACAACAAAAACCAGTACTACGGCATCACTGCTGGTCCGGCTTACCGCATTAACGACTGGGCAAGCATCTACGGTGTAGTGGGTGTGGGTTATGGTAAATTCCAGACCACTGAATACCCGACCTACAAACACGACACCAGCGACTACGGTT  >  NC_000913/850647‑850837

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: ‑