Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,026,812	(C)6→7	coding (204/459 nt)	mgsA ←	methylglyoxal synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,026,806	1	.	C	98.9%	331.3 / ‑4.3	91	coding (210/459 nt)	mgsA	methylglyoxal synthase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (47/43);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑A‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCACCCAGCTCATCA  >  NC_000913/1026669‑1026939                                                                                                                                            |                                                                                                                                      aCCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑A‑cc                                                                                                                                     >  1:480277/1‑140 (MQ=255)   cGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑A‑cccc                                                                                                                                   >  2:878373/1‑140 (MQ=255)        tGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑A‑cccccccat                                                                                                                              <  2:798469/140‑4 (MQ=255)          ttCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑A‑cccccccattg                                                                                                                            <  1:1267205/140‑6 (MQ=255)            ccATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTggg                                                                                                                           <  1:356338/140‑1 (MQ=255)            ccATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTggg                                                                                                                           <  1:782854/140‑1 (MQ=255)            ccATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTggg                                                                                                                           <  1:331562/140‑1 (MQ=255)            ccATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTggg                                                                                                                           <  1:236352/140‑1 (MQ=255)                aCCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCAc                                                                                                                       >  2:914155/1‑140 (MQ=255)                aCCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCAc                                                                                                                       >  2:1033546/1‑140 (MQ=255)                aCCGCCGCCAGACGCAGCAACGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTCCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCAc                                                                                                                       >  2:1252376/1‑140 (MQ=255)                    tCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCaa                                                                                                                   >  1:130471/1‑140 (MQ=255)                    tCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCaa                                                                                                                   >  2:1082960/1‑140 (MQ=255)                    tCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCaa                                                                                                                   >  2:1100486/1‑140 (MQ=255)                    tCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCaa                                                                                                                   >  2:563333/1‑140 (MQ=255)                    tCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCaa                                                                                                                   >  2:713693/1‑140 (MQ=255)                    tCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCaa                                                                                                                   >  2:1323862/1‑140 (MQ=255)                    tCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCaa                                                                                                                   >  2:898496/1‑140 (MQ=255)                    tCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCaa                                                                                                                   >  2:904524/1‑140 (MQ=255)                      gCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAGCCTGCTGGTCAA‑CCCCCCATTGGGCCACTCAACa                                                                                                                 <  1:1389574/140‑1 (MQ=255)                      gCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACa                                                                                                                 >  2:890006/1‑140 (MQ=255)                         aGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATcg                                                                                                              >  1:430321/1‑140 (MQ=255)                           aCGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATcgcg                                                                                                            >  2:744942/1‑140 (MQ=255)                            cGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGt                                                                                                           <  1:563333/140‑1 (MQ=255)                               agcaAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATcg                                                                                                              <  2:582014/134‑1 (MQ=255)                               agcaAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATcg                                                                                                              >  1:582014/1‑134 (MQ=255)                                      cTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCAt                                                                                                 >  2:642644/1‑140 (MQ=255)                                      cTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCAt                                                                                                 <  1:878373/140‑1 (MQ=255)                                                aGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGACCACTCAACATCGCGTTGACGTTCATGCCGGTcgcg                                                                                       >  1:222110/1‑140 (MQ=255)                                                 ggATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTcgcgc                                                                                      <  1:890006/140‑1 (MQ=255)                                                 ggATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTcgcgc                                                                                      <  1:85814/140‑1 (MQ=255)                                                 ggATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTcgcgc                                                                                      <  1:1100486/140‑1 (MQ=255)                                                  gATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTcgcgcg                                                                                     <  1:1249166/140‑1 (MQ=255)                                                     cGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGa                                                                                  >  1:130065/1‑140 (MQ=255)                                                       tGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaaa                                                                                >  2:529150/1‑140 (MQ=255)                                                       tGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaaa                                                                                >  2:389857/1‑140 (MQ=255)                                                       tGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaaa                                                                                >  2:161995/1‑140 (MQ=255)                                                       tGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaaa                                                                                >  2:1370640/1‑140 (MQ=255)                                                       tGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaaa                                                                                >  1:1303995/1‑140 (MQ=255)                                                       tGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGaaa                                                                                >  2:1074974/1‑140 (MQ=255)                                                         cggcacggcaTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAAtt                                                                              <  1:1033546/140‑1 (MQ=255)                                                         cggcacggcaTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAAtt                                                                              <  1:725629/140‑1 (MQ=255)                                                           gcacggcaTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTaa                                                                            <  1:351078/140‑1 (MQ=255)                                                            cacggcaTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAg                                                                           >  2:1221938/1‑140 (MQ=255)                                                              cggcaTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGtt                                                                         >  2:104487/1‑140 (MQ=255)                                                               ggcaTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTa                                                                        <  2:222110/140‑1 (MQ=255)                                                                 caTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTAcc                                                                      >  1:1001942/1‑140 (MQ=255)                                                                 caTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTAcc                                                                      >  2:182753/1‑140 (MQ=255)                                                                  aTTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCg                                                                     <  2:1076772/140‑1 (MQ=255)                                                                    ttAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGt                                                                   <  1:744942/140‑1 (MQ=255)                                                                      aGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAg                                                                 >  2:740884/1‑140 (MQ=255)                                                                      aGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAg                                                                 <  1:1082960/140‑1 (MQ=255)                                                                       gTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGt                                                                <  1:1298901/140‑1 (MQ=255)                                                                         ggATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTaa                                                                            >  1:611114/1‑126 (MQ=255)                                                                         ggATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTaa                                                                            <  2:611114/126‑1 (MQ=255)                                                                         ggATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGc                                                              <  2:230005/140‑1 (MQ=255)                                                                         ggATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGc                                                              <  1:1221938/140‑1 (MQ=255)                                                                         ggATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGc                                                              <  1:740884/140‑1 (MQ=255)                                                                         ggATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGc                                                              <  2:1016043/140‑1 (MQ=255)                                                                           aTCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCt                                                            >  1:797368/1‑140 (MQ=255)                                                                             cccAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTa                                                                        >  1:216307/1‑126 (MQ=255)                                                                             cccAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTa                                                                        <  2:216307/126‑1 (MQ=255)                                                                              ccAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGtt                                                         <  2:1001942/140‑1 (MQ=255)                                                                              ccAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGtt                                                         >  1:317719/1‑140 (MQ=255)                                                                                agaagaAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTAc                                                                       >  1:842195/1‑124 (MQ=255)                                                                                agaagaAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTAc                                                                       <  2:842195/124‑1 (MQ=255)                                                                                agaagaAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCgg                                                                    >  1:1329322/1‑127 (MQ=255)                                                                                agaagaAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCgg                                                                    <  2:1329322/127‑1 (MQ=255)                                                                                 gaagaaAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCa                                                      <  1:1189193/140‑1 (MQ=255)                                                                                  aagaaAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGTCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCAt                                                     >  1:867766/1‑140 (MQ=255)                                                                                  aagaaAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCAt                                                     <  1:643483/140‑1 (MQ=255)                                                                                  aagaaAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCAt                                                     <  1:666866/140‑1 (MQ=255)                                                                                   agaaAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATa                                                    <  1:104487/140‑1 (MQ=255)                                                                                   agaaAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATa                                                    <  2:1303995/140‑1 (MQ=255)                                                                                    gaaAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGTCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATAc                                                   <  2:867766/140‑1 (MQ=255)                                                                                    gaaAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATAc                                                   >  2:1393507/1‑140 (MQ=255)                                                                                        aTCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTa                                               <  1:904524/140‑1 (MQ=255)                                                                                         tCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCt                                                            >  2:662567/1‑126 (MQ=255)                                                                                         tCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCt                                                            <  1:662567/126‑1 (MQ=255)                                                                                            aTACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACgtg                                           <  1:898496/140‑1 (MQ=255)                                                                                               cATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGtgttg                                        >  1:1058286/1‑140 (MQ=255)                                                                                                     ttttCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAg                                  >  2:491602/1‑140 (MQ=255)                                                                                                          ccTTCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACg                             >  2:262451/1‑140 (MQ=255)                                                                                                             tCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTaa                               <  2:295181/135‑1 (MQ=255)                                                                                                             tCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTaa                               >  1:295181/1‑135 (MQ=255)                                                                                                             tCTGAGATCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGtt                          >  1:763051/1‑140 (MQ=255)                                                                                                                 agaTCAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATg                      <  1:642644/140‑1 (MQ=255)                                                                                                                     cAATGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGt                  <  1:504490/140‑1 (MQ=255)                                                                                                                        tGCGCCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTcc               >  1:1081669/1‑140 (MQ=255)                                                                                                                         gcgcCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCa              >  1:666640/1‑140 (MQ=255)                                                                                                                           gcCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCAcc            >  1:1171711/1‑140 (MQ=255)                                                                                                                           gcCAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCAcc            <  2:130471/140‑1 (MQ=255)                                                                                                                            ccAACCTGCTGGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCTGTTCCAccc           <  1:713693/140‑1 (MQ=255)                                                                                                                                   gctgGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCACCCAGCtcat    <  1:529150/140‑1 (MQ=255)                                                                                                                                    ctgGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCACCCAGCtcatc   >  1:491408/1‑140 (MQ=255)                                                                                                                                    ctgGTC‑ACCCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCC‑GTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCACCCAGCtcatca  >  1:32120/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                      ACCGGAATGTTCCATACCGTCGCCAGACGCAGCAAGGCTTTCACGTCAGGATCGTGCGGCACGGCATTTAGTGGATCCCAGAAGAAAATCAATACATCAATTTTCCCTTCTGAGATCAATGCGCCAACCTGCTGGTC‑A‑CCCCCCATTGGGCCACTCAACATCGCGTTGACGTTCATGCCGGTCGCGCGGGAAATTAAGTTACCGGTAGTGCCTGTTGCATACAGTACGTGTTGTTCCAGTAACGGTTGATGCCGTTCCACCCAGCTCATCA  >  NC_000913/1026669‑1026939

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