Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,177,332	C→T	5.8%	intergenic (‑128/+180)	gatY ← / ← fbaB	D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,177,332	0	C	T	5.8%	287.1 / 5.7	86	intergenic (‑128/+180)	gatY/fbaB	D‑tagatose 1,6‑bisphosphate aldolase 2, catalytic subunit/fructose‑bisphosphate aldolase class I
Reads supporting (aligned to +/- strand):  ref base C (43/38);  new base T (3/2);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTC  >  NC_000913/2177184‑2177480                                                                                                                                                      |                                                                                                                                                     cTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGcc                                                                                                                                                      >  1:340231/1‑149 (MQ=255) cTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGcc                                                                                                                                                      >  4:1012296/1‑149 (MQ=255)    tGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTcc                                                                                                                                                   <  2:1075091/149‑1 (MQ=255)       cGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGa                                                                                                                                                <  3:311118/149‑1 (MQ=255)       cGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGa                                                                                                                                                <  4:171611/149‑1 (MQ=255)          tACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAg                                                                                                                                             <  3:130119/149‑1 (MQ=255)           aCCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAgg                                                                                                                                            <  2:51467/149‑1 (MQ=255)              aCGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTCGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGaa                                                                                                                                         <  4:590096/149‑1 (MQ=255)               cGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAg                                                                                                                                        <  4:582725/149‑1 (MQ=255)                gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCTTCCCTATCGGGAGGc                                                                                                                                       >  2:1113262/1‑149 (MQ=255)                gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGc                                                                                                                                       >  2:695267/1‑149 (MQ=255)                gTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGc                                                                                                                                       <  3:63232/149‑1 (MQ=255)                 tACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCTTCCCTATCGGGAGGcc                                                                                                                                      <  1:1113262/149‑1 (MQ=255)                  aCATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCg                                                                                                                                     >  2:970896/1‑149 (MQ=255)                    aTTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATAGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                                                                   <  4:190589/149‑1 (MQ=255)                     ttttCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAg                                                                                                                                  <  4:1124422/149‑1 (MQ=255)                           tatCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGt                                                                                                                            <  1:357037/149‑1 (MQ=255)                                tGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGa                                                                                                                                                <  1:773423/124‑1 (MQ=255)                                tGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGa                                                                                                                                                >  4:853048/1‑124 (MQ=255)                                tGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGa                                                                                                                                                <  3:853048/124‑1 (MQ=255)                                tGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGa                                                                                                                                                >  2:773423/1‑124 (MQ=255)                                    gtttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGtt                                                                                                                   <  1:507141/149‑1 (MQ=255)                                     tttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAgg                                                                                                                                            >  3:716463/1‑123 (MQ=255)                                     tttgtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAgg                                                                                                                                            <  4:716463/123‑1 (MQ=255)                                        gtttTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAATTTCGCGTGTTTTTa                                                                                                               <  3:85260/149‑1 (MQ=255)                                         ttttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTAt                                                                                                              >  2:964939/1‑149 (MQ=255)                                           ttCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGc                                                                                                            <  2:340231/149‑1 (MQ=255)                                                  tCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGc                                                                                                     >  3:638579/1‑149 (MQ=255)                                                           aaTGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAAc                                                                                            >  2:504654/1‑149 (MQ=255)                                                             tGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACAt                                                                                          >  3:874079/1‑149 (MQ=255)                                                              gAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATa                                                                                         >  2:325600/1‑149 (MQ=255)                                                               aaaTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                                                                        >  2:111001/1‑149 (MQ=255)                                                               aaaTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                                                                        <  1:695267/149‑1 (MQ=255)                                                                aaTTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACg                                                                                       >  4:922274/1‑149 (MQ=255)                                                                 attattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGt                                                                                      <  3:922274/149‑1 (MQ=255)                                                                  ttattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                                                                   <  3:948202/103‑1 (MQ=255)                                                                  ttattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTa                                                                                                                                   >  4:948202/1‑103 (MQ=255)                                                                  ttattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTg                                                                                     >  3:409295/1‑149 (MQ=255)                                                                   tattTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTgg                                                                                    >  4:1125029/1‑149 (MQ=255)                                                                     tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCa                                                                                  <  4:874079/149‑1 (MQ=255)                                                                     tttGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCa                                                                                  <  3:1012296/149‑1 (MQ=255)                                                                     tttGGTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCa                                                                                  <  3:1125029/149‑1 (MQ=255)                                                                       tGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGCGGCAGt                                                                                >  2:238395/1‑149 (MQ=255)                                                                          ttttAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCt                                                                             <  2:1145573/149‑1 (MQ=255)                                                                                         tAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATaa                                                              >  1:1013647/1‑149 (MQ=255)                                                                                          aaCGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAg                                                             <  1:970896/149‑1 (MQ=255)                                                                                                 acaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCtt                                                       <  3:1146738/148‑1 (MQ=255)                                                                                                 acaAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCtt                                                       >  4:1146738/1‑148 (MQ=255)                                                                                                   aaaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGCAAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCa                                                     >  1:899862/1‑149 (MQ=255)                                                                                                   aaaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCtgaag                                                    >  2:734833/1‑145 (MQ=255)                                                                                                    aaaaCGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGc                                                   >  1:336995/1‑149 (MQ=255)                                                                                                          aCAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCa                                             <  1:964939/149‑1 (MQ=255)                                                                                                          aCAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTGAAGCGTctta                                             <  1:734833/149‑4 (MQ=255)                                                                                                           cAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCAt                                            >  4:299314/1‑149 (MQ=255)                                                                                                            aaTATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATc                                           >  2:24326/1‑149 (MQ=255)                                                                                                              tatGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAg                                         <  4:409295/149‑1 (MQ=255)                                                                                                               atGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                                                                        >  4:918326/1‑101 (MQ=255)                                                                                                               atGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATAc                                                                                        <  3:918326/101‑1 (MQ=255)                                                                                                               atGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGCTTTTATGCCGGATGCGGTGTAAACATACGTGGCACTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAgg                                        >  4:938114/1‑149 (MQ=255)                                                                                                                  aaaaTTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGTCATGATAAGACGCTTCAGCGTCGCATCAGGCAt                                     >  1:903422/1‑149 (MQ=255)                                                                                                                  aaaaTTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGa                                                            >  3:306643/1‑126 (MQ=255)                                                                                                                  aaaaTTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGa                                                            <  4:306643/126‑1 (MQ=255)                                                                                                                  aaaaTTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCAt                                     >  1:1068838/1‑149 (MQ=255)                                                                                                                   aaaTTATTCGAGGAGTGAAAGGCAAAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCAt                                     <  2:899862/149‑1 (MQ=255)                                                                                                                   aaaTTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATa                                    >  2:1094412/1‑149 (MQ=255)                                                                                                                         ttCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGc                              <  2:1068838/149‑1 (MQ=255)                                                                                                                            gaggagTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg                           >  3:817613/1‑149 (MQ=255)                                                                                                                            gaggagTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCgg                           >  1:68304/1‑149 (MQ=255)                                                                                                                             aggagTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGa                          >  3:974308/1‑149 (MQ=255)                                                                                                                                  tGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGGGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAg                     >  2:982065/1‑149 (MQ=255)                                                                                                                                   gAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGc                    >  2:875334/1‑149 (MQ=255)                                                                                                                                     aaGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGt                  >  1:917727/1‑149 (MQ=255)                                                                                                                                     aaGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGt                  >  4:534236/1‑149 (MQ=255)                                                                                                                                      aGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAg                                         >  1:1084391/1‑125 (MQ=255)                                                                                                                                      aGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAg                                         <  2:1084391/125‑1 (MQ=255)                                                                                                                                        gC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTaaa               >  2:583418/1‑149 (MQ=255)                                                                                                                                         c‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAAc              <  2:903422/149‑1 (MQ=255)                                                                                                                                         c‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGCGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAAc              <  1:238395/149‑1 (MQ=255)                                                                                                                                           aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACg             >  4:235748/1‑149 (MQ=255)                                                                                                                                           aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACg             >  4:1037349/1‑149 (MQ=255)                                                                                                                                           aaaaaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACg             >  3:775196/1‑149 (MQ=255)                                                                                                                                              aaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCGGGCATAGGTTGCCGGATGCAGCGTAAACGCCt          >  4:947729/1‑149 (MQ=255)                                                                                                                                              aaaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCt          >  2:199098/1‑149 (MQ=255)                                                                                                                                               aaaCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCtt         <  2:68304/149‑1 (MQ=255)                                                                                                                                                taCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGc            <  4:1139069/144‑1 (MQ=255)                                                                                                                                                taCGGCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGc            >  3:1139069/2‑145 (MQ=39)                                                                                                                                                 aCGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTAt       <  1:504654/149‑1 (MQ=255)                                                                                                                                                   ggCTTCCCTATCGGGAGGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTaaa               >  3:875519/1‑139 (MQ=255)                                                                                                                                                      cTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTc  >  1:836384/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CTTTGTCGATACCACGTACATTTTCATATCCTGTCGTTTGTTTTCGATTTCAAAATATAATGAAATTATTTGTTTTTAAATATCGAGATAACGATCACAAAAACGACAATATGAAAATTATTCGAGGAGTGAAAGGC‑AAAAAAACGGCCTCCCGATAGGGAAGCCGTAGCAAAGTGCGCGTGTTTTTATGCCGGATGCGGTGTAAACATACGTGGCAGTGCTATGTAGGCATGATAAGACGCTTCAGCGTCGCATCAGGCATAGGTTGCCGGATGCAGCGTAAACGCCTTATCCGTC  >  NC_000913/2177184‑2177480

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