Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,243,989	T→C	5.0%	E251E (GAA→GAG)	uxaC ←	uronate isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,243,989	0	T	C	5.0%	262.6 / 2.9	82	E251E (GAA→GAG)	uxaC	uronate isomerase
Reads supporting (aligned to +/- strand):  ref base T (43/34);  new base C (2/2);  total (45/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.01e-01

GCGCCAATGTGGTACTGCTGTACCCAGCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTC  >  NC_000913/3243843‑3244136                                                                                                                                                   |                                                                                                                                                    gcgcCAATGTGGTACTGCTGTACCCAGCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCttc                                                                                                                                                   <  8:73606/149‑1 (MQ=255)    ccAATGTGGTACTGCTGTACCCAGCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCttcgtt                                                                                                                                                <  3:60503/149‑1 (MQ=255)     cAATGTGGTACTGCTGTACCCAGCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCttcgttc                                                                                                                                               >  5:287991/1‑149 (MQ=255)         gtgGTACTGCTGTACCCAGCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCtt                                                                                                                                           >  7:41898/1‑149 (MQ=255)               ctgctgTACCCAGCCGCGACGTGCAAAATCGGCACCAAGGAACACCAGTACCGCAGTAATGCACTGCGCCACTTCGAGCTCCCTCAGGGTTTCGCCACCCAGACGTCGCGCCAGAATGCTCTCCAGTACCGCTTCGTTCGCTTCAGCa                                                                                                                                      >  3:24104/1‑148 (MQ=255)                tgctgTACCCAGCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCttcgtt                                                                                                                                                <  6:84929/137‑1 (MQ=255)                tgctgTACCCAGCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCttcgtt                                                                                                                                                >  5:84929/1‑137 (MQ=255)                tgctgTACCCAGCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCaaa                                                                                                                                    >  1:10937/1‑149 (MQ=255)                tgctgTACCCAGCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCaaa                                                                                                                                    >  2:136004/1‑149 (MQ=255)                     tACCCAGCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATca                                                                                                                               >  3:375030/1‑149 (MQ=255)                       cccAGCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATcacc                                                                                                                             <  2:412535/149‑1 (MQ=255)                        ccAGCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATcacca                                                                                                                            >  7:212303/1‑149 (MQ=255)                        ccAGCCGCGAAGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATcaac                                                                                                                             >  8:191624/1‑146 (MQ=255)                           gCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACAt                                                                                                                         >  4:235358/1‑149 (MQ=255)                                   gTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGcgcg                                                                                                                 >  7:401366/1‑149 (MQ=255)                                    tGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGt                                                                                                                >  1:267941/1‑149 (MQ=255)                                     gCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTg                                                                                                               >  6:344443/1‑149 (MQ=255)                                      cATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTgg                                                                                                              <  7:302500/149‑1 (MQ=255)                                      cATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTg                                                                                                               <  8:212303/148‑1 (MQ=255)                                       atatTCGGCACCAAGGAACACCTGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGt                                                                                                             <  2:10937/149‑1 (MQ=255)                                        tatTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGt                                                                                                             >  5:349717/1‑148 (MQ=255)                                        tatTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGt                                                                                                             >  1:289948/1‑148 (MQ=255)                                         atTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTc                                                                                                                                          >  2:397181/1‑118 (MQ=255)                                         atTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTc                                                                                                                                          <  1:397181/118‑1 (MQ=255)                                          ttCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGt                                                                                                             >  7:283187/1‑146 (MQ=255)                                          ttCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCCTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAg                                                                                                          >  7:398071/1‑149 (MQ=255)                                              gCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGacac                                                                                                      >  7:300413/1‑149 (MQ=255)                                                aCCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACtt                                                                                                    >  1:228860/1‑149 (MQ=255)                                                aCCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACtt                                                                                                    >  1:142781/1‑149 (MQ=255)                                                 cctaGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCttcgttc                                                                                                                                               >  8:363110/4‑105 (MQ=255)                                                 ccAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACtt                                                                                                    >  1:28180/1‑148 (MQ=255)                                                  cAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTa                                                                                                  <  4:375030/149‑1 (MQ=255)                                                  cAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGAAACTTTa                                                                                                  <  6:287991/149‑1 (MQ=255)                                                  cAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTAATGGTTTCGCCAGCCAGACGGCGCGCCAGAATGATGTCCAGTTCAGCTTCGTTCGCTTCAGCAAAAATCAACAAATACAGCGCGTGGGCAGACACttt                                                                                                   >  2:308068/1‑148 (MQ=255)                                                    aGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTa                                                                                                  >  6:28169/1‑147 (MQ=255)                                                       aaCACCCGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGACAGAATGCTTTCCAGTTCCGCTTGGTTCGCTTCAGCGAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGcc                                                                                             <  6:265148/149‑1 (MQ=255)                                                         cacCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGt                                                                                           >  6:312185/1‑149 (MQ=255)                                                             aGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAgc                                                                                       <  2:289948/149‑1 (MQ=255)                                                               tACCGCCGTGTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAgcgg                                                                                     <  4:29451/149‑1 (MQ=255)                                                               tACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAgcgg                                                                                     <  5:28169/149‑1 (MQ=255)                                                                 ccGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAgcggc                                                                                    >  4:259915/1‑148 (MQ=255)                                                                  cGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCttcg                                                                                                                                                  <  8:375531/85‑1 (MQ=255)                                                                  cGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCttcg                                                                                                                                                  >  7:375531/1‑85 (MQ=255)                                                                  cGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAgc                                                                                                                    <  1:135668/115‑1 (MQ=255)                                                                  cGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAgc                                                                                                                    >  2:135668/1‑115 (MQ=255)                                                                   gCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGaa                                                                                 <  8:283187/149‑1 (MQ=255)                                                                   gCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGaa                                                                                 <  2:267941/149‑1 (MQ=255)                                                                   gCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGaa                                                                                 <  6:215990/149‑1 (MQ=255)                                                                    cAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCa                                                                                                           >  1:140685/1‑122 (MQ=255)                                                                    cAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCa                                                                                                           <  2:140685/122‑1 (MQ=255)                                                                       ttttGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGa                                                                             >  1:348082/1‑149 (MQ=255)                                                                           gAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCCTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATc                                                                           >  8:78193/1‑147 (MQ=255)                                                                            aaCTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCa                                                                         <  3:235358/148‑1 (MQ=255)                                                                            aaCTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAg                                                                        <  8:401366/149‑1 (MQ=255)                                                                                  gcCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGttt                                                                  <  5:287256/149‑1 (MQ=255)                                                                                   ccACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACa                                                                                                >  7:169988/1‑118 (MQ=255)                                                                                   ccACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACa                                                                                                <  8:169988/118‑1 (MQ=255)                                                                                   ccACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGttt                                                                  >  2:246885/1‑148 (MQ=255)                                                                                   ccACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGttt                                                                  <  1:246885/148‑1 (MQ=255)                                                                                    cACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTgg                                                                >  4:250841/1‑149 (MQ=255)                                                                                    cACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCAACACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTgg                                                                <  7:191624/149‑1 (MQ=255)                                                                                         cGTGCTCGCTCAGGTTTTCTCCAGCGAGTCGGCTCGCCAGAATTCTTTCCAGTTCAGCTTCGTTCGATTCAGCAAACATCACCACATCAAGCGCGTGGTCAGACACTTTAGAGACGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAg                                                            <  2:348082/148‑1 (MQ=255)                                                                                         cGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCa                                                             <  2:28180/147‑1 (MQ=255)                                                                                                 ctcAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCCTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTg                                                   <  7:78193/149‑1 (MQ=255)                                                                                                  tcAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGc                                                  <  3:259915/149‑1 (MQ=255)                                                                                                    aGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGc                                                  <  5:344443/147‑1 (MQ=255)                                                                                                    aGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCa                                                 <  2:228860/148‑1 (MQ=255)                                                                                                    aGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAg                                                <  6:349717/149‑1 (MQ=255)                                                                                                      ggTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGGGAGCGGCGAAGTGCTCCAGACGTTTGGTCAGGGCAGTTTGCAGGt                                              >  1:129265/1‑149 (MQ=255)                                                                                                        tttCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCa                                                 >  7:106138/1‑144 (MQ=255)                                                                                                        tttCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCa                                                 <  8:106138/144‑1 (MQ=255)                                                                                                        tttCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGt                                              >  7:129485/1‑147 (MQ=255)                                                                                                            gccagccaGACGGCGCGCCAGAATTCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCa                                         >  6:236181/1‑148 (MQ=255)                                                                                                                 ccagACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGATGTGATCCAGACGtt                                                                   >  2:184974/1‑117 (MQ=255)                                                                                                                 ccagACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGATGTGATCCAGACGtt                                                                   <  1:184974/117‑1 (MQ=255)                                                                                                                          gcgcCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCg                          <  5:312185/149‑1 (MQ=255)                                                                                                                              cAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGCGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTAt                      >  5:211551/1‑149 (MQ=255)                                                                                                                              cAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTAt                      >  8:42007/1‑149 (MQ=255)                                                                                                                               aGAATGCTGTCCAGTTCCGCCTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATc                     <  8:398071/149‑1 (MQ=255)                                                                                                                                     cTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGaaa                <  5:236181/148‑1 (MQ=255)                                                                                                                                     cTGGCCAGTTACGCATCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTAGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGGATCGGAAAc               <  1:308068/149‑1 (MQ=255)                                                                                                                                        tCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTc            >  5:203134/1‑149 (MQ=255)                                                                                                                                        tCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTc            >  4:74541/1‑149 (MQ=255)                                                                                                                                        tCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCAGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTc            >  3:175454/1‑149 (MQ=255)                                                                                                                                         ccAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCg           >  8:143149/1‑149 (MQ=255)                                                                                                                                            gTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGccc        <  3:250841/149‑1 (MQ=255)                                                                                                                                            gTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGcc         >  5:145385/1‑148 (MQ=255)                                                                                                                                                cGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCGGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGc     <  8:300413/148‑1 (MQ=255)                                                                                                                                                cGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCt    >  1:119279/1‑149 (MQ=255)                                                                                                                                                  cTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTATACAGCCCTGAGCGGCGAAGTGATCCAGACCTATGGTCAGGGCAGATAGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTc  >  3:137913/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GCGCCAATGTGGTACTGCTGTACCCAGCCGCGACGTGCATATTCGGCACCAAGGAACACCAGTACCGCAGTTTTGAACTGCGCCACTTCGTGCTCGCTCAGGGTTTCGCCAGCCAGACGGCGCGCCAGAATGCTGTCCAGTTCCGCTTCGTTCGCTTCAGCAAACATCACCACATCCAGCGCGTGGTCAGACACTTTACAGCCGTGAGCGGCGAAGTGATCCAGACGTTTGGTCAGGGCAGTTTGCAGGTCAGCAAAGCGGCGAATGTCGGTATCGGAAACTTCGCCCAGCTTC  >  NC_000913/3243843‑3244136

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