Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,042,457	C→A	G30G (GGC→GGA)	yodB →	cytochrome b561 homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,042,457	0	C	A	95.7%	227.4 / 1.4	69	G30G (GGC→GGA)	yodB	cytochrome b561 homolog
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base A (28/38);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

AGGAATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACC  >  NC_000913/2042311‑2042597                                                                                                                                                   |                                                                                                                                             aGGAATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATg                                                                                                                                            <  8:184219/149‑1 (MQ=255)    aaTTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGtt                                                                                                                                         <  2:329339/149‑1 (MQ=255)    aaTTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGtt                                                                                                                                         >  2:132064/1‑149 (MQ=255)        tttAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTcc                                                                                                                                      <  1:132064/148‑1 (MQ=255)           aaGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTaa                                                                                                                                   <  1:288939/148‑1 (MQ=255)            aGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTaaa                                                                                                                                  <  6:355152/148‑1 (MQ=255)               ttCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgta                                                                                                                              <  1:394480/149‑1 (MQ=255)                  gTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtagta                                                                                                                           <  8:168895/149‑1 (MQ=255)                     ggATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtagtagta                                                                                                                        <  6:224817/149‑1 (MQ=255)                     ggATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATGCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtagtagta                                                                                                                        <  2:383960/149‑1 (MQ=255)                      gATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTAc                                                                                                                       <  2:357902/149‑1 (MQ=255)                       aTAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACt                                                                                                                      >  7:226105/1‑149 (MQ=255)                         aGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTa                                                                                                                    >  1:23343/1‑149 (MQ=255)                          gCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTAt                                                                                                                   <  5:133512/149‑1 (MQ=255)                             tGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtag                                                                                                                             <  1:384718/136‑1 (MQ=255)                             tGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGgtag                                                                                                                             >  2:384718/1‑136 (MQ=255)                              gACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATctgc                                                                                                               >  6:363798/1‑149 (MQ=255)                                    gAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGc                                                                                                         >  5:241615/1‑149 (MQ=255)                                      aaaTCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGc                                                                                                         >  6:300475/1‑147 (MQ=255)                                      aaaTCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCga                                                                                                       <  6:396040/149‑1 (MQ=255)                                      aaaTCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTACTAGTACTTATCTGCTGATGCga                                                                                                       <  8:209176/149‑1 (MQ=255)                                         tCACAGGCGCTACAAAAATGAACCGAGTCTCAAAAACTCAAATTTATTTACAATGGATAGCGCTGCTATTCGTTGCAATACCCTATCCCGCGATGGACCGCCGTGGATGGTTTCCTAAAGGTAGTAGTACTGATCTGCTGATGagagc                                                                                                     >  2:351841/1‑147 (MQ=255)                                                 gCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTa                                                                                            <  3:14594/149‑1 (MQ=255)                                                    aCAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACaa                                                                                         >  7:400546/1‑149 (MQ=255)                                                      aaaaaTGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATg                                                                                       <  5:296911/149‑1 (MQ=255)                                                           tGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGt                                                                                  <  5:17671/149‑1 (MQ=255)                                                           tGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGt                                                                                  <  2:23343/149‑1 (MQ=255)                                                             aaCCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGta                                                                                 >  7:110033/1‑148 (MQ=255)                                                             aaCCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGta                                                                                 >  3:316588/1‑148 (MQ=255)                                                                   ttCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGc                                                                                                         >  7:394355/1‑118 (MQ=255)                                                                   ttCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGc                                                                                                         <  8:394355/118‑1 (MQ=255)                                                                   ttCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGt                                                                          <  3:363135/149‑1 (MQ=255)                                                                        aaaaaCTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCgg                                                                                    >  5:119201/1‑134 (MQ=255)                                                                        aaaaaCTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGttt                                                                        <  5:300475/146‑1 (MQ=255)                                                                        aaaaaCTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGttt                                                                        <  8:226105/146‑1 (MQ=255)                                                                        aaaaaCTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGttt                                                                        <  6:101584/146‑1 (MQ=255)                                                                        aaaaaCTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGGTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCgg                                                                                    <  6:119201/134‑1 (MQ=255)                                                                            aCTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGgtgt                                                                  >  6:236279/1‑148 (MQ=255)                                                                              tCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTaa                                                               >  7:168385/1‑149 (MQ=255)                                                                              tCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCggg                                                                                   <  3:152702/129‑1 (MQ=255)                                                                              tCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCggg                                                                                   >  4:152702/1‑129 (MQ=255)                                                                              tCAAATTTATTTACATTGGATAACGCTGATTTTCGTTGCAATAACATATCCCGCGATGGAACTCAGTGGATGGTTTCCTAAACGTAGTATTACTTATCTTCTGATGCGAGAAACACATTACAATGAGGGTATATTAGTTTGGGTGTTaa                                                               >  8:123883/1‑149 (MQ=255)                                                                                aaaTTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATg                                                             <  6:241615/149‑1 (MQ=255)                                                                                  atttatttACATTGGATAACGCTGCTTTTCGTTGAAATAACCTATGCAACGATGGAACTCCGTTGATGGTTTCCTAAAGGTAGTAGTACTTATCTCCTCATGCGATAAACACAATACAATGCGTGTATATTCTTTTGGGTGTATATGtt                                                           >  1:242803/1‑149 (MQ=255)                                                                                       tttaCATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCAc                                                      <  8:360236/149‑1 (MQ=255)                                                                                        ttaCATTGGATAACGCTGCTTTTCGTTGCAAAAACCAATGCCGCGATGGAACTACGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACg                                                     >  8:33375/1‑149 (MQ=255)                                                                                           cATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCt                                                  >  2:378748/1‑149 (MQ=255)                                                                                            aTTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCt                                                  >  3:261216/1‑148 (MQ=255)                                                                                             ttGGATATCGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGATCTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGa                                                <  1:392879/149‑1 (MQ=255)                                                                                                  tAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATta                                             <  5:129205/147‑1 (MQ=255)                                                                                                  tAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATta                                             >  6:129205/1‑147 (MQ=255)                                                                                                    aagctgctTTTCGTTGCAATAACCTATGCAGCGATGGAACTCCGTGGATTGTTTCCTAAAGGTAGTAGTACTTATCTGGTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATaa                                          <  2:182594/146‑1 (MQ=255)                                                                                                    aCGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATaaa                                         <  8:110033/149‑1 (MQ=255)                                                                                                    aCGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATaa                                          <  7:135969/148‑1 (MQ=255)                                                                                                      gctgctTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATaaaa                                        >  6:317081/1‑148 (MQ=255)                                                                                                       ctgctTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAaca                                      <  5:363798/149‑1 (MQ=255)                                                                                                          ctTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACAcc                                    <  8:400546/148‑1 (MQ=255)                                                                                                           ttttCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGt                                  >  6:13710/1‑149 (MQ=255)                                                                                                             ttCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTta                                >  3:110168/1‑149 (MQ=255)                                                                                                                gTTGCAATAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAg                             >  7:87413/1‑149 (MQ=255)                                                                                                                    cAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGAt                         <  8:168385/149‑1 (MQ=255)                                                                                                                     aaTAACCTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCATGCCTGATTATAAAACACCGTTATAGTGATc                        >  8:300807/1‑149 (MQ=255)                                                                                                                           cTATGCCGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTa                  >  4:370313/1‑149 (MQ=255)                                                                                                                                ccGCGATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATtgt              <  5:236279/148‑1 (MQ=255)                                                                                                                                   cgATGGAACTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTgcca          <  5:317081/149‑1 (MQ=255)                                                                                                                                       ggAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTgccaccgc      >  3:77931/1‑149 (MQ=255)                                                                                                                                           cTCCGTGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTgccaccgccacc  <  7:33375/149‑1 (MQ=255)                                                                                                                                                tGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATaa                                          >  3:141841/1‑104 (MQ=255)                                                                                                                                                tGGATGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATaa                                          <  4:141841/104‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                             AGGAATTTTTAAGGTTCGTGGGATAGCTTGACTGTGAAAATCACAGGAGCTACAAAAATGAACCGATTCTCAAAAACTCAAATTTATTTACATTGGATAACGCTGCTTTTCGTTGCAATAACCTATGCCGCGATGGAACTCCGTGGCTGGTTTCCTAAAGGTAGTAGTACTTATCTGCTGATGCGAGAAACACATTACAATGCGGGTATATTCGTTTGGGTGTTAATGTTTTCACGCCTGATTATAAAACACCGTTATAGTGATCCTTCTATTGTGCCACCGCCACC  >  NC_000913/2042311‑2042597

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