Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,866,469	T→C	6.9%	*143* (TAA→TAG)	ibpB ←	heat shock chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,866,469	0	T	C	6.9%	222.5 / 6.0	72	*143* (TAA→TAG)	ibpB	heat shock chaperone
Reads supporting (aligned to +/- strand):  ref base T (33/34);  new base C (2/3);  total (35/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 = 6.58e-01

GAATCGGTGAAGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAG  >  NC_000913/3866321‑3866615                                                                                                                                                     |                                                                                                                                                   gAATCGGTGAAGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGt                                                                                                                                                    <  2:445552/149‑1 (MQ=255)     cGGTGAAGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGc                                                                                                                                                <  3:379431/149‑1 (MQ=255)      ggTGAAGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGGTAGTTAGCt                                                                                                                                               <  1:512592/149‑1 (MQ=255)       gTGAAGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTa                                                                                                                                              <  1:310262/149‑1 (MQ=255)           aGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAACTAGCTAttct                                                                                                                                          <  4:463031/149‑3 (MQ=255)            gCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTaa                                                                                                                                         >  1:573255/1‑149 (MQ=255)            gCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTaa                                                                                                                                         >  2:325919/1‑149 (MQ=255)             ctgATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAAc                                                                                                                                        <  1:456233/147‑1 (MQ=255)              aGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAACTAGCTAttctctt                                                                                                                                       <  2:164458/149‑6 (MQ=255)               gATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGACGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATCGCTAGTTAGCTATTTAAcgc                                                                                                                                      >  4:563629/1‑149 (MQ=255)                 tttGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCgg                                                                                                                                    <  2:504343/149‑1 (MQ=255)                   tGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGa                                                                                                                                  >  4:151682/1‑149 (MQ=255)                      cTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGt                                                                                                                               >  2:174693/1‑149 (MQ=255)                          taGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGc                                                                                                                           <  4:261433/149‑1 (MQ=255)                            gcgcACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTg                                                                                                                         >  4:537869/1‑149 (MQ=255)                            gcacacaTTATCACGGGAAAGAACGGGCCGAAATATCTTAAACAGTCGAAAGTTATCGCATTGTTGAGTAATTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTg                                                                                                                         >  4:491938/4‑149 (MQ=255)                                acaTTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAg                                                                                                                     <  2:422607/149‑1 (MQ=255)                                      tCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATAc                                                                                                               >  2:91018/1‑149 (MQ=255)                                                 aaCGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAACTAGCTAttctct                                                                                                                                        <  3:141153/113‑5 (MQ=255)                                                 aaCGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAACTAGCTAttctct                                                                                                                                        >  4:141153/1‑109 (MQ=255)                                                   cGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGAt                                                                                                  <  1:70727/149‑1 (MQ=255)                                                       ccGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGt                                                                                              >  2:11876/1‑149 (MQ=255)                                                          aaaTATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCa                                                                                           <  2:9218/149‑1 (MQ=255)                                                          aaaTATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCa                                                                                           <  2:573255/149‑1 (MQ=255)                                                          aaaTATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCa                                                                                           <  4:159886/149‑1 (MQ=255)                                                          aaaTATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCa                                                                                           >  2:337389/1‑149 (MQ=255)                                                               tCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTc                                                                                      >  1:517904/1‑149 (MQ=255)                                                                 ttAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCAt                                                                                    >  3:377447/1‑149 (MQ=255)                                                                 ttAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCAt                                                                                    >  4:2658/1‑149 (MQ=255)                                                                   aaaCAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTa                                                                                  >  2:105502/1‑149 (MQ=255)                                                                     aCAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACg                                                                                >  4:486644/1‑149 (MQ=255)                                                                       aGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGaa                                                                              <  1:109086/149‑1 (MQ=255)                                                                           gAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTaa                                                                          <  1:91018/149‑1 (MQ=255)                                                                             aaTTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAAt                                                                        <  4:215824/149‑1 (MQ=255)                                                                                ttAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCaa                                                                     <  1:174693/149‑1 (MQ=255)                                                                                   aCGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAatat                                                                  <  1:337389/149‑1 (MQ=255)                                                                                        ttgttgAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGt                                                             >  1:275955/1‑149 (MQ=255)                                                                                           ttgAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTaaa                                                          <  3:53302/149‑1 (MQ=255)                                                                                             gAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAAcc                                                        >  4:442619/1‑149 (MQ=255)                                                                                               gtaagtaaGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGt                                                      >  2:415574/1‑149 (MQ=255)                                                                                                 aagtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGttt                                                    >  4:63659/1‑149 (MQ=255)                                                                                                  agtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCa                                                               <  2:354083/137‑1 (MQ=255)                                                                                                  agtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCa                                                               >  1:354083/1‑137 (MQ=255)                                                                                                  agtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTa                                                   <  1:15554/149‑1 (MQ=255)                                                                                                   gtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTAc                                                  <  3:151682/149‑1 (MQ=255)                                                                                                   gtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTAc                                                  >  1:362432/1‑149 (MQ=255)                                                                                                   gtaagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTAc                                                  <  4:226526/149‑1 (MQ=255)                                                                                                    taagACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACg                                                 >  3:262809/1‑149 (MQ=255)                                                                                                       gACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAg                                              <  3:442619/149‑1 (MQ=255)                                                                                                       gACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAg                                              <  2:362432/149‑1 (MQ=255)                                                                                                        aCAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAgg                                             >  2:186681/1‑149 (MQ=255)                                                                                                         cAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGt                                            >  4:157813/1‑149 (MQ=255)                                                                                                         cAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGt                                            >  2:414935/1‑149 (MQ=255)                                                                                                          aaaaaaaGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGtt                                           <  1:186681/149‑1 (MQ=255)                                                                                                            aaaaaGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTgc                                         >  1:28274/1‑149 (MQ=255)                                                                                                                     caccaGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGAc                                <  3:305580/149‑1 (MQ=255)                                                                                                                            aTGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCata                         <  3:537869/149‑1 (MQ=255)                                                                                                                            aTGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCata                         <  2:275955/149‑1 (MQ=255)                                                                                                                            aTGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCata                         >  1:503323/1‑149 (MQ=255)                                                                                                                             tGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCatat                        <  1:536492/149‑1 (MQ=255)                                                                                                                                 ggggCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTc                    <  3:486644/149‑1 (MQ=255)                                                                                                                                  gggCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTca                   >  4:148895/1‑149 (MQ=255)                                                                                                                                  gggCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTca                   <  1:105502/149‑1 (MQ=255)                                                                                                                                      aaaGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTcagcca               >  3:62428/1‑149 (MQ=255)                                                                                                                                      aaaGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTcagcca               >  3:267414/1‑149 (MQ=255)                                                                                                                                         ttaaaTAGCTAACTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCg            >  2:271971/4‑149 (MQ=255)                                                                                                                                          agaATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGt           <  3:2658/149‑1 (MQ=255)                                                                                                                                              tAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAg       <  1:11876/149‑1 (MQ=255)                                                                                                                                               aGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGc      >  2:123728/1‑149 (MQ=255)                                                                                                                                                gCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCt     >  4:110182/1‑149 (MQ=255)                                                                                                                                                  tAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCa   <  3:252450/149‑1 (MQ=255)                                                                                                                                                  tAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCa   >  1:76368/1‑149 (MQ=255)                                                                                                                                                   aGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAg  <  4:76368/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GAATCGGTGAAGCAGATTTGCCTTATAGCGCACATTATCACGGGAAAGAACGTGCCGAAATATCTTAAACAGTCGAAATTTAACGCATTGTTGAGTAAGTAAGACAAAAAAAGCCCCACCAGAATGGCGGGGCAAAGAGAATAGCTAGTTAGCTATTTAACGCGGGACGTTCGCTGATAGCGATACGCTGCGCTGCGATGGGTTCAGGCTCATTACGAATTAAATCAATATGCAGTAAACCGTTTACGAAGGTTGCGCCAGAGACTTCCATATTTTCAGCCAGCGTAAAGCTCAG  >  NC_000913/3866321‑3866615

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