Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,366,855	C→A	8.3%	intergenic (‑82/+36)	dcuA ← / ← aspA	C4‑dicarboxylate antiporter/aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,366,855	0	C	A	8.3%	199.7 / 7.6	60	intergenic (‑82/+36)	dcuA/aspA	C4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base A (2/3);  total (29/31)
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

CTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAA  >  NC_000913/4366707‑4367003                                                                                                                                                      |                                                                                                                                                     cTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATc                                                                                                                                                      >  4:91362/1‑149 (MQ=255) cTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATc                                                                                                                                                      >  1:215393/1‑149 (MQ=255)  tATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCt                                                                                                                                                     >  1:431083/1‑149 (MQ=255)   aTTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTg                                                                                                                                                    >  1:295961/1‑149 (MQ=255)      ccccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACg                                                                                                                                                 >  4:41214/1‑149 (MQ=255)       cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGt                                                                                                                                                >  3:522231/1‑149 (MQ=255)       cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGt                                                                                                                                                <  3:36676/149‑1 (MQ=255)       cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGt                                                                                                                                                >  3:168592/1‑149 (MQ=255)       cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGt                                                                                                                                                >  1:178005/1‑149 (MQ=255)       cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGt                                                                                                                                                <  1:458832/149‑1 (MQ=255)        ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                               >  2:495475/1‑149 (MQ=255)        ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                               >  4:144099/1‑149 (MQ=255)         cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGc                                                                                                                                              <  4:176500/149‑1 (MQ=255)          ccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGcc                                                                                                                                             <  2:248554/149‑1 (MQ=255)          ccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATAGTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGcc                                                                                                                                             >  2:349089/1‑149 (MQ=255)           cAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAATAAAAAAGGCACGTCAGATGACGTGcc                                                                                                                                             <  3:412150/149‑1 (MQ=255)            aaTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTCTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCtt                                                                                                                                           <  2:456588/149‑1 (MQ=255)                    gcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAttt                                                                                                                                   <  3:416114/149‑1 (MQ=255)                    gcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAttt                                                                                                                                   >  2:518788/1‑149 (MQ=255)                        aagaagATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCt                                                                                                                                            >  1:468128/1‑136 (MQ=255)                        aagaagATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCAGATGACGTGCCt                                                                                                                                            <  2:468128/136‑1 (MQ=255)                            agaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAcc                                                                                                                           <  1:495475/149‑1 (MQ=255)                            agaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAcc                                                                                                                           <  3:144099/149‑1 (MQ=255)                                 gCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTa                                                                                                                      <  3:127882/149‑1 (MQ=255)                                    aGCAAATCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGa                                                                                                                   <  1:373290/149‑1 (MQ=255)                                    aGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGa                                                                                                                   >  4:366943/1‑149 (MQ=255)                                     gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAt                                                                                                                  <  3:366943/149‑1 (MQ=255)                                                 tgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCtt                                                                                                      <  2:351933/149‑1 (MQ=255)                                                  gaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCttt                                                                                                     >  1:302962/1‑149 (MQ=255)                                                  gaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCttt                                                                                                     <  4:522231/149‑1 (MQ=255)                                                    gTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                                   >  4:150227/1‑149 (MQ=255)                                                    gTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                                   <  3:150227/149‑1 (MQ=255)                                                              tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGc                                                                                         >  1:492305/1‑149 (MQ=255)                                                                  atatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGttt                                                                                     >  4:468511/1‑149 (MQ=255)                                                                       aGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCtt                                                                                >  2:231828/1‑149 (MQ=255)                                                                       aGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCtt                                                                                <  1:349089/149‑1 (MQ=255)                                                                         ccttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTg                                                                              <  2:178005/149‑1 (MQ=255)                                                                          cttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGt                                                                             <  3:114493/149‑1 (MQ=255)                                                                            tccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTaa                                                                           <  1:231828/149‑1 (MQ=255)                                                                             ccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAg                                                                          <  2:215393/149‑1 (MQ=255)                                                                               ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGcc                                                                        >  1:225185/1‑149 (MQ=255)                                                                                   tttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGGTTTGTAAGCCGGGt                                                                    <  3:41214/149‑1 (MQ=255)                                                                                             gTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTGAGCCGGGTGCATCAGAtt                                                          >  3:138007/1‑149 (MQ=255)                                                                                             gTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGAtt                                                          >  4:450265/1‑149 (MQ=255)                                                                                              ttGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTc                                                         >  1:64732/1‑149 (MQ=255)                                                                                              ttGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTc                                                         <  4:168592/149‑1 (MQ=255)                                                                                                gATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAATGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTg                                                       >  4:141083/1‑149 (MQ=255)                                                                                                  tatTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTa                                                     <  1:448373/149‑1 (MQ=255)                                                                                                  tatTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTa                                                     <  2:393975/149‑1 (MQ=255)                                                                                                        aTTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGgaga                                               >  2:575018/1‑149 (MQ=255)                                                                                                             tATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAata                                          <  2:302962/149‑1 (MQ=255)                                                                                                                  ttAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTc                                     <  1:272783/149‑1 (MQ=255)                                                                                                                      gTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGt                                 >  4:476181/1‑149 (MQ=255)                                                                                                                               cac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCa                        <  1:575018/149‑1 (MQ=255)                                                                                                                               cac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCa                        <  3:450265/149‑1 (MQ=255)                                                                                                                                ac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAata                                          >  3:435423/1‑130 (MQ=255)                                                                                                                                ac‑aAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAata                                          <  4:435423/130‑1 (MQ=255)                                                                                                                                   aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                     >  1:556760/1‑149 (MQ=255)                                                                                                                                         aaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGc                            >  4:254815/1‑136 (MQ=255)                                                                                                                                         aaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGc                            <  3:254815/136‑1 (MQ=255)                                                                                                                                            ggCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcg            >  1:366762/1‑149 (MQ=255)                                                                                                                                             gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCGTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc           >  2:439936/1‑149 (MQ=255)                                                                                                                                                      cTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGaa  >  1:184817/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAC‑AAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAA  >  NC_000913/4366707‑4367003

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