Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	346,022	G→A	7.0%	R106C (CGC→TGC)	yahN ←	amino acid exporter for proline, lysine, glutamate, homoserine

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	346,022	0	G	A	7.0%	171.9 / 4.1	57	R106C (CGC→TGC)	yahN	amino acid exporter for proline, lysine, glutamate, homoserine
Reads supporting (aligned to +/- strand):  ref base G (30/23);  new base A (2/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.78e-02

CCAGATAATTGATGCGAGCACAATCCCCGCCCAGGCCATTAAACGTGCCCATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTT‑GTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAA  >  NC_000913/345810‑346289                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                             ccAGATAATTGATGCGAGCACAATCCCCGCCCAGGCCATTAAACGTGCCCATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGaaaaaa                                                                                                                                                                                                                <  2:413084/275‑1 (MQ=255)                            cGCCCAGGCCATTAAACGTGCCCATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAACAATACTCTCACCATGCGTAATTAGCGTTg                                                                                                                                                                                     >  2:306964/1‑275 (MQ=255)                                                 ccATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTaa                                                                                                                                                                                              >  2:224257/1‑245 (MQ=255)                                                        ggTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACAGACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGAATCACGCAAACCATAAGAGATAAGCGCAGACGACGATTCTGATAAGCGaa                                                                                                                                                                                                                    >  4:145782/1‑216 (MQ=255)                                                            tttCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTTTTCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGaaaaaa                                                                                                                                                                                                                <  2:145783/216‑1 (MQ=255)                                                                ggCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGTCGGCGCATGCTGCACCACGCAAACCATACGAGATAAGCGCTGCCGACGATTCNGACAAGCGAAAACACCCCCACACACCGCGTAATTCTCGTTGCAAGCCCAAACAAACCCAACCCGCAACAAAATGcc                                                                                                                                                  >  3:109689/1‑273 (MQ=255)                                                                    tttAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCAt                                                                                                                                                                                                                                                                           <  7:63623/149‑1 (MQ=255)                                                                      tAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCAtgc                                                                                                                                                                                                                                                                         <  5:23009/149‑1 (MQ=255)                                                                          gTTACTGAGAAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGca                                                                                                                                                                                                                                                                     >  7:48462/1‑149 (MQ=255)                                                                                 agaAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCa                                                                                                                                                                                                                                                              >  7:215797/1‑149 (MQ=255)                                                                                 agaAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGcc                                                                                                                                 >  2:320126/1‑274 (MQ=255)                                                                                   aaaaTACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCaaa                                                                                                                                                                                                                                                            >  7:453306/1‑149 (MQ=255)                                                                                         cTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATaa                                                                                                                                                                                                                                                      <  6:211364/149‑1 (MQ=255)                                                                                           gATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAaga                                                                                                                                                                                                                                                    <  5:215663/149‑1 (MQ=255)                                                                                             tAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTCGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAANCCATACGAGATAAGCGCCGTCGACGATTCTGTTAAGCGACAAAATccccccc                                                                                                                                                                                                        >  3:132868/1‑188 (MQ=255)                                                                                                aaaaTAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATaa                                                                                                                                                                                                                                               <  6:437172/149‑1 (MQ=255)                                                                                                   aTAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAgcg                                                                                                                                                                                                                                            >  8:214601/1‑149 (MQ=255)                                                                                                     aaaaCGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCAtgc                                                                                                                                                                                                                                                                         >  5:188702/1‑118 (MQ=255)                                                                                                     aaaaCGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCAtgc                                                                                                                                                                                                                                                                         <  6:188702/118‑1 (MQ=255)                                                                                                      aaaCGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCg                                                                                                                                                                                        >  4:327058/1‑198 (MQ=255)                                                                                                        aCGGTTTGCGGGGTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCgccgcc                                                                                                                                                                                                                                       >  8:373928/1‑149 (MQ=255)                                                                                                          ggTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGc                                                                                                                                                                                    <  2:327062/198‑1 (MQ=255)                                                                                                               gCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGAtt                                                                                                                                                                                                                                <  8:48462/149‑1 (MQ=255)                                                                                                                cGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTc                                                                                                                                                                                                                               <  7:214601/149‑1 (MQ=255)                                                                                                                  ggTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTg                                                                                                                                                                                                                             <  8:453306/149‑1 (MQ=255)                                                                                                                    ttAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGAt                                                                                                                                                                                                                           <  7:174582/149‑1 (MQ=255)                                                                                                                       gagagaTCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCACATGCTGCACCACGCAAACCATaa                                                                                                                                                                                                                                                      >  7:207193/1‑119 (MQ=255)                                                                                                                       gagagaTCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCACATGCTGCACCACGCAAACCATaa                                                                                                                                                                                                                                                      <  8:207193/119‑1 (MQ=255)                                                                                                                       gagagaTCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCACATGCTGCACCACGCAAACCATaa                                                                                                                                                                                                                                                      <  6:279175/119‑1 (MQ=255)                                                                                                                       gagagaTCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCACATGCTGCACCACGCAAACCATaa                                                                                                                                                                                                                                                      >  5:279175/1‑119 (MQ=255)                                                                                                                                gTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAAt                                                                                                                                                                                                               >  5:156116/1‑149 (MQ=255)                                                                                                                                   attaatCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTTACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGcgccg                                                                                                                                                                                                                                         >  2:491103/1‑120 (MQ=255)                                                                                                                                       atCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTTACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCcga                                                                                                                                                                                                                                     <  4:491098/120‑1 (MQ=255)                                                                                                                                           gcggcgAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCCCGCAAACCATAAGAGATAAGCGCCGCCGACGCTTTCGATAAGCGAACAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAACCAAACCCAACCCGGAATAAAATGCACCGCCCAGCGCCCCGCCCAGACCGGTCAGCACCCCTGCGCGTCCCACGGAAGCCAGGCTGTGTT‑Gtaccgggcgc                                                                      >  4:149120/1‑267 (MQ=255)                                                                                                                                               cgAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAgg                                                                                         <  2:149121/252‑1 (MQ=255)                                                                                                                                                 aaaaaaGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTaa                                                                                                                                                                                              >  8:187420/1‑149 (MQ=255)                                                                                                                                                   aaaaGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAAtt                                                                                                                                                                                            >  5:306022/1‑149 (MQ=255)                                                                                                                                                   aaaaGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAAtt                                                                                                                                                                                            >  7:365927/1‑149 (MQ=255)                                                                                                                                                      aGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCt                                                                                                                                                                                                                              >  8:394076/1‑112 (MQ=255)                                                                                                                                                      aGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCt                                                                                                                                                                                                                              <  7:394076/112‑1 (MQ=255)                                                                                                                                                         cATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGtt                                                                                                                                                                                      <  8:215797/149‑1 (MQ=255)                                                                                                                                                                     cgcTAATCGTTTGTTGTAGTGTGCTCATTGGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGcc                                                                                                                                 <  1:143180/190‑1 (MQ=255)                                                                                                                                                                     cgcTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTT‑GTACTACCACAAAGAGATTGGCTCCCGGATTaaaaaaa                                           >  2:473789/1‑275 (MQ=255)                                                                                                                                                                           tCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAaaccc                                                                                                                                                                    <  5:95147/149‑1 (MQ=255)                                                                                                                                                                              gttgttGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAaacccaac                                                                                                                                                                 >  6:317409/1‑149 (MQ=255)                                                                                                                                                                               ttgttgTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAaacccaacc                                                                                                                                                                >  6:441223/1‑149 (MQ=255)                                                                                                                                                                                 gttgtGGTGGGGTCATTTGCGGGGTTTACTGGCAGAGGAGGCAGAACCACGGAAAACAGAAGAGAGAAGGGGCGCCGACGATTCTTATAATCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGAGCAAACAAACCCAACCCGGAATAAAATGCTGCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTT‑GTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTATTCACGaa                                <  4:320122/274‑1 (MQ=255)                                                                                                                                                                                         gtgCTCATTTGCGGTGTTGACTGGCGGCGCATTCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCCCCTCAAACTGCGTAAGTAGCGTTGGAATACAAAACAAAACCAACCCAGAATAACAAGCATCGCCCAGAGCCACGCACAGCGCAGTCAGGAACCCGCAGGCGCGAAGGGAAGCACGGCGGGTTTGGTAC‑ACCACAAAGAGAATAGGACCAGTGTAAAAAAAAGTATACACGAAAAgacaaac                       >  4:329394/1‑269 (MQ=255)                                                                                                                                                                                         gtgCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCCGGTTT‑GTACTACCACAAAGAGCTTGGCTCCCGGATAAAAAAAAGTAATCACGAACAGTCCTAc                       >  2:307830/1‑275 (MQ=255)                                                                                                                                                                                          tgCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAAt                                                                                                                                                     >  5:117614/1‑149 (MQ=255)                                                                                                                                                                                            cTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCCGCGCCACGCCCAGCCCGGTCAACACCCCTGCGCCTCCACCGGAAGCCAGGCCGGTTT‑GTCCTACCACAAAGAGATTGGCCCCAGGCTAATACAAAAGACTCACGAAAAGACCTAcag                     >  4:60453/1‑272 (MQ=255)                                                                                                                                                                                            cTCATTTGCGGTGTTGACTGGCGGCGCATGCTGAACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATGAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACTCGGAATAAAATGCAACCCCCAGCGCTACACCCAGCCCGGTCCGCACCACTTCGCGTCGACCGGAAGCCAGGCTCGTTT‑GTACTACCCTAAAGCGATTGGCTCCCGGATTAAAAAACATAAGCACGCACAGCCCTACggc                    >  4:329701/1‑274 (MQ=255)                                                                                                                                                                                                tttGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATGAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTT‑GTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGaa                                <  2:329705/259‑1 (MQ=255)                                                                                                                                                                                                         ttGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGGCCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGcc                                                                                                                                      <  7:373928/149‑1 (MQ=255)                                                                                                                                                                                                           gACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCAc                                                                                                                                    >  6:194065/1‑149 (MQ=255)                                                                                                                                                                                                           gACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTT‑GTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATg     >  2:2902/1‑275 (MQ=255)                                                                                                                                                                                                              tGGCGGCGCATGCTGCACCACGAAAACCATAAGGGATAATCGCCGCCGACGATGCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTT‑GTACGACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCaa  <  1:109689/275‑1 (MQ=255)                                                                                                                                                                                                                   gcgcATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGcc                                                                                                                            >  8:355923/1‑149 (MQ=255)                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                             CCAGATAATTGATGCGAGCACAATCCCCGCCCAGGCCATTAAACGTGCCCATGTTGGTGTTTCGGCATTTAATGTTACTGAGAAAATACTGATAAAAAATAAAACGGTTTGCGGGTTAGAGAGATCGGTAATTAATCCGCGGCGAAAAAAGACATACCAGGGGGCGCTAATCGGTTGTTGTAGTGTGCTCATTTGCGGTGTTGACTGGCGGCGCATGCTGCACCACGCAAACCATAAGAGATAAGCGCCGCCGACGATTCTGATAAGCGAAAAAATCTCCTCACACTGCGTAATTAGCGTTGCAAGACCAAACAAACCCAACCCGGAATAAAATGCATCGCCCAGCGCCACGCCCAGCCCGGTCAGCACCCCTGCGCGTCGACCGGAAGCCAGGCTGGTTT‑GTACTACCACAAAGAGATTGGCTCCCGGATTAAAAAAAGTAATCACGAACAGTCCTACGGTCAGGTAAACGGCATGCAA  >  NC_000913/345810‑346289

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