Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,001,626	G→A	6.4%	G428S (GGC→AGC)	xdhA →	xanthine dehydrogenase, molybdenum binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,001,626	0	G	A	6.4%	203.8 / 3.5	62	G428S (GGC→AGC)	xdhA	xanthine dehydrogenase, molybdenum binding subunit
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base A (2/2);  total (28/34)
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.46e-01

GTTTACGCAACGCCGCCCGCGAAGGAGATGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAGCGGCGCGACGGAAATCGGTCAGG  >  NC_000913/3001483‑3001773                                                                                                                                                |                                                                                                                                                    gTTTACGCCACGCCGCCCGCGAAGGAGATGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAAt                                                                                                                                                <  1:614535/149‑1 (MQ=255) gTTTACGCAACGCCGCCCGCGAAGGAGATGCTAATCCGCTCACGGGCAAATGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAAt                                                                                                                                                >  4:416788/1‑149 (MQ=255)  tttACGCAACGCCGCCCGCGAAGGAGATGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGCATGGGAAAAACGCCGTGCAGGATGCCAGAACCAGCAAGGCAAtt                                                                                                                                               >  4:543909/1‑149 (MQ=255)  tttACGCAACGCCGCCCGCGAAGGAGATGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAAtt                                                                                                                                               <  3:90295/149‑1 (MQ=255)  tttACGCAACGCCGCCCGCGAAGGAGATGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAAtt                                                                                                                                               <  3:171632/149‑1 (MQ=255)       gCAACGCCGCCCGCGCAGGAGATGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCTGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTtcgc                                                                                                                                          >  3:569494/1‑149 (MQ=255)             ccgcccgcGAAGGAGATGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCgcggc                                                                                                                                    <  3:537643/149‑1 (MQ=255)                cccgcGAAGGAGATGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGtt                                                                                                                                 <  2:157774/149‑1 (MQ=255)                     gAAGGAGATGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGt                                                                                                                            >  1:523896/1‑149 (MQ=255)                     gAAGGAGATGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGt                                                                                                                            <  3:246037/149‑1 (MQ=255)                     gAAGGAGATGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGt                                                                                                                            <  3:532691/149‑1 (MQ=255)                        ggAGATGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGc                                                                                                                         <  2:201651/149‑1 (MQ=255)                         gagaTGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGcc                                                                                                                        <  4:569494/149‑1 (MQ=255)                           gaTGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTg                                                                                                                      >  2:172166/1‑149 (MQ=255)                              gCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAAGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGttt                                                                                                                   >  2:72394/1‑149 (MQ=255)                              gCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAAGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGttt                                                                                                                   >  2:26180/1‑149 (MQ=255)                                        tCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACAcc                                                                                                         >  1:193571/1‑149 (MQ=255)                                           cGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAAGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACctct                                                                                                      <  1:72394/149‑1 (MQ=255)                                           cGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAAGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACctct                                                                                                      <  1:26180/149‑1 (MQ=255)                                            gggCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTa                                                                                                     >  3:2287/1‑149 (MQ=255)                                            gggCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTa                                                                                                     >  3:386366/1‑149 (MQ=255)                                             ggCAAACGTATTTACAGCGCAGGGCTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGCGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTaa                                                                                                    <  3:408770/149‑1 (MQ=255)                                                aaaCGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAacac                                                                                                 >  3:175097/1‑149 (MQ=255)                                                  aCGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCt                                                                                               >  3:521810/1‑149 (MQ=255)                                                  aCGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCt                                                                                               <  3:439573/149‑1 (MQ=255)                                                       tttACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCt                                                                                          >  3:347209/1‑149 (MQ=255)                                                          aCAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTc                                                                                       >  3:602030/1‑149 (MQ=255)                                                           cAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCg                                                                                      <  2:448339/149‑1 (MQ=255)                                                             gcgcAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGc                                                                                    <  4:347209/149‑1 (MQ=255)                                                                        ccGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGc                                                                         >  3:90287/1‑149 (MQ=255)                                                                        ccGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGc                                                                         <  4:560089/149‑1 (MQ=255)                                                                         cGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCa                                                                        <  2:202264/149‑1 (MQ=255)                                                                          ggAGTGTCTTGGAAAAGGCCGGAAAATCTTAGAACGGGAAAAACGCCGTGCAGAATTCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTACCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAg                                                                       <  4:518654/149‑1 (MQ=255)                                                                                      aaaaGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTc                                                           <  2:522514/149‑1 (MQ=255)                                                                                      aaaaGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTc                                                           >  2:501915/1‑149 (MQ=255)                                                                                       aaaGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCt                                                          <  4:90287/149‑1 (MQ=255)                                                                                       aaaGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCt                                                          >  3:380868/1‑149 (MQ=255)                                                                                         aGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCtga                                                        <  2:523896/149‑1 (MQ=255)                                                                                         aGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTAGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCCACCACCTGGCCTGTTGGAGTAGAAATAGCAGGCGCGCGCCTTCtga                                                        <  4:380831/149‑1 (MQ=255)                                                                                           gCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCtgatg                                                      <  1:562689/149‑1 (MQ=255)                                                                                              ggAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAAt                                                   <  3:273616/149‑1 (MQ=255)                                                                                                aaaaTCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGGACCAGCAAGGCAATTTACGCCGTGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGTCTTCTAATGAATCa                                                 >  3:383359/1‑149 (MQ=255)                                                                                                aaaaTCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGTCTTCTGATGAATCa                                                 <  1:54334/149‑1 (MQ=255)                                                                                                    tCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGAt                                             >  3:417142/1‑149 (MQ=255)                                                                                                    tCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGAt                                             >  4:267056/1‑149 (MQ=255)                                                                                                         gAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAAc                                        >  1:510214/1‑149 (MQ=255)                                                                                                                    aCGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGc                             <  2:193571/149‑1 (MQ=255)                                                                                                                       ccGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCaaa                          >  4:170369/1‑149 (MQ=255)                                                                                                                         gTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAgc                        >  1:209151/1‑149 (MQ=255)                                                                                                                            cAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAgcggc                     <  4:211843/149‑1 (MQ=255)                                                                                                                            cAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAgcg                       >  2:232908/1‑147 (MQ=255)                                                                                                                            cAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAgcg                       <  1:232908/147‑1 (MQ=255)                                                                                                                                aTGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAGCGGCGCGa                 <  4:521810/149‑1 (MQ=255)                                                                                                                                 tGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAGCGGCGCGAc                <  3:162327/149‑1 (MQ=255)                                                                                                                                      gAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAGCGGCGCGACGGaaa           >  4:460148/1‑149 (MQ=255)                                                                                                                                         ccAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAGCGGCGCGACGGAAATCg        <  3:460148/149‑1 (MQ=255)                                                                                                                                           aGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAGCGGCGCGACGGAAATCGGt      >  2:310695/1‑149 (MQ=255)                                                                                                                                            gCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAGCGGCGCGACGGAAATCGGTc     <  1:130107/149‑1 (MQ=255)                                                                                                                                            gCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAGCGGCGCGACGGAAATCGGTc     <  3:416788/149‑1 (MQ=255)                                                                                                                                              aaGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAGCGGCGCGACGGAAATCGGTCAg   >  2:73361/1‑149 (MQ=255)                                                                                                                                               aGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAGCGGCGCGACGGAAATCGGTCAgg  <  4:175097/149‑1 (MQ=255)                                                                                                                                               aGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAGCGGCGCGACGGAAATCGGTCAgg  >  2:597962/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    GTTTACGCAACGCCGCCCGCGAAGGAGATGCTAATCCGCTCACGGGCAAACGTATTTACAGCGCAGGGTTGCCGGAGTGTCTTGAAAAAGGCCGGAAAATCTTTGAATGGGAAAAACGCCGTGCAGAATGCCAGAACCAGCAAGGCAATTTGCGCCGCGGCGTTGGCGTCGCCTGTTTTAGCTACACCTCTAACACCTGGCCTGTCGGCGTAGAAATAGCAGGCGCGCGCCTTCTGATGAATCAGGATGGAACCATCAACGTGCAAAGCGGCGCGACGGAAATCGGTCAGG  >  NC_000913/3001483‑3001773

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