Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,967,136	G→A	A101V (GCG→GTG)	cheY ←	chemotaxis protein CheY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,967,136	0	G	A	98.7%	255.0 / NA	78	A101V (GCG→GTG)	cheY	chemotaxis protein CheY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (42/35);  minor base C (0/1);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTC  >  NC_000913/1967004‑1967269                                                                                                                                     |                                                                                                                                        gCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTgcgc                                                                                                                                  <  2:1404344/140‑1 (MQ=255) gCTCGTCAGCAGGTTCGATTGATGGATGCTTCATAGTAGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTgcgc                                                                                                                                  <  1:798077/140‑1 (MQ=255)     gTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCgca                                                                                                                              >  2:102099/1‑140 (MQ=255)      tCAGCAGGTTTGATTGATGGTTCCATCATAGCCGCCTCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCgcag                                                                                                                             <  2:486249/140‑1 (MQ=255)        agcagGTTTGATTGATGGTTGCATCATAGTCGCATCCTCGCATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGTGCCgcagca                                                                                                                           >  2:1091709/1‑140 (MQ=255)        agcagGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCgcagca                                                                                                                           <  1:1693884/140‑1 (MQ=255)                  attgatGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATATCCACTGGACCCCACGTGGGCCCCAGCAATGATGTtct                                                                                                                 >  1:607490/1‑140 (MQ=255)                       tGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCtctttc                                                                                                            >  2:152768/1‑140 (MQ=255)                       tGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCtctttc                                                                                                            >  1:1262093/1‑140 (MQ=255)                       tGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCtctttc                                                                                                            >  1:1324149/1‑140 (MQ=255)                        ggTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCtctttct                                                                                                           >  2:1472104/1‑140 (MQ=255)                         gTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCtctttctt                                                                                                          <  1:391216/140‑1 (MQ=255)                         gTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCtctttctt                                                                                                          <  1:152768/140‑1 (MQ=255)                           tGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCGCCCACTTACGCGGCGGCAATGATGTTCTCTTTCTTCg                                                                                                        >  1:489757/1‑140 (MQ=255)                           tGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCg                                                                                                        >  1:1007680/1‑140 (MQ=255)                             catcatAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCt                                                                                                      >  1:40217/1‑140 (MQ=255)                              atcatAGTCGCATCCTCCCATGTCCAGTTTCACAATGATCTTGTTGAGCTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCtt                                                                                                     <  2:489757/140‑1 (MQ=255)                              atcatAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCtt                                                                                                     <  2:1007680/140‑1 (MQ=255)                                  tAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTg                                                                                                  >  2:1257505/1‑139 (MQ=255)                                  tAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTg                                                                                                  <  1:1257505/139‑1 (MQ=255)                                   aGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCtt                                                                                                     <  1:1715974/135‑1 (MQ=255)                                   aGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCtt                                                                                                     >  2:1715974/1‑135 (MQ=255)                                     tCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGt                                                                                              >  2:1310987/1‑140 (MQ=255)                                     tCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGt                                                                                              <  1:1091709/140‑1 (MQ=255)                                      cGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTc                                                                                             <  1:102099/140‑1 (MQ=255)                                      cGCATCCTCACATGCCCAGTTTATCCAAGTTTTTGTTCAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTc                                                                                             <  1:567444/140‑1 (MQ=255)                                        cATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTc                                                                                             <  2:1635390/138‑1 (MQ=255)                                        cATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTc                                                                                             >  1:1635390/1‑138 (MQ=255)                                          tCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAg                                                                                               <  1:400615/134‑1 (MQ=255)                                          tCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAg                                                                                               >  2:400615/1‑134 (MQ=255)                                            cTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCACTGAAGTTCTCTTGCTTCCCTTCTGCAGTCACCAtt                                                                                       >  2:612111/1‑140 (MQ=255)                                             tCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTc                                                                                             <  1:1028039/133‑1 (MQ=255)                                             tCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTc                                                                                             >  2:1028039/1‑133 (MQ=255)                                                   gCCCAGTTTCTCAAAAATTTTGCTGAGTTGTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTg                                                                              <  1:716956/140‑1 (MQ=255)                                                         tttCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCATTGATGTTCACTTGCTTCGCTTCTGCGGTCACCATTA‑‑ACACTGGCAATg                                                                        >  2:716956/1‑140 (MQ=255)                                                         tttCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATg                                                                        >  2:345500/1‑140 (MQ=255)                                                         tttCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATg                                                                        >  2:1069776/1‑140 (MQ=255)                                                                   aTTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATcgc                                                              >  1:477223/1‑140 (MQ=255)                                                                     tttGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCgcgc                                                            <  1:1069776/140‑1 (MQ=255)                                                                       tgttgAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCa                                                          >  1:1633549/1‑140 (MQ=255)                                                                       tgttgAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGA‑GTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCAt                                                         >  2:23650/1‑140 (MQ=255)                                                                            aGTTTTTCGTCCAGCGCGGCGGGGGTAAATGGCTTCACCGCATAGCCAGTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGc                                                     <  2:885231/140‑1 (MQ=255)                                                                             gTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGcc                                                    >  2:249128/1‑139 (MQ=255)                                                                             gTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCa                                                    >  2:253112/1‑140 (MQ=255)                                                                              tttttCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCAc                                                   <  1:1310987/140‑1 (MQ=255)                                                                                   cctccAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACCGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAAtt                                              >  2:1247352/1‑140 (MQ=255)                                                                                        aGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTTGCAATGCCGACGTCTCGCCATGCGCACGAATTGtttt                                         >  1:884191/1‑140 (MQ=255)                                                                                          cGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCa                                       >  2:948610/1‑140 (MQ=255)                                                                                           gTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAg                                      >  2:72952/1‑140 (MQ=255)                                                                                            tCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATGA‑‑ACACTGGAAATTCCGGCATCGCGCAATCCGCACGAATTGTTTTCAGc                                     >  2:575451/1‑140 (MQ=255)                                                                                                    gTAAATGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAg                             >  1:995929/1‑140 (MQ=255)                                                                                                    gTAAATGGCTTCACCACATAGCCACCGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAg                             <  1:1247352/140‑1 (MQ=255)                                                                                                     tAAATGGCTTCACCACATAGCCCCTGGCCCCCACTAGCGACGCAGCAACGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAgg                            <  2:607490/140‑1 (MQ=255)                                                                                                      aaaTGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAg                             <  1:1472104/140‑1 (MQ=255)                                                                                                      aaaTGGCTTCACCACATAGCCACTGG‑CCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGcc                          <  1:304941/140‑1 (MQ=255)                                                                                                        aTGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGcca                         <  2:40217/140‑1 (MQ=255)                                                                                                        aTGGCTTCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGcca                         <  1:253112/140‑1 (MQ=255)                                                                                                         tGGCTTCACCACATAGCCCCTGGCCCCCACTTGAGCCGCAGTAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGGGCCATCCGCACGAATTGTTTTCAGCATTTCCAGGccat                        <  2:1859336/140‑1 (MQ=255)                                                                                                          ggCTTCACCACATAGTCACTGGCCCCCACCTGCGCCGTAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGccatc                       <  1:575451/140‑1 (MQ=255)                                                                                                              tCACCACATAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCata                   <  2:1262093/140‑1 (MQ=255)                                                                                                              cagtcacaTAGCCACCGGCCCCTACCTGCTCCGCAGCCATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCata                   <  1:574112/136‑1 (MQ=255)                                                                                                               agccacATAGCCACTCGCCCCCCCTGGCGCCGCAGCACTGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCatat                  <  1:549750/138‑1 (MQ=255)                                                                                                                  cacaTAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTgg               >  2:1849092/1‑140 (MQ=255)                                                                                                                   acaTAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTggg              <  1:72952/140‑1 (MQ=255)                                                                                                                     aTAGCCACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCa            >  2:1876245/1‑140 (MQ=255)                                                                                                                         ccACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGtt        >  2:217688/1‑140 (MQ=255)                                                                                                                          cACTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTc       <  2:1324149/140‑1 (MQ=255)                                                                                                                            cTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGc             >  2:1155734/1‑132 (MQ=255)                                                                                                                            cTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGc             <  1:1155734/132‑1 (MQ=255)                                                                                                                            cTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCa            >  1:1065359/1‑133 (MQ=255)                                                                                                                            cTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCa            <  2:1065359/133‑1 (MQ=255)                                                                                                                            cTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGAACGAATTGGTTTGAGCAATTCCAGGCCATCGATATTGGGCATGTTcct     >  2:549750/1‑139 (MQ=255)                                                                                                                            cTGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGAACGAATTGATTTGAGCGATTCCAGGACATCAATATTGGGCATGTTctg     >  2:574112/1‑138 (MQ=255)                                                                                                                             tGGCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAg    >  2:1680632/1‑140 (MQ=255)                                                                                                                              ggCCCCCACTTGCGCCGTAGCACTGATGTTCCCTTACTTCGCTTCTGCTGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGt   <  2:728873/140‑1 (MQ=255)                                                                                                                               gCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGc             <  2:1908308/129‑1 (MQ=255)                                                                                                                               gCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGc             >  1:1908308/1‑129 (MQ=255)                                                                                                                               gCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTc  >  1:295750/1‑140 (MQ=255)                                                                                                                               gCCCCCACTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTc  >  1:85733/1‑140 (MQ=255)                                                                                                                                     |                                                                                                                                        GCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTA‑‑ACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGTTTTCAGCAATTCCAGGCCATCCATATTGGGCATGTTCCAGTC  >  NC_000913/1967004‑1967269

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