Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,678,937	T→A	Q84L (CAG→CTG)	yhjH ←	cyclic‑di‑GMP phosphodiesterase, FlhDC‑regulated

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,678,937	0	T	A	100.0%	311.0 / NA	87	Q84L (CAG→CTG)	yhjH	cyclic‑di‑GMP phosphodiesterase, FlhDC‑regulated
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (50/37);  total (50/37)

CTGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGA  >  NC_000913/3678840‑3679030                                                                                                  |                                                                                              cTGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGc                                                                                             <  2:671559/100‑1 (MQ=255)  tGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGct                                                                                            >  1:1116016/1‑100 (MQ=255)     cGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCtt                                                                                         >  2:1920823/1‑100 (MQ=255)     cGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCtt                                                                                         >  1:1656326/1‑100 (MQ=255)      gCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGGAGGCCGTGGTCTATAGAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCttt                                                                                        >  2:1951235/1‑100 (MQ=255)       naGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTc                                                                                       >  2:2319316/2‑100 (MQ=255)       cAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTc                                                                                       <  1:1979619/100‑1 (MQ=255)           gCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACaa                                                                                   <  1:1358844/100‑1 (MQ=255)           gCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACaa                                                                                   <  1:1699088/100‑1 (MQ=255)           gCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACaa                                                                                   >  1:1186115/1‑100 (MQ=255)              aTGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCt                                                                                >  1:411042/1‑100 (MQ=255)                gAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTcc                                                                              >  2:1388078/1‑100 (MQ=255)                   cGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATc                                                                           >  1:618529/1‑100 (MQ=255)                       gggCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCAGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGAt                                                                       <  1:1405797/100‑1 (MQ=255)                       atctcATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGAt                                                                       <  2:301010/96‑1 (MQ=255)                         gCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATgg                                                                     >  1:576694/1‑100 (MQ=255)                         gCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATgg                                                                     >  2:1704060/1‑100 (MQ=255)                         gCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATgg                                                                     >  1:2151711/1‑100 (MQ=255)                           cATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGAtagat                                                                   >  1:301012/1‑96 (MQ=255)                           cATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCt                                                                   >  1:2310845/1‑100 (MQ=255)                             tCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGa                                                                 <  2:1416080/100‑1 (MQ=255)                             tCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGa                                                                 >  2:1404048/1‑100 (MQ=255)                             tCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGa                                                                 >  2:1871551/1‑100 (MQ=255)                              cAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGAc                                                                >  2:568330/1‑100 (MQ=255)                                  atTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTg                                                            >  2:522392/1‑100 (MQ=255)                                  atTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTg                                                            <  2:576693/100‑1 (MQ=255)                                   ttGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGa                                                           <  2:618528/100‑1 (MQ=255)                                     gACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGAtt                                                         >  1:2076129/1‑100 (MQ=255)                                     gACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGAtt                                                         >  2:914292/1‑100 (MQ=255)                                      aCCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGAttt                                                        <  1:277360/100‑1 (MQ=255)                                       ctgATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTc                                                       <  1:1920207/98‑1 (MQ=255)                                       ccGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTc                                                       <  2:1677566/100‑1 (MQ=255)                                       ccGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTc                                                       <  1:2308376/100‑1 (MQ=255)                                         gATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAg                                                     >  2:1920200/1‑100 (MQ=255)                                            gCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTaa                                                  >  2:1371735/1‑100 (MQ=255)                                             ccAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTaaa                                                 >  2:2090084/1‑100 (MQ=255)                                             ccAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTaaa                                                 >  1:200942/1‑100 (MQ=255)                                             ccAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTaaa                                                 <  2:955729/100‑1 (MQ=255)                                             ccAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATAAGCTCTTTCACAACCTCCATCCGATGGCTGACGGGGATTTCAGTaaa                                                 >  2:1158372/1‑100 (MQ=255)                                                gcagGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATa                                              >  1:140118/1‑100 (MQ=255)                                                 tccgatctGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAg                                             <  1:422353/93‑1 (MQ=255)                                                  agGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGc                                            <  1:522394/100‑1 (MQ=255)                                                  agGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGc                                            <  1:149233/100‑1 (MQ=255)                                                   ggCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCg                                           >  1:969661/1‑100 (MQ=255)                                                    gCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGa                                          >  1:979258/1‑100 (MQ=255)                                                     ccGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGAt                                         <  2:276311/100‑1 (MQ=255)                                                      cGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATc                                        <  1:426673/100‑1 (MQ=255)                                                         gCTCTATAAAGAAGTCGGCTTTTTTTTTTTTCATATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGAGATcgga                                     >  2:422351/1‑97 (MQ=255)                                                         gCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCcgg                                     >  2:1102479/1‑100 (MQ=255)                                                           tctATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCcggcg                                   <  2:1856099/100‑1 (MQ=255)                                                            ctATAAAGAAGTCGGCTTTTTGTTCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGAAGGTGATTTCAGTAAAATAGCGATCcggcgg                                  >  2:1638325/1‑100 (MQ=255)                                                            ctATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGCTTTCAGTAAAATAGCGATCcggcgg                                  >  2:2209796/1‑100 (MQ=255)                                                            ctATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCcggcgg                                  <  2:2310836/100‑1 (MQ=255)                                                            ctATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCcgg                                     >  2:340177/1‑97 (MQ=255)                                                             tataAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCggcggc                                 >  1:2115425/1‑100 (MQ=255)                                                             tataAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCggcggc                                 >  2:854915/1‑100 (MQ=255)                                                             tataAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGGGATTTCAGTAAAATATCGATCCggcggc                                 >  2:2182995/1‑100 (MQ=255)                                                              ataAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCa                                >  1:52147/1‑100 (MQ=255)                                                               taAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTATTTCACAACCTCCATCCGATGGCTGAAGGTGATTTCAGTAAAATAGCGATCCGGCGGCAg                               >  2:1906690/1‑100 (MQ=255)                                                                aaaGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAgg                              >  2:1855043/1‑100 (MQ=255)                                                                     agTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTg                         <  1:2288646/100‑1 (MQ=255)                                                                     agTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTg                         <  1:2319325/100‑1 (MQ=255)                                                                      gTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGc                        >  1:1492065/1‑100 (MQ=255)                                                                        cGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGa                      >  2:1215396/1‑100 (MQ=255)                                                                        cGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGa                      <  2:2051303/100‑1 (MQ=255)                                                                        cGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGa                      <  2:2051284/100‑1 (MQ=255)                                                                        cGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGa                      >  1:820159/1‑100 (MQ=255)                                                                          gCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAg                    <  1:2090091/100‑1 (MQ=255)                                                                          gCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAg                    >  1:725600/1‑100 (MQ=255)                                                                            tttttGCTCCAGCAAATCAATCATCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAggg                  <  1:2183002/100‑1 (MQ=255)                                                                            tttttGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAggg                  <  2:2151704/100‑1 (MQ=255)                                                                               ttGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGt                 <  1:1228692/98‑1 (MQ=255)                                                                               ttGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTc               <  2:494311/100‑1 (MQ=255)                                                                                tGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCa              <  1:1413909/100‑1 (MQ=255)                                                                                tGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCa              <  2:411040/100‑1 (MQ=255)                                                                                  cgcCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAg            >  1:855953/1‑100 (MQ=255)                                                                                   gcCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAgg           >  2:2116882/1‑100 (MQ=255)                                                                                   gcCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAgg           <  2:1116014/100‑1 (MQ=255)                                                                                   gcCAGCAAATCAATCAGCTCTTTCACAACCTCCAGCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAgg           <  2:1188395/100‑1 (MQ=255)                                                                                    ccAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAggg          >  1:548460/1‑100 (MQ=255)                                                                                     cAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGa         >  1:2033934/1‑100 (MQ=255)                                                                                       gCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATg       <  1:1404051/100‑1 (MQ=255)                                                                                        cAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAGGGGATGc      >  1:1056193/1‑100 (MQ=255)                                                                                        cAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGc      >  1:1137035/1‑100 (MQ=255)                                                                                         aaaTCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCg     <  2:969659/100‑1 (MQ=255)                                                                                          aatcaatcAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGt    >  1:757594/1‑100 (MQ=255)                                                                                           atcaatcaGCTCTTTCACAACCTCCATTCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTg   <  2:184900/100‑1 (MQ=255)                                                                                            tcaatcAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGa  <  1:360680/100‑1 (MQ=255)                                                                                                  |                                                                                              CTGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGA  >  NC_000913/3678840‑3679030

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