Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	652,430	T→G	S66A (TCC→GCC)	citA →	sensory histidine kinase in two‑component regulatory system with CitB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	652,430	0	T	G	100.0%	235.0 / NA	63	S66A (TCC→GCC)	citA	sensory histidine kinase in two‑component regulatory system with CitB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (30/33);  total (30/33)

GGCATTTCCGCTGCGTATCTTTTTGCTGATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGATCCGCCTTTACCATCCTAATCCGG  >  NC_000913/652282‑652574                                                                                                                                                     |                                                                                                                                                 ggCATTTCCGCTGCGTATCTTTTTGCTGATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTgccgcc                                                                                                                                                >  1:462170/1‑151 (MQ=255) ggCATTTCCGCTGCGTATCTTTTTGCTGATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTgccgcc                                                                                                                                                >  2:570029/1‑151 (MQ=255)  gCATTTCCGCTGCGTATCTTTTTGCTGATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCa                                                                                                                                               <  1:308028/151‑1 (MQ=255)        ccGCTGCGTATCTTTTTGCTGATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGAc                                                                                                                                          <  1:710167/150‑1 (MQ=255)        ccGCTGCGTATCTTTTTGCTGATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGAc                                                                                                                                          >  1:351278/1‑150 (MQ=255)         cGCTGCGTATCTTTTTGCTGATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAg                                                                                                                                        <  2:72437/151‑1 (MQ=255)            tGCGTATCTTTTTGCTGATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAgt                                                                                                                                       <  2:856269/149‑1 (MQ=255)             gCGTATCTTTTTGCTGATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAgtgt                                                                                                                                     <  2:233123/150‑1 (MQ=255)             gCGTATCTTTTTGCTGATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAgtg                                                                                                                                      <  2:43094/149‑1 (MQ=255)             gCGTATCTTTTTGCTGATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTc                                                                                                                                    >  1:731080/1‑151 (MQ=255)              cGTATCTTTTTGCTGATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCa                                                                                                                                   <  1:80765/151‑1 (MQ=255)                            gATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGaa                                                                                                                      <  1:710688/150‑1 (MQ=255)                             aTTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGaaa                                                                                                                     >  1:734836/1‑150 (MQ=255)                             aTTCTGGTGTTCTCAATAGTTGACATTGCAGCCCTAGAGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAAGCAGGCCAAAATTATTGCCGCCAATGACATTGTCATCTCTGCGGTGa                                                                                                                       <  2:734240/148‑1 (MQ=255)                               tCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAAc                                                                                                                   <  2:759194/150‑1 (MQ=255)                                  ggTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAAcgcg                                                                                                                >  1:850479/1‑150 (MQ=255)                                       tctcAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACt                                                                                                           >  1:265025/1‑150 (MQ=255)                                          cAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACa                                                                                                        <  1:496688/150‑1 (MQ=255)                                            atatTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAAc                                                                                                     >  2:438563/1‑151 (MQ=255)                                             tatTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACg                                                                                                    >  1:563923/1‑151 (MQ=255)                                              atTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACgg                                                                                                   >  2:42957/1‑151 (MQ=255)                                              atTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACg                                                                                                    <  2:634615/150‑1 (MQ=255)                                              atTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACg                                                                                                    <  2:185917/150‑1 (MQ=255)                                               tttGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGc                                                                                                  >  2:31227/1‑151 (MQ=255)                                                 tGTCATTGCAGCCCTGGCACAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTg                                                                                                >  2:197006/1‑151 (MQ=255)                                                  gTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTg                                                                                                <  2:396768/150‑1 (MQ=255)                                                              cTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCt                                                                                    >  2:785490/1‑150 (MQ=255)                                                                  cgcAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACaa                                                                               <  1:344067/151‑1 (MQ=255)                                                                     aaTATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAAt                                                                             <  1:598393/150‑1 (MQ=255)                                                                      atatTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAAtt                                                                            >  2:102826/1‑150 (MQ=255)                                                                       tatTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTTTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTAc                                                                          >  1:27046/1‑151 (MQ=255)                                                                                  cAGTTTTGAGGACTATTTAGCGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATAcc                                                               <  2:259766/151‑1 (MQ=255)                                                                                    gTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCg                                                              >  2:330965/1‑150 (MQ=255)                                                                                     ttttGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGAt                                                            >  1:630768/1‑151 (MQ=255)                                                                                       ttGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGAttt                                                          >  2:778410/1‑151 (MQ=255)                                                                                        tGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGTGACCATCGCTAACAAATTACAAAGAGATACCGAtttt                                                         >  2:579122/1‑151 (MQ=255)                                                                                        tGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGAtttt                                                         >  2:831063/1‑151 (MQ=255)                                                                                        tGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGAttt                                                          <  1:539628/150‑1 (MQ=255)                                                                                         gAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTg                                                        <  1:189185/151‑1 (MQ=255)                                                                                         gAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTg                                                        <  2:211213/151‑1 (MQ=255)                                                                                           gaaCTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGAGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGa                                                       <  2:84660/148‑1 (MQ=255)                                                                                                   tAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATgtgg                                               <  1:38062/150‑1 (MQ=255)                                                                                                     aCGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATgtggtg                                             >  1:264573/1‑150 (MQ=255)                                                                                                     aCGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGa                                            <  1:665786/151‑1 (MQ=255)                                                                                                      cGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGAt                                           <  2:528638/151‑1 (MQ=255)                                                                                                       gCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGAtt                                          <  2:375632/151‑1 (MQ=255)                                                                                                         ttCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTgg                                        >  2:778934/1‑151 (MQ=255)                                                                                                           cATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTgggg                                      <  2:816311/151‑1 (MQ=255)                                                                                                           cATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTgggg                                      >  2:160317/1‑151 (MQ=255)                                                                                                                aCGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCg                                  <  1:495619/150‑1 (MQ=255)                                                                                                                aCGCGACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCg                                  <  1:381992/150‑1 (MQ=255)                                                                                                                    gACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACt                             >  1:125140/1‑151 (MQ=255)                                                                                                                    gACATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACt                             >  2:733416/1‑151 (MQ=255)                                                                                                                     aCATGGCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTc                            >  2:605770/1‑151 (MQ=255)                                                                                                                         ggCAATGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGAt                         <  1:455427/150‑1 (MQ=255)                                                                                                                            aaTGAATCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGATCCGc                     >  2:129478/1‑151 (MQ=255)                                                                                                                                aaTCAGGGGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGATCCGCCttt                 <  2:302478/151‑1 (MQ=255)                                                                                                                                aaTCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGATCCGCCttt                 >  1:758654/1‑151 (MQ=255)                                                                                                                                aaTCAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGATCCGCCtt                  >  1:703500/1‑150 (MQ=255)                                                                                                                                   cAGGCGAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGATCCGCCTTTAc               >  1:539756/1‑150 (MQ=255)                                                                                                                                        gAAAATTATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGATCCGCCTTTACCATCCt         >  1:519015/1‑151 (MQ=255)                                                                                                                                         aaaaTAATTGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGATCCGCCTTTACCATCCTa        <  2:200265/151‑1 (MQ=255)                                                                                                                                            attattGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGATCCGCCTTTACCATCCTAATc     >  2:682426/1‑151 (MQ=255)                                                                                                                                              tattGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGATCCGCCTTTACCATCCTAATcc    <  2:357567/150‑1 (MQ=255)                                                                                                                                               attGCCGCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGATCCGCCTTTACCATCCTAATCCgg  <  1:401973/151‑1 (MQ=255)                                                                                                                                                  gccgccAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGATCCGCCTTTACCATCCTAATCCgg  >  2:553091/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GGCATTTCCGCTGCGTATCTTTTTGCTGATTCTGGTGTTCTCAATATTTGTCATTGCAGCCCTGGCGCAATATTTTACGGCCAGTTTTGAGGACTATTTAACGCTTCATGTACGCGACATGGCAATGAATCAGGCGAAAATTATTGCCTCCAATGACAGTGTCATCTCTGCGGTGAAAACGCGTGACTACAAACGGCTGGCGACCATCGCTAACAAATTACAAAGAGATACCGATTTTGATTATGTGGTGATTGGGGACCGGCACTCGATCCGCCTTTACCATCCTAATCCGG  >  NC_000913/652282‑652574

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑