Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,892,935	+T	coding (301/1911 nt)	yoaA ←	putative ATP‑dependent helicase, DinG family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,892,932	1	.	T	94.6%	247.4 / 6.1	74	coding (304/1911 nt)	yoaA	putative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base . (1/3);  new base T (35/35);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCC‑T‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCAACGCTTTTGAGC  >  NC_000913/1892833‑1893024                                                                                                     |                                                                                               aGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGccc                                                                                                 <  1:1828328/100‑1 (MQ=255)  gcgcAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCC‑t                                                                                               >  2:2254180/1‑100 (MQ=255)   cgccgcAGGAACACAGCACTTAAGAATTGTACCGGCAGATCGCATCCCGCCAGCGCCTGCTGGTCGAAACGTTCGAGGCAGAGGTCGTTTGAGCGCCC‑t‑t                                                                                             >  2:1968739/1‑100 (MQ=255)      agcagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCC‑t‑tttc                                                                                          <  2:1202567/100‑3 (MQ=255)       gcagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCC‑t‑tttca                                                                                         >  1:1623638/1‑97 (MQ=255)       gcagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCC‑t‑tttca                                                                                         <  1:2070637/100‑4 (MQ=255)       gcagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCC‑t‑tttca                                                                                         >  1:2236037/1‑97 (MQ=255)       gcagGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCC‑t‑tttca                                                                                         <  1:2615750/100‑4 (MQ=255)         agGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCC‑t‑tttcagc                                                                                       <  1:1201848/100‑6 (MQ=255)          ggATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑cagca                                                                                      >  1:1648772/1‑100 (MQ=255)          ggATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑cagca                                                                                      >  1:1096207/1‑100 (MQ=255)          ggATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑cagca                                                                                      <  2:668319/100‑1 (MQ=255)             tCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAgcg                                                                                   >  2:692308/1‑100 (MQ=255)             tCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAgcg                                                                                   <  1:883864/100‑1 (MQ=255)                caTCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCa                                                                                <  2:513155/100‑1 (MQ=255)                caTCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCa                                                                                >  1:1028243/1‑100 (MQ=255)                    gCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGtt                                                                            >  1:1326499/1‑100 (MQ=255)                    gCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCC‑T‑TT‑CAGCAGCGCCACGTTg                                                                           >  2:2400146/1‑100 (MQ=255)                      ttAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGc                                                                          <  1:67577/100‑1 (MQ=255)                      ttAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGc                                                                          <  1:67606/100‑1 (MQ=255)                       tAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGcc                                                                         <  2:1343402/100‑1 (MQ=255)                        aaGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGccc                                                                        <  1:197210/100‑1 (MQ=255)                        aaGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCGGCTGTGCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGccc                                                                        <  2:127293/100‑1 (MQ=255)                          gATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGt                                                                      >  2:1894599/1‑100 (MQ=255)                          gATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGt                                                                      >  1:1013583/1‑100 (MQ=255)                           aTTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGta                                                                     <  1:772469/100‑1 (MQ=255)                            tttGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGtat                                                                    >  1:1220286/1‑100 (MQ=255)                              tGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTatat                                                                  >  1:1271373/1‑100 (MQ=255)                               gTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATAtt                                                                 <  1:1591076/100‑1 (MQ=255)                                  ccGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCa                                                              >  1:2514277/1‑100 (MQ=255)                                  ccGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCa                                                              >  2:58165/1‑100 (MQ=255)                                  ccGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGATTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCa                                                              >  2:1556604/1‑100 (MQ=255)                                      cAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGc                                                          >  1:2317252/1‑100 (MQ=255)                                       aGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGcc                                                         <  2:1859144/100‑1 (MQ=255)                                        gATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCt                                                        >  1:1539683/1‑100 (MQ=255)                                           cGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTg                                                     <  2:2331707/100‑1 (MQ=255)                                            gCCCCCCGCCATCGCCGTCTGTTCGATACGTTCTGTGCATAGGTAGTTTGACCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTACATTTCAATGCCTTTga                                                    >  2:2085212/1‑100 (MQ=255)                                            gCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTga                                                    <  2:2066833/100‑1 (MQ=255)                                            gCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTga                                                    >  2:149139/1‑100 (MQ=255)                                             ccccccGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTgag                                                   <  2:840208/100‑1 (MQ=255)                                             ccccccGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTgag                                                   <  1:2454063/100‑1 (MQ=255)                                                    ccAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTg                                            <  2:307718/100‑1 (MQ=255)                                                    ccAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTg                                            <  1:1629170/100‑1 (MQ=255)                                                      aGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGc                                          >  2:798475/1‑100 (MQ=255)                                                      aGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCC‑TCTT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGc                                          <  2:2306039/100‑1 (MQ=255)                                                       gcgcCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCa                                         <  1:2435946/100‑1 (MQ=255)                                                        cgcCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCaa                                        <  2:973597/100‑1 (MQ=255)                                                          ccTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAAt                                      >  2:2060656/1‑100 (MQ=255)                                                             gctgTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATcgc                                   >  1:2140631/1‑100 (MQ=255)                                                              ctgTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATcgcg                                  >  1:1600872/1‑100 (MQ=255)                                                                gTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCC‑T‑TTGCAGCAGCGCCACGTTGCCCGGATAGTTCAAGGCCTTTGCGACTGTTGGCAAATCGCGGc                                <  2:2337442/100‑1 (MQ=255)                                                                   cGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGt                             >  1:1086096/1‑100 (MQ=255)                                                                     agaCGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTag                           <  1:2268319/100‑1 (MQ=255)                                                                      gaCGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTaga                          >  2:816898/1‑100 (MQ=255)                                                                        cGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGc                        >  2:2239986/1‑100 (MQ=255)                                                                        cGTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGc                        >  1:1504057/1‑100 (MQ=255)                                                                         gTTCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCt                       >  2:528378/1‑100 (MQ=255)                                                                          ttCGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTg                      <  2:1984009/100‑1 (MQ=255)                                                                            cGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGAt                    <  2:256874/100‑1 (MQ=255)                                                                            cGAGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGAt                    >  2:769233/1‑100 (MQ=255)                                                                              aGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAATTCGCGGCTGTAGCGCTGATcc                  >  2:1634764/1‑100 (MQ=255)                                                                              aGGCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATcc                  <  1:1713836/100‑1 (MQ=255)                                                                               ggCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCt                 >  1:1001492/1‑100 (MQ=255)                                                                               ggCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCt                 >  1:1199485/1‑100 (MQ=255)                                                                               ggCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCt                 <  1:1015393/100‑1 (MQ=255)                                                                               ggCAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCt                 >  1:434877/1‑100 (MQ=255)                                                                                 cAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGc               <  2:1252927/100‑1 (MQ=255)                                                                                 cAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGc               <  2:948869/100‑1 (MQ=255)                                                                                 cAGAGGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGc               >  2:1007837/1‑100 (MQ=255)                                                                                  agagGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCa              <  1:1078207/100‑1 (MQ=255)                                                                                  agagGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCa              <  1:1221639/100‑1 (MQ=255)                                                                                  agagGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCa              <  1:1221651/100‑1 (MQ=255)                                                                                   gagGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCaa             >  1:1003019/1‑100 (MQ=255)                                                                                    agGTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCAAc            >  1:1582972/1‑100 (MQ=255)                                                                                      gTAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCAACGc          >  2:913805/1‑100 (MQ=255)                                                                                       tAGTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCAACGCt         <  2:1626208/100‑1 (MQ=255)                                                                                         gTTTGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCAACGCttt       <  1:2625458/100‑1 (MQ=255)                                                                                          tttGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCAACGCtttt      >  1:517122/1‑100 (MQ=255)                                                                                          tttGAGCGCCC‑T‑TTGCAGCAGCGCCACGTGGCCCGTATATTTCAAAGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCAACACtttt      <  2:137652/100‑1 (MQ=255)                                                                                           ttGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCAACGCTTTTg     <  2:1202285/100‑1 (MQ=255)                                                                                            tGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCAACGCTTTTGa    <  2:472611/100‑1 (MQ=255)                                                                                            tGAGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCAACGCTTTTGa    <  1:1691232/100‑1 (MQ=255)                                                                                              aGCGCCCTT‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCAACGCTTTTGAGc  <  2:790397/100‑1 (MQ=255)                                                                                                     |                                                                                               AGCGCAGCAGGATCACATCGCTTAAGATTTGTACCGGCAGATCGCCCCCCGCCAGCGCCTGCTGTTCGAGACGTTCGAGGCAGAGGTAGTTTGAGCGCCC‑T‑TT‑CAGCAGCGCCACGTTGCCCGTATATTTCAATGCCTTTGAGACTGTTGGCAAATCGCGGCTGTAGAGCTGATCCTGCAACGCTTTTGAGC  >  NC_000913/1892833‑1893024

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