Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,831,457	(G)5→4	coding (1001/1710 nt)	yicH →	putative inner membrane‑anchored periplasmic AsmA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,831,453	0	G	.	100.0%	296.4 / NA	71	coding (997/1710 nt)	yicH	putative inner membrane‑anchored periplasmic AsmA family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (36/35);  total (36/35)

GCAGACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCG  >  NC_000913/3831317‑3831585                                                                                                                                         |                                                                                                                                     gCAGACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAggg                                                                                                                                    >  2:648901/1‑139 (MQ=255)    gACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAggggta                                                                                                                                 <  1:648901/139‑3 (MQ=255)    gACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAggggta                                                                                                                                 <  1:1324397/139‑3 (MQ=255)        caGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGt                                                                                                                            <  1:652669/139‑1 (MQ=255)           gAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAg                                                                                                                         <  1:15902/139‑1 (MQ=255)           gAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAg                                                                                                                         >  1:8255/1‑139 (MQ=255)                   aCTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAg                                                                                                                 >  1:648658/1‑139 (MQ=255)                     tGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTc                                                                                                                   >  2:832442/1‑135 (MQ=255)                     tGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTc                                                                                                                   <  1:832442/135‑1 (MQ=255)                                   aGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGa                                                                                                 <  2:665955/139‑1 (MQ=255)                                        gTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGa                                                                                            <  1:1265972/139‑1 (MQ=255)                                         ttCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGa                                                                                                 <  2:714169/133‑1 (MQ=255)                                         ttCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGa                                                                                                 >  1:714169/1‑133 (MQ=255)                                         ttCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGCGACGGGaa                                                                                           >  1:226832/1‑139 (MQ=255)                                               tATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGtt                                                                                     <  2:1336676/139‑1 (MQ=255)                                               tATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGtt                                                                                     <  1:264176/139‑1 (MQ=255)                                                aTGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCTGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGCGACGGGAAAACGTTg                                                                                    <  2:1008178/139‑1 (MQ=255)                                                 tGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGCGACGGGAAAACGTTGa                                                                                   <  2:226832/139‑1 (MQ=255)                                                      cGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTg                                                                                                  <  2:1355585/119‑1 (MQ=255)                                                      cGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTg                                                                                                  >  1:1355585/1‑119 (MQ=255)                                                      cGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCtt                                                                              >  1:1210326/1‑139 (MQ=255)                                                        cTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTtga                                                                            >  2:215365/1‑139 (MQ=255)                                                                atttGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCgg                                                                    >  2:112494/1‑139 (MQ=255)                                                                     aCCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGaa                                                                                                                       >  1:1347388/1‑83 (MQ=255)                                                                     aCCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGaa                                                                                                                       <  2:1347388/83‑1 (MQ=255)                                                                     aCCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAAtt                                                               >  2:218321/1‑139 (MQ=255)                                                                      cccGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTg                                                              >  2:1205019/1‑139 (MQ=255)                                                                       ccGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGc                                                             <  1:215365/139‑1 (MQ=255)                                                                          aTTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAg                                                                                                                         <  2:802532/76‑1 (MQ=255)                                                                          aTTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAg                                                                                                                         >  1:802532/1‑76 (MQ=255)                                                                          aTTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCgg                                                          >  2:656800/1‑139 (MQ=255)                                                                           ttATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCggg                                                         <  1:617340/139‑1 (MQ=255)                                                                           ttATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCggg                                                         <  1:1205019/139‑1 (MQ=255)                                                                           ttATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCggg                                                         >  2:1087726/1‑139 (MQ=255)                                                                           ttATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCggg                                                         <  1:218321/139‑1 (MQ=255)                                                                            tataAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGc                                                        <  2:8255/139‑1 (MQ=255)                                                                               aaaCACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTc                                                                                                                   >  1:195923/1‑77 (MQ=255)                                                                               aaaCACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTc                                                                                                                   >  1:878706/1‑77 (MQ=255)                                                                               aaaCACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTc                                                                                                                   <  2:195923/77‑1 (MQ=255)                                                                               aaaCACGCACTTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTc                                                                                                                   <  2:878706/77‑1 (MQ=255)                                                                                aaCACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGa                                                    >  2:1296498/1‑139 (MQ=255)                                                                                     ggAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAg                                                                                                                         <  2:188978/65‑1 (MQ=255)                                                                                     ggAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAg                                                                                                                         >  1:188978/1‑65 (MQ=255)                                                                                     ggAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGa                                                                                            >  1:60908/1‑94 (MQ=255)                                                                                     ggAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGa                                                                                            <  2:60908/94‑1 (MQ=255)                                                                                     ggAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAAtata                                               >  2:748036/1‑139 (MQ=255)                                                                                     ggAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAAtata                                               >  1:1476435/1‑139 (MQ=255)                                                                                     ggAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAAtata                                               >  2:193310/1‑139 (MQ=255)                                                                                     ggAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAAtata                                               >  2:1318677/1‑139 (MQ=255)                                                                                     ggAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAAtata                                               >  2:1106684/1‑139 (MQ=255)                                                                                     ggAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAAtata                                               >  1:51795/1‑139 (MQ=255)                                                                                     ggAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAAtata                                               >  1:119670/1‑139 (MQ=255)                                                                                        aTTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCt                                            <  2:1210326/139‑1 (MQ=255)                                                                                             ccccGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCg                                       <  2:51795/139‑1 (MQ=255)                                                                                                       cGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGc                             >  1:1311526/1‑139 (MQ=255)                                                                                                         tAGCGCTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCaa                           <  1:1087726/139‑1 (MQ=255)                                                                                                            cgcTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCggg                                                         <  2:657016/106‑1 (MQ=255)                                                                                                            cgcTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCggg                                                         >  1:657016/1‑106 (MQ=255)                                                                                                              cTGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGtt                      <  1:1133171/139‑1 (MQ=255)                                                                                                               tGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGc                                                        >  2:1271956/1‑104 (MQ=255)                                                                                                               tGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGc                                                        <  1:1271956/104‑1 (MQ=255)                                                                                                               tGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTtg                     >  2:1468470/1‑139 (MQ=255)                                                                                                               tGCGCCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTtg                     >  2:1368435/1‑139 (MQ=255)                                                                                                                 cgcCAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTtgtg                   >  1:736461/1‑139 (MQ=255)                                                                                                                   ccAGTTCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGtgga                 <  2:1347013/139‑1 (MQ=255)                                                                                                                        tCACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGaa            <  2:1236806/139‑1 (MQ=255)                                                                                                                         cACCAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGaaa           <  1:619784/139‑1 (MQ=255)                                                                                                                            cAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCaa                           <  1:956502/120‑1 (MQ=255)                                                                                                                            cAGCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCaa                           >  2:956502/1‑120 (MQ=255)                                                                                                                              gCCGCTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGaca      >  1:80840/1‑139 (MQ=255)                                                                                                                                  cTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCg  <  2:119670/139‑1 (MQ=255)                                                                                                                                  cTGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCg  <  1:1296498/139‑1 (MQ=255)                                                                                                                                   tGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGc                                                        <  2:727085/84‑1 (MQ=255)                                                                                                                                   tGGGAA‑GGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGc                                                        >  1:727085/1‑84 (MQ=255)                                                                                                                                         |                                                                                                                                     GCAGACACAGGAAGGCAAACTGTCGATGAACGCTAGCGAGTTCATTTATGGTTCGCTGCATTTATTTGACCCGATTATAAACACGGAATTTTCCCCGCAGGGCGTAGCGCTGCGCCAGTTCACCAGCCGCTGGGAAGGGGGTATGGTCAGAACGTCAGGGAACTGGCTGCGTGACGGGAAAACGTTGATCCTTGATGATGCGGCAATTGCCGGGCTGGAATATACCTTGCCGAAAAACTGGCAACAGTTGTGGATGGAAACGACACCCG  >  NC_000913/3831317‑3831585

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