Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,019,460	G→A	V359I (GTC→ATC)	fliK →	flagellar hook‑length control protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,019,460	0	G	A	100.0%	229.3 / NA	63	V359I (GTC→ATC)	fliK	flagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (28/35);  total (28/35)

GGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGACACAGGA  >  NC_000913/2019319‑2019605                                                                                                                                               |                                                                                                                                                  ggCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTa                                                                                                                                           <  1:532715/150‑1 (MQ=255)      tCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCCACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAgg                                                                                                                                      <  1:794158/150‑1 (MQ=255)      tCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAgg                                                                                                                                      >  2:1131976/1‑150 (MQ=255)        aGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGc                                                                                                                                    <  4:38728/150‑1 (MQ=255)              ggCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTaa                                                                                                                              <  4:128854/150‑1 (MQ=255)                 aaaGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAg                                                                                                                           >  1:1038344/1‑150 (MQ=255)                        cATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAg                                                                                                                    <  1:662842/150‑1 (MQ=255)                         aTCAGTGGCGAAAGTTTTAGTGGTCCGCAGAAGAACGCTTCCAAGCAATAGTAAAGCCAACGCACAGCAAACCATGAACCTATAGCGGGGGAAGAGGACGATACGCTTTCGGTTA‑CCATCTCTTTAAAAGGGCGTGTAGCAGGCAACAgc                                                                                                                   >  2:914657/1‑150 (MQ=255)                         aTCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAgc                                                                                                                   >  2:577360/1‑150 (MQ=255)                          tCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAgcg                                                                                                                  <  2:260813/150‑1 (MQ=255)                             gTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAgcggcg                                                                                                               >  2:164641/1‑150 (MQ=255)                                gCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTg                                                                                                            >  2:907905/1‑150 (MQ=255)                                  gAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGat                                                                                                          >  1:326001/1‑150 (MQ=255)                                       cTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt                                                                                                     >  2:1129666/1‑150 (MQ=255)                                       cTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATAtttt                                                                                                     >  3:221565/1‑150 (MQ=255)                                        tttAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTc                                                                                                    >  3:292546/1‑150 (MQ=255)                                          tAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGc                                                                                                  <  3:179855/150‑1 (MQ=255)                                          tAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGc                                                                                                  <  1:963389/150‑1 (MQ=255)                                          tAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGc                                                                                                  <  1:929212/150‑1 (MQ=255)                                               gTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                                                                                             >  1:343759/1‑150 (MQ=255)                                               gTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                                                                                             <  1:559346/150‑1 (MQ=255)                                               gTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAAc                                                                                             <  1:187168/150‑1 (MQ=255)                                                 cagcagcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGt                                                                                           >  4:353615/1‑150 (MQ=255)                                                      gcagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCagag                                                                                      <  1:1101286/150‑1 (MQ=255)                                                       cagGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAgg                                                                                     <  1:929338/150‑1 (MQ=255)                                                           ccGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGc                                                                                 <  2:57364/150‑1 (MQ=255)                                                            cGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCa                                                                                >  2:303423/1‑150 (MQ=255)                                                             gCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCAc                                                                               >  1:711479/1‑150 (MQ=255)                                                              cTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCAcc                                                                              <  3:164505/150‑1 (MQ=255)                                                               ttCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCg                                                                             >  1:79331/1‑150 (MQ=255)                                                                   cagcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCttt                                                                                                                                            <  1:410168/83‑1 (MQ=255)                                                                   cagcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCttt                                                                                                                                            >  2:410168/1‑83 (MQ=255)                                                                     gcaacagcaaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATcc                                                                       >  3:340216/1‑150 (MQ=255)                                                                            caaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACCGGCAACAGCGGCGTTGATTGTTGCGCCTAACGTCTGAGGGGGCACCGTAAGCCGCGTCtt                                                                >  1:398597/1‑150 (MQ=255)                                                                            caaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                                                                                             <  2:158434/105‑1 (MQ=255)                                                                            caaAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGtt                                                                                                             >  1:158434/1‑105 (MQ=255)                                                                                 ccAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTccc                                                           <  4:340216/150‑1 (MQ=255)                                                                                   aaCGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCg                                                         >  2:191121/1‑150 (MQ=255)                                                                                   aaCGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCg                                                         >  2:191037/1‑150 (MQ=255)                                                                                      gCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCtt                                                      >  1:864301/1‑150 (MQ=255)                                                                                        acaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATcc                                                                       <  4:216813/131‑1 (MQ=255)                                                                                        acaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATcc                                                                       >  3:216813/1‑131 (MQ=255)                                                                                        acaGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTttg                                                    <  1:1131976/150‑1 (MQ=255)                                                                                          aGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTttgtt                                                  <  4:221565/150‑1 (MQ=255)                                                                                          aGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTttgtt                                                  <  2:819171/150‑1 (MQ=255)                                                                                                  aTGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCaa                                          <  2:63672/150‑1 (MQ=255)                                                                                                  aTGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCaa                                          <  2:28433/150‑1 (MQ=255)                                                                                                    gAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGa                                        <  1:602857/150‑1 (MQ=255)                                                                                                       ccTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGc                                     <  1:1076675/150‑1 (MQ=255)                                                                                                         tctGGCGGGGCAAGACGACGATAC‑AGTCGTGTGCTCCACCGCTCTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAg                                   <  1:905948/150‑1 (MQ=255)                                                                                                            ggCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACNGTANTCCGCGTCTTTTCCNCGCTTTGTTGCGCTCAAGACGCAGGAt                                <  3:413962/150‑1 (MQ=255)                                                                                                               gggggAAGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAAt                             <  1:459583/150‑1 (MQ=255)                                                                                                                    aaGACGACGATACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGcc                        <  2:326001/150‑1 (MQ=255)                                                                                                                        cgacgaTACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATa                    <  2:680922/150‑1 (MQ=255)                                                                                                                         gacgaTACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATaa                   <  1:164641/150‑1 (MQ=255)                                                                                                                              tACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGt              >  3:208433/1‑150 (MQ=255)                                                                                                                              tACGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGt              >  3:234348/1‑150 (MQ=255)                                                                                                                                cGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAg            <  2:386280/150‑1 (MQ=255)                                                                                                                                cGCTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAg            <  1:577360/150‑1 (MQ=255)                                                                                                                                  cTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCg          >  3:146150/1‑150 (MQ=255)                                                                                                                                  cTTCCGGTTC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCg          >  3:353483/1‑150 (MQ=255)                                                                                                                                         ttC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGACACAgg   >  2:76454/1‑150 (MQ=255)                                                                                                                                          tC‑CCATCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGTAGTAGCGACACAGGa  <  2:711479/150‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                  GGCATTCAGTTAGGGCAAAGCAACATCAGTGGCGAAAGCTTTAGTGGTCAGCAGCAGGCCGCTTCCCAGCAACAGCAAAGCCAACGCACAGCAAACCATGAACCTCTGGCGGGGGAAGACGACGATACGCTTCCGGTTC‑CCGTCTCTTTACAAGGGCGTGTAACAGGCAACAGCGGCGTTGATATTTTCGCCTAACGTCAGAGGTAGCACCGTAATCCGCGTCTTTTCCCCGCTTTGTTGCGCTCAAGACGCAGGATAATTAGCCGATAAGCAGTAGCGACACAGGA  >  NC_000913/2019319‑2019605

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