Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,139,301	A→G	5.6%	E467G (GAA→GGA)	yegH →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,139,301	0	A	G	5.6%	219.2 / 3.0	71	E467G (GAA→GGA)	yegH	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (33/34);  new base G (2/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

CTGAAACCATTGCCGGTAACTTACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTGCAGAAGGTGCAGATTATACC  >  NC_000913/2139154‑2139448                                                                                                                                                    |                                                                                                                                                    cTGAAACCATTGCCGGTAACTTACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGa                                                                                                                                                     >  4:124981/1‑148 (MQ=255)       ccATTGCCGGTAACTTACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGt                                                                                                                                              >  1:105001/1‑149 (MQ=255)       ccATTGCCGGTAACTTACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGt                                                                                                                                              <  3:280019/149‑1 (MQ=255)         aTTGCCGGTAACTTACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGa                                                                                                                                            <  3:124981/149‑1 (MQ=255)              cGGTAACTTACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAatat                                                                                                                                        >  1:149728/1‑148 (MQ=255)               ggTAACTTACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATca                                                                                                                                      <  2:252380/149‑1 (MQ=255)               ggTAACTTACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGAGCCGCTGCAGCTGGATGAAAAACGTGAATATca                                                                                                                                      >  6:71446/1‑149 (MQ=255)                gTAACTTACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATcac                                                                                                                                     <  5:10585/149‑1 (MQ=255)                gTAACTTACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATcac                                                                                                                                     <  6:285290/149‑1 (MQ=255)                 tAACTTACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCATAAGAATGCGGACGGTTCCTGGACGGCGAATGTTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATcaca                                                                                                                                    <  2:173632/149‑1 (MQ=255)                    cTTACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCa                                                                                                                                 >  1:54393/1‑149 (MQ=255)                     ttACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCAt                                                                                                                                >  2:298777/1‑149 (MQ=255)                        ccGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCAGTGc                                                                                                                             >  8:77829/1‑149 (MQ=255)                          gaacgaaGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGGAAAACGTGAATATCACACCATTGCCg                                                                                                                           >  1:28723/1‑149 (MQ=255)                          gaacgaaGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGGAAAACGTGAATATCACACCATTGCCg                                                                                                                           >  1:23629/1‑149 (MQ=255)                                aGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGt                                                                                                                     <  8:190819/149‑1 (MQ=255)                                  tGGAAGAGATCGACGCCCGGCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTg                                                                                                                   <  3:290023/149‑1 (MQ=255)                                     aaGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATg                                                                                                                <  2:105001/149‑1 (MQ=255)                                       gagaTCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGGAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGa                                                                                                              <  2:23629/149‑1 (MQ=255)                                       gagaTCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGGAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGa                                                                                                              <  2:28723/149‑1 (MQ=255)                                       gagaTCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGa                                                                                                              <  4:242760/149‑1 (MQ=255)                                          aTCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGa                                                                                                              >  2:75865/1‑146 (MQ=255)                                          aTCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGa                                                                                                              <  1:75865/146‑1 (MQ=255)                                             gACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTAtt                                                                                                         <  4:150645/148‑1 (MQ=255)                                              aCGCCCGCCCTGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTTGAGGATCTGGTGCAATATGTGGCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTg                                                                                                       <  5:71446/149‑1 (MQ=255)                                                   cgccATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCg                                                                                                  >  7:54236/1‑149 (MQ=255)                                                    gccATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCg                                                                                                  >  6:226091/1‑148 (MQ=255)                                                     ccATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTa                                                                                                <  2:54393/149‑1 (MQ=255)                                                        tGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTc                                                                                             <  2:149728/149‑1 (MQ=255)                                                                agaagaATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCt                                                                                     >  4:306092/1‑149 (MQ=255)                                                                agaagaATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCt                                                                                     <  4:157697/149‑1 (MQ=255)                                                                 gaagaaTGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTg                                                                                    <  4:123150/149‑1 (MQ=255)                                                                  aagaaTGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTgg                                                                                   <  5:226091/149‑1 (MQ=255)                                                                    gaaTGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGa                                                                                                                                            <  5:299517/90‑1 (MQ=255)                                                                    gaaTGCGGAAGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGa                                                                                                                                            >  6:299517/1‑90 (MQ=255)                                                                        gCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCgaag                                                                              >  1:275475/1‑148 (MQ=255)                                                                         cGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTTCCGCTGGATGAAAAACCTTAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCgaagaa                                                                            <  5:57286/149‑1 (MQ=255)                                                                               gTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGa                                                                                                              <  1:280003/109‑1 (MQ=255)                                                                               gTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGa                                                                                                              >  2:280003/1‑109 (MQ=255)                                                                                ttCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCAGCTGCCGCTGGATGAAACACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAgg                                                                     <  6:283047/149‑1 (MQ=255)                                                                                          gCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGatta                                                           >  2:261114/1‑149 (MQ=255)                                                                                             aaTGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTAttt                                                        >  7:266530/1‑149 (MQ=255)                                                                                               tGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGc                                                      >  6:151594/1‑149 (MQ=255)                                                                                                 gTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCTTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCt                                                     >  8:170652/1‑148 (MQ=255)                                                                                                 gTCATATGCCGATGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCt                                                     >  5:69543/1‑148 (MQ=255)                                                                                                  tCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCtt                                                    >  4:309321/1‑148 (MQ=255)                                                                                                  ccATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGGTGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCtt                                                    >  3:260631/2‑148 (MQ=255)                                                                                                     tatGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTaaaa                                                <  2:275475/149‑1 (MQ=255)                                                                                                       tGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACg                                              >  2:305768/1‑149 (MQ=255)                                                                                                            cTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAAAGTTg                                           >  1:64758/1‑147 (MQ=255)                                                                                                             tGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTTCACCGTATTACAAAACCTGGCGAAGAAGTTAAGGTGGGGGATTATTTTCTTAAAACGTTGCAg                                        >  7:165330/1‑149 (MQ=255)                                                                                                             tGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAg                                        >  8:25023/1‑149 (MQ=255)                                                                                                               gaggaTCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGt                                      <  1:261114/149‑1 (MQ=255)                                                                                                                aggaTCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTa                                     <  1:298777/149‑1 (MQ=255)                                                                                                                      tGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGaaa                                 <  5:69619/147‑1 (MQ=255)                                                                                                                      tGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGaaa                                 >  6:69619/1‑147 (MQ=255)                                                                                                                          gCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATc                           >  8:85408/1‑149 (MQ=255)                                                                                                                              tatGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGt                       <  8:54236/149‑1 (MQ=255)                                                                                                                               atGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTg                      >  1:339500/1‑149 (MQ=255)                                                                                                                                tgtgCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTg                      <  2:339500/148‑1 (MQ=255)                                                                                                                                    ccgctgccgctgGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTGCAGaa                 <  3:309321/149‑1 (MQ=255)                                                                                                                                     cgctgccgctgGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTGCAGAAg                <  6:69543/149‑1 (MQ=255)                                                                                                                                      gctgccgctgGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTCGCAGCGTATTCAAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTGCAGAAgg               >  2:193832/1‑149 (MQ=255)                                                                                                                                          ccgctgGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTGCAGAAGGTGc            <  8:266530/148‑1 (MQ=255)                                                                                                                                            gctgGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTGCAGAAGGTg             >  7:136824/1‑145 (MQ=255)                                                                                                                                            gctgGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTGCAGAAGGTGCAGa         <  5:151594/149‑1 (MQ=255)                                                                                                                                            gctgGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTGCAGAAGGTGCAGa         <  1:243399/149‑1 (MQ=255)                                                                                                                                                gATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTGCAGAAGGTGCAGATta      >  3:148397/1‑148 (MQ=255)                                                                                                                                                  tGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTaaaa                                                <  6:321653/104‑1 (MQ=255)                                                                                                                                                  tGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTaaaa                                                >  5:321653/1‑104 (MQ=255)                                                                                                                                                  tGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTGCAGAAGGTGCAGATTATAc   >  3:280463/1‑149 (MQ=255)                                                                                                                                                   gAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTGCAGAAGGTGCAGATTATAcc  >  4:131013/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CTGAAACCATTGCCGGTAACTTACCGAACGAAGTGGAAGAGATCGACGCCCGCCATGATATTCAGAAGAATGCGGACGGTTCCTGGACGGCGAATGGTCATATGCCGCTGGAGGATCTGGTGCAATATGTGCCGCTGCCGCTGGATGAAAAACGTGAATATCACACCATTGCCGGGCTGTTGATGGAGTATTTGCAGCGTATTCCAAAACCTGGCGAAGAAGTTCAGGTGGGGGATTATTTGCTTAAAACGTTGCAGGTAGAAAGCCATCGCGTGCAGAAGGTGCAGATTATACC  >  NC_000913/2139154‑2139448

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