Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,842,330	T→C	R25G (AGA→GGA)	yicN ←	DUF1198 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,842,330	0	T	C	100.0%	232.4 / NA	65	R25G (AGA→GGA)	yicN	DUF1198 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (30/35);  total (30/35)

GCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTC  >  NC_000913/3842183‑3842476                                                                                                                                                    |                                                                                                                                                   gcgcGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCgg                                                                                                                                                  <  1:877837/150‑1 (MQ=255)     ggcgTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGcccca                                                                                                                                             <  2:921883/151‑2 (MQ=255)            tGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTc                                                                                                                                       <  1:15273/150‑1 (MQ=255)               aGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGc                                                                                                                                    >  1:964202/1‑150 (MQ=255)               aGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGc                                                                                                                                    >  2:601072/1‑150 (MQ=255)                ggTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCAc                                                                                                                                  <  1:650377/151‑1 (MQ=255)                   gCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCg                                                                                                                               <  1:88866/151‑1 (MQ=255)                     gACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGaa                                                                                                                             <  1:855746/151‑1 (MQ=255)                        tCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAaac                                                                                                                          >  1:906570/1‑151 (MQ=255)                            ccGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAaaccaac                                                                                                                      >  1:371757/1‑151 (MQ=255)                                     gCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAaca                                                                                                              <  2:1004781/150‑1 (MQ=255)                                       aaaTAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACt                                                                                                            >  1:62938/1‑150 (MQ=255)                                       aaaTAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACt                                                                                                            >  1:1032555/1‑150 (MQ=255)                                       aaaTAACGGAAAAATTCGACACAAGCTGAATTTCCCATTTGATCGACCATCGACTCCACGGGCACAACATCGATGTTCAGCCGATAGCTGAGACGGCGAATCGCTTTCCCGATCCCGGATGACAGCACCCGAAAACCAACCACAAaca                                                                                                              >  2:149894/1‑148 (MQ=255)                                           aaCGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGc                                                                                                       >  2:869775/1‑151 (MQ=255)                                             cGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTCCGACATCGATGTTCAGCCGATCGCCGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGcc                                                                                                      >  1:297923/1‑150 (MQ=255)                                              gTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGcc                                                                                                      >  1:669430/1‑149 (MQ=255)                                                aaaaaTTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCg                                                                                                  <  2:835864/151‑1 (MQ=255)                                                aaaaaTTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCg                                                                                                  <  1:921791/151‑1 (MQ=255)                                                       cGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCg                                                                                             <  1:663266/149‑1 (MQ=255)                                                        gTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGc                                                                                          >  1:808898/1‑151 (MQ=255)                                                              accTGACTTTCCCATTTGATCGACCATCGACTCAACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTCCCGGCCCCGGAGGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTAt                                                                                     <  2:758333/148‑1 (MQ=255)                                                              aGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATc                                                                                    <  1:156571/151‑1 (MQ=255)                                                                 tGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCCACCACAAACACTACCGCCAGCGCGGCGCGCATTATcc                                                                                   >  2:435144/1‑149 (MQ=255)                                                                   aCTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCa                                                                                  <  2:505723/148‑1 (MQ=255)                                                                     tttCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATc                                                                              >  1:739483/1‑150 (MQ=255)                                                                      ttCCCATTTGATCGACCATCGACTCCCCGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCa                                                                             <  1:263930/150‑1 (MQ=255)                                                                      ttCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCAt                                                                            >  1:427689/1‑151 (MQ=255)                                                                      ttCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCAt                                                                            >  2:332197/1‑151 (MQ=255)                                                                        cccATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGa                                                                          <  2:488039/151‑1 (MQ=255)                                                                         ccATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATg                                                                           >  2:1011325/1‑149 (MQ=255)                                                                           aTTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGtt                                                                       >  1:10468/1‑151 (MQ=255)                                                                                 acGACCATCGACTCCCGGGGTACGACATCGATGTTCAGCCGATCGCTCAGACGGCGAATCGCTTTTCCGTCCCCGGTTGTCAGCACCCGAAAACCAACCAAAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTg                                                                 <  2:438139/150‑1 (MQ=255)                                                                                     ccATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGt                                                              <  2:973035/150‑1 (MQ=255)                                                                                     ccATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGt                                                              <  1:13097/150‑1 (MQ=255)                                                                                        tCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATc                                                          >  1:271242/1‑151 (MQ=255)                                                                                         cGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATc                                                          <  1:285649/150‑1 (MQ=255)                                                                                         cGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCa                                                         <  1:404482/151‑1 (MQ=255)                                                                                         cGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCa                                                         <  1:834542/151‑1 (MQ=255)                                                                                          gacgacACGGGTACGACATCGATGGCAAGACGATCGCTGAGACGGCGAAAGGGCAGGCCGGCAACGGAGGGCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATAAGGAGTTCTCCCGCGTAATc                                                          >  1:780347/6‑149 (MQ=255)                                                                                          gACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGCCGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCCCAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAg                                                        <  1:96310/151‑1 (MQ=255)                                                                                           aCTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATc                                                          >  2:460164/1‑148 (MQ=255)                                                                                                  gggTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATg                                                >  1:686177/1‑151 (MQ=255)                                                                                                    gTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGat                                              <  1:658976/151‑1 (MQ=255)                                                                                                    gTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGat                                              >  2:417124/1‑151 (MQ=255)                                                                                                      aCGACATCGATGTTCAGCCGATCGCTGAGACGGCTAATCGCTTTTCCGGCCCCGGTTGTCAGCTCCCGAAAACCAACCACAACCACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGat                                              <  2:775688/149‑1 (MQ=255)                                                                                                             cGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTCCCGGCCCCGGATGTCAGCACCCGAAAACAACCCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGaaa                                     <  1:383367/151‑1 (MQ=255)                                                                                                              gATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGaaa                                     <  1:645618/150‑1 (MQ=255)                                                                                                                tGTTCAGCCGATCGCTGGGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAACCCAACCACAAACCCTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAAc                                    <  2:412007/149‑1 (MQ=255)                                                                                                                 gTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCt                                 <  1:960583/151‑1 (MQ=255)                                                                                                                  ttCAGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTc                                >  2:983291/1‑151 (MQ=255)                                                                                                                   tCAGCCGATCGCTGAGACTGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGCTCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTc                                >  2:880144/1‑150 (MQ=255)                                                                                                                     aGCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGcc                             >  1:519669/1‑151 (MQ=255)                                                                                                                      gCCGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGcc                             >  2:767913/1‑150 (MQ=255)                                                                                                                        cGATCGCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGt                          <  1:478707/151‑1 (MQ=255)                                                                                                                             gCTGAGACGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTa                     <  2:19829/151‑1 (MQ=255)                                                                                                                                   aCGGCGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATccac               >  2:658503/1‑151 (MQ=255)                                                                                                                                       cGAATCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATccacc              >  1:816096/1‑148 (MQ=255)                                                                                                                                          aTCGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGaata        <  1:806505/151‑1 (MQ=255)                                                                                                                                            cGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGaataat      >  1:282739/1‑151 (MQ=255)                                                                                                                                            cGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGaataat      <  1:683209/151‑1 (MQ=255)                                                                                                                                            cGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGaataat      <  1:982889/151‑1 (MQ=255)                                                                                                                                            cGCTTTTCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGaataa       <  2:739442/150‑1 (MQ=255)                                                                                                                                                tttCCGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTc  >  2:308857/1‑151 (MQ=255)                                                                                                                                                 ttCCGGCTCCGGATTTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCt    <  2:247685/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCGCGGCGTTTTGCAGGTGCGACTCATCCGGACGATGCAAATAACGTAAAAATTCGTCACCAGCTGACTTTCCCATTTGATCGACCATCGACTCCACGGGTACGACATCGATGTTCAGCCGATCGCTGAGACGGCGAATCGCTTTTCTGGCCCCGGATGTCAGCACCCGAAAACCAACCACAAACACTACCGCCAGCGTGGCGAGCATTATCCAGATCATGAGTTCTCCTGCGTAATCAGAACTCATGATATCTGGAAACCCTCGCCAGTTACTAATCCACCAGAATAATCTTC  >  NC_000913/3842183‑3842476

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