Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,428,923	G→A	pseudogene (128/189 nt)	lomR →	Rac prophage; protein LomR_1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,428,923	0	G	A	100.0%	255.2 / NA	80	pseudogene (128/189 nt)	lomR	Rac prophage; protein LomR_1
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (38/42);  total (38/42)

CCTTCCCAAGGGGAAAACGCACGACGTGCTGACCGGAAGTGATGACGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTA  >  NC_000913/1428789‑1429057                                                                                                                                       |                                                                                                                                       ccTTCCCAAGGGGAAAGCGCACTACGTGCTGCCCGAAAGAGATGACGGTCGCCACAGCAACACGTCTCTGGCGT‑GGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCgg                                                                                                                                  <  2:668280/140‑1 (MQ=255)  cTTCCCAAGGGGAAAACGCACGACGTGCTGACCGGAAGTGATGACGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCgg                                                                                                                                  >  2:836928/1‑140 (MQ=255)    tCCCAAGGGGAAAACGCACGACGTGCTGACCGGAAGTGATGACGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCa                                                                                                                                >  1:287948/1‑140 (MQ=255)          ggggAAAACGCACGACGTGCTGACCGGAAGTGATGACGGTCGCCACAGCAACACGTCTCTGGCGT‑GGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGACATTCCTCATGAAGCCCCCGGGAGTGGCgg                                                                                                                         >  1:668280/1‑139 (MQ=255)           gggAAAACGCACGACGCTCTGACCGGAAGTTATTAGGGTTGCCACAGCGACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGa                                                                                                                         <  2:712197/140‑1 (MQ=255)            ggAAAACGCACGACGTGCTGACCGGAAGTGATGACGGTCGCCACAGCAACACGTCTCTAGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGAc                                                                                                                        >  1:307200/1‑140 (MQ=255)              aaaaCGCACGACGTGCTGACCGGAAGTGATGACGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTg                                                                                                                      >  2:1240438/1‑140 (MQ=255)               aaaCGCACGACGTGCTGACCGGAAGTGATGACGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTgg                                                                                                                     <  2:1060092/140‑1 (MQ=255)                       gacgTGCTGACCGGAAGTGATGACGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGa                                                                                                             <  1:794362/140‑1 (MQ=255)                       gacgTGCTGACCGGAAGGGATGCCGGTCGCCACAGCAACGCGTCTCTGCCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGa                                                                                                             <  2:664691/140‑1 (MQ=255)                            gCTGACCGGAAGTGATGACGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCGCGGCAGCGGCGACTGGCGCACTGACGGtt                                                                                                        >  1:567821/1‑140 (MQ=255)                             cTGACCGGAAGTGATGACGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGttt                                                                                                       <  1:341416/140‑1 (MQ=255)                                       gTGATGACGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGgt                                                                                             <  1:1240438/140‑1 (MQ=255)                                        tgatgaCGGTCGCCACGCCAACGCGTCTCTTGCTTTGGGCGCTGGCGTCCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGgtg                                                                                            <  2:756321/140‑1 (MQ=255)                                        tgatgaCGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGgtg                                                                                            <  2:287948/140‑1 (MQ=255)                                          atgaCGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTc                                                                                          >  2:1165508/1‑140 (MQ=255)                                             aCGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGt                                                                                       <  2:1279865/140‑1 (MQ=255)                                                 tCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCg                                                                                         <  2:1143691/134‑1 (MQ=255)                                                 tCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCg                                                                                         >  1:1143691/1‑134 (MQ=255)                                                           aaCACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGCTTCATCGTGGGTGTCGGTAATAAGTTGTGGtt                                                                         >  2:1425717/1‑140 (MQ=255)                                                           aaCACGTCTCTGGCGGGGGGAGCTGACGTCCAGTTTATCCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGAtt                                                                         <  1:446314/140‑1 (MQ=255)                                                             cacGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAg                                                                       <  2:139593/140‑1 (MQ=255)                                                             cacGTCTCTGGCGTGGGAAGCTTGCGTGCAGATTAACCCAACCGAATCTGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAg                                                                       <  2:567821/140‑1 (MQ=255)                                                             cacGTCTCTAGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAg                                                                       <  2:307200/140‑1 (MQ=255)                                                               cGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGcc                                                                     <  2:1419006/140‑1 (MQ=255)                                                               cGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGcc                                                                     <  1:1501544/140‑1 (MQ=255)                                                               cGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGcc                                                                     <  2:30300/140‑1 (MQ=255)                                                                 tctcTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAg                                                                   <  1:846282/140‑1 (MQ=255)                                                                  ctctGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGGACTGACGGCTTCTTCGAGGGTGTCGGTTATAAGTTCTGATTAGCCAgg                                                                  >  2:446314/1‑140 (MQ=255)                                                                  ctctGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGTGTGCCGGTTATAAGTTCTGATTAGCCAgg                                                                  >  1:505071/1‑140 (MQ=255)                                                                       gCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAaca                                                             <  1:1046484/140‑1 (MQ=255)                                                                       gCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAaca                                                             <  1:271419/140‑1 (MQ=255)                                                                        ggTGGGGAGCTGGCGTGCAGTTGAACCCGACTGAATCCGGGGCCATTGCTATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAacac                                                            <  2:538051/139‑1 (MQ=255)                                                                        cGTGGGGAGCTGGCGTGGAGTTTAACCCGACCGAACCCGCGGCCATTGGTATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAacac                                                            <  2:473061/140‑1 (MQ=255)                                                                        cGTGGGGAGCTGGCGTGCAGTTTAATCCGTCCGAATCAGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAacac                                                            <  1:523748/140‑1 (MQ=255)                                                                        cGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAacac                                                            <  1:836928/140‑1 (MQ=255)                                                                        cGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAacac                                                            <  1:1165508/140‑1 (MQ=255)                                                                        cGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCAATGACGGTTTCATCGTGGGTGTCGGTTCTAAGTTCTGATTAGCCAGGTAacac                                                            <  2:1054741/140‑1 (MQ=255)                                                                        cGTGGGGAGCTGGCGTGCAGTTCAACCCGACTGAACCCGTGGCCATTGCTATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAacac                                                            <  2:505071/140‑1 (MQ=255)                                                                          tGGGGAGCTGGCGTGCAGTTTATCCCGACCGTATCCGTGACCATTGATGTTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAg                                                          <  2:490775/140‑1 (MQ=255)                                                                            gggAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAg                                                                       <  2:1495587/125‑1 (MQ=255)                                                                            gggAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAg                                                                       >  1:1495587/1‑125 (MQ=255)                                                                            gggAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAgg                                                                  <  2:1309623/130‑1 (MQ=255)                                                                            gggAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAgg                                                                  >  1:1309623/1‑130 (MQ=255)                                                                            gggAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAgtg                                                        >  2:1175159/1‑140 (MQ=255)                                                                            gggAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAgtg                                                        >  2:243793/1‑140 (MQ=255)                                                                            gggAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAgtg                                                        >  2:1532558/1‑140 (MQ=255)                                                                            gggAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGGGTCGGTTGTAACTTGTGATTATCCAGGTAACACAgtg                                                        >  1:728898/1‑140 (MQ=255)                                                                               aGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAacac                                                            <  1:977123/133‑1 (MQ=255)                                                                               aGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAacac                                                            >  2:977123/1‑133 (MQ=255)                                                                                 cTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATg                                                   >  1:294773/1‑140 (MQ=255)                                                                                      gTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGc                                              >  1:760221/1‑140 (MQ=255)                                                                                                      ccgaATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCGGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCgg                              >  2:24771/1‑140 (MQ=255)                                                                                                      ccgaATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTc                                    >  1:815474/1‑134 (MQ=255)                                                                                                      ccgaATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTc                                    <  2:815474/134‑1 (MQ=255)                                                                                                      ccgaATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCgg                              >  1:1117702/1‑140 (MQ=255)                                                                                                      ccgaATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCgg                              >  2:206930/1‑140 (MQ=255)                                                                                                      ccgaATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCgg                              >  2:1491510/1‑140 (MQ=255)                                                                                                      ccgaATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCgg                              >  1:312072/1‑140 (MQ=255)                                                                                                      ccgaATCCGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTATGATTAGCGAGGGAACACAGTGTTATGACAGCCCGCCGGTTCAGGCgg                              >  2:403368/1‑140 (MQ=255)                                                                                                            ccGTGGCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCtttt                        >  2:1115547/1‑140 (MQ=255)                                                                                                                ggCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGTGGGCTTTTTtgt                    <  2:760221/140‑1 (MQ=255)                                                                                                                ggCCATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTtgt                    >  2:1083210/1‑140 (MQ=255)                                                                                                                   cATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                 <  1:1083210/140‑1 (MQ=255)                                                                                                                   cATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                 <  2:1486252/140‑1 (MQ=255)                                                                                                                   cATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                 <  1:1115547/140‑1 (MQ=255)                                                                                                                   cATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                 <  1:1175159/140‑1 (MQ=255)                                                                                                                   cATTGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTAATGACAGGCCGGCGGTTCAGGTGGGCTTTTTTGTggg                 >  1:500603/1‑140 (MQ=255)                                                                                                                      tGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                 >  2:1151287/1‑137 (MQ=255)                                                                                                                      tGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                 <  1:1151287/137‑1 (MQ=255)                                                                                                                      tGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTg              >  2:12957/1‑140 (MQ=255)                                                                                                                      tGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTg              >  1:767427/1‑140 (MQ=255)                                                                                                                      tGATATTGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATAACAGCCCTCCGGTTCAGGCGGTCTTTTTTGTGGGGTg              >  1:635949/1‑140 (MQ=255)                                                                                                                            tGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                 <  2:1360591/131‑1 (MQ=255)                                                                                                                            tGCTTATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                 >  1:1360591/1‑131 (MQ=255)                                                                                                                                tATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                 <  2:1402971/127‑1 (MQ=255)                                                                                                                                tATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTggg                 >  1:1402971/1‑127 (MQ=255)                                                                                                                                tATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAg    >  1:992454/1‑140 (MQ=255)                                                                                                                                tATGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAg    >  2:1152947/1‑140 (MQ=255)                                                                                                                                  tGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGCCAGCCCGCAGGTCCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTa  >  2:613174/1‑140 (MQ=255)                                                                                                                                  tGAAGACCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTa  <  2:992454/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                       CCTTCCCAAGGGGAAAACGCACGACGTGCTGACCGGAAGTGATGACGGTCGCCACAGCAACACGTCTCTGGCGTGGGGAGCTGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTA  >  NC_000913/1428789‑1429057

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