Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	136,457	(T)8→9	coding (8/867 nt)	speE ←	spermidine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	136,449	1	.	T	98.4%	220.2 / ‑3.6	63	coding (16/867 nt)	speE	spermidine synthase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (26/36);  total (26/37)
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.59e-01

CGCCATCCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGAC  >  NC_000913/136311‑136578                                                                                                                                            |                                                                                                                                    cGCCATCCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑‑t                                                                                                                                   >  2:123019/1‑140 (MQ=255) cGCCATCCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑‑t                                                                                                                                   >  1:272718/1‑140 (MQ=255) cGCCATCCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑‑t                                                                                                                                   >  1:687709/1‑140 (MQ=255)      tCCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑‑tttt‑tt                                                                                                                             >  2:222889/1‑140 (MQ=255)         aGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑‑tttt‑ttttt                                                                                                                          <  1:75742/140‑2 (MQ=255)         aGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑‑tttt‑ttttt                                                                                                                          <  2:552180/140‑2 (MQ=255)           cgcCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑‑tttt‑tttttcg                                                                                                                        >  2:762477/1‑137 (MQ=255)           cgcCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑‑tttt‑tttttcg                                                                                                                        >  1:402589/1‑137 (MQ=255)                ttACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCAt                                                                                                                    <  2:402589/140‑1 (MQ=255)                ttACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCAt                                                                                                                    <  1:762477/140‑1 (MQ=255)                       aCCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTT‑TTTTCGGCCATGGGTTg                                                                                                              <  2:568216/140‑1 (MQ=255)                       aCCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGa                                                                                                             <  2:20055/140‑1 (MQ=255)                              gCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATAcctcc                                                                                                      <  1:735468/140‑1 (MQ=255)                                  cGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTT‑TTTTCGGCCATGGGTTGATACCTCCttt                                                                                                   >  2:270160/1‑140 (MQ=255)                                      ctcAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTaa                                                                                              >  2:949961/1‑140 (MQ=255)                                                   aGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTaa                                                                                              >  1:399284/1‑127 (MQ=255)                                                   aGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTaa                                                                                              <  2:399284/127‑1 (MQ=255)                                                       ccTTGTGATCGGTCTTTTCATGATACAGAACGTTATCTGCCACAAAGTACTGCCCAAACTGGTCGTGTATCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAggg                                                                             <  1:795797/140‑1 (MQ=255)                                                       ccTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAggg                                                                             <  2:485946/140‑1 (MQ=255)                                                         tGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGgcg                                                                           <  1:502811/140‑1 (MQ=255)                                                         tGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGgcg                                                                           <  1:330811/140‑1 (MQ=255)                                                          ggTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGgcgc                                                                          >  2:341571/1‑140 (MQ=255)                                                           gTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCa                                                                         <  1:124437/140‑1 (MQ=255)                                                            tGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCaa                                                                        <  1:295664/140‑1 (MQ=255)                                                               tCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACat                                                                     >  1:513276/1‑140 (MQ=255)                                                               tCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACat                                                                     <  2:513276/140‑1 (MQ=255)                                                               tCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACat                                                                     <  2:209237/140‑1 (MQ=255)                                                               tCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACat                                                                     >  2:901911/1‑140 (MQ=255)                                                               tCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACat                                                                     >  1:129208/1‑140 (MQ=255)                                                               tCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACat                                                                     >  1:87734/1‑140 (MQ=255)                                                                cGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACata                                                                    >  2:539697/1‑140 (MQ=255)                                                                cGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACata                                                                    >  2:99002/1‑140 (MQ=255)                                                                cGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACa                                                                      >  2:560766/1‑138 (MQ=255)                                                                cGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACa                                                                      <  1:560766/138‑1 (MQ=255)                                                                 ggTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGgcgc                                                                          >  1:965482/1‑133 (MQ=255)                                                                 ggTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGgcgc                                                                          <  2:965482/133‑1 (MQ=255)                                                                  gTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACataat                                                                  >  2:852874/1‑140 (MQ=255)                                                                    tttttCATGATACAGAACGTTATCTTCCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAg                                                                <  2:25451/139‑1 (MQ=255)                                                                    cTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAg                                                                <  2:848696/140‑1 (MQ=255)                                                                    cTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAg                                                                <  2:836193/140‑1 (MQ=255)                                                                    cTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAg                                                                <  2:129208/140‑1 (MQ=255)                                                                          aTGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTAGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACa                                                          <  2:140023/140‑1 (MQ=255)                                                                             aTACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTg                                                       >  1:745285/1‑140 (MQ=255)                                                                               aCAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGAc                                                     <  1:341571/140‑1 (MQ=255)                                                                                 aGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACata                                                                    >  1:750738/1‑123 (MQ=255)                                                                                 aGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACcg                                                   >  2:303309/1‑140 (MQ=255)                                                                                 aGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACcg                                                   <  1:539697/140‑1 (MQ=255)                                                                                 aGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑‑TTTTCTTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACata                                                                    <  2:750738/123‑1 (MQ=255)                                                                                   aaCGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACcgcg                                                 >  1:449519/1‑140 (MQ=255)                                                                                      gTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCAGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGAt                                              >  2:511855/1‑140 (MQ=255)                                                                                       ttATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATg                                             <  1:403970/140‑1 (MQ=255)                                                                                         aTCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCa                                           <  2:745285/140‑1 (MQ=255)                                                                                         aTCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCa                                           <  1:99002/140‑1 (MQ=255)                                                                                         aTCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCa                                           <  1:123019/140‑1 (MQ=255)                                                                                              ccGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTc                                      <  2:272718/140‑1 (MQ=255)                                                                                               cGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCa                                     >  1:559831/1‑140 (MQ=255)                                                                                                gCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCaa                                    >  1:369021/1‑140 (MQ=255)                                                                                                 cAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAat                                   >  1:774699/1‑140 (MQ=255)                                                                                                   aaGTACTGCCCAAACTGGTCGGGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAatat                                 >  2:769281/1‑140 (MQ=255)                                                                                                     gTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCAGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATAttt                               <  1:511855/140‑1 (MQ=255)                                                                                                         tGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTAc                           >  1:685164/1‑140 (MQ=255)                                                                                                           cccAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACaa                         >  1:611299/1‑140 (MQ=255)                                                                                                             cAAACTGGTCGTGTAGCGTTTCATGCCACTGTTTTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACaaa                        <  1:270160/140‑1 (MQ=255)                                                                                                             cAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACaaaa                       <  2:685164/140‑1 (MQ=255)                                                                                                             cAAACTGGTCGTGTAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACaaaa                       <  1:769281/140‑1 (MQ=255)                                                                                                              aaaCTGGTCGTGTAGCGTTTCATGCCACTGTTTTTT‑TTTTCGGCCCTGGGTTGATACCTCCTTTGTTAACACCCTTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCCCGGATGCACGGTCAATATTTCTACaaaa                       >  2:952856/1‑140 (MQ=255)                                                                                                                         gtAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGa           <  2:369021/140‑1 (MQ=255)                                                                                                                         gtAGCGTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGa           <  2:774699/140‑1 (MQ=255)                                                                                                                              gTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGattatt      <  1:101541/140‑1 (MQ=255)                                                                                                                              gTTTCATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGattatt      <  2:611299/140‑1 (MQ=255)                                                                                                                                  cATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGAc  >  1:42883/1‑140 (MQ=255)                                                                                                                                  cATGCCACTGT‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGAc  >  1:77360/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                    CGCCATCCAGCGCCATTACGCGACCAAATGCAGCGTTCTCAAAAATGATCAGATCCTGGTGATCGGTCTTTTCATGATACAGAACGTTATCTACCGCAAAGTACTGCCCAAACTGGTCGTGTAGCGTTTCATGCCACTG‑‑TTTT‑TTTTCGGCCATGGGTTGATACCTCCTTTGTTAACACCCGTAAAAAAAGGGCGCAACATAATAGCTAACATTGACCGCGGATGCACGGTCAATATTTCTACAAAAAGGTATCAGGGATTATTTGAC  >  NC_000913/136311‑136578

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: ‑