Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,623,523	G→T	5.1%	Q199K (CAG→AAG)	lplA ←	lipoate‑protein ligase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,623,523	0	G	T	5.1%	275.9 / 4.2	78	Q199K (CAG→AAG)	lplA	lipoate‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base G (35/39);  new base T (2/2);  total (37/41)
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.70e-01

ACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAG  >  NC_000913/4623377‑4623662                                                                                                                                                   |                                                                                                                                            aCCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGc                                                                                                                                            <  2:613716/148‑1 (MQ=255)  ccGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTc                                                                                                                                          >  4:204494/1‑149 (MQ=255)     aaGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCAt                                                                                                                                        <  1:18482/148‑1 (MQ=255)     aaGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCAt                                                                                                                                        >  2:68495/1‑148 (MQ=255)      aattCCACTCCCATCTACTCGGGCGGGCAAAGGTTTCGGCGAAGTTTGGAAAGTCTGGCGTTTTGTTCGGGGAGATCATTTCCACTCCCTCACGCTCGCAATAATGGGCGAAAAAGGACTCGGTTATGGCCTCGCAAACCTGCTCATg                                                                                                                                       <  2:172289/146‑1 (MQ=255)      aGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATg                                                                                                                                       <  1:261029/148‑1 (MQ=255)        ttCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATggg                                                                                                                                     <  2:133420/148‑1 (MQ=255)         tCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGt                                                                                                                                    >  1:376893/1‑148 (MQ=255)                ccAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCg                                                                                                                            >  3:750689/1‑149 (MQ=255)                 cAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATg                                                                                                                                       <  1:689013/137‑1 (MQ=255)                 cAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATg                                                                                                                                       >  2:689013/1‑137 (MQ=255)                     tACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCaa                                                                                                                        <  1:79768/148‑1 (MQ=255)                        tctGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGc                                                                                                                    >  3:20520/1‑149 (MQ=255)                          tGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCt                                                                                                                   >  2:226870/1‑148 (MQ=255)                           ggCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTc                                                                                                                  <  2:13912/148‑1 (MQ=255)                           ggCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTc                                                                                                                  <  1:224900/148‑1 (MQ=255)                               ggCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGa                                                                                                             <  4:643346/149‑1 (MQ=255)                               ggCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGa                                                                                                             <  4:876759/149‑1 (MQ=255)                                  aaaGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAgg                                                                                                           <  1:226870/148‑1 (MQ=255)                                      gTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGAt                                                                                                                                 <  4:687442/122‑1 (MQ=255)                                      gTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGAt                                                                                                                                 >  3:687442/1‑122 (MQ=255)                                              gAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCggg                                                                                              >  4:327195/1‑149 (MQ=255)                                               aaGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATg                                                                                                                                       <  1:882401/107‑1 (MQ=255)                                               aaGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATg                                                                                                                                       >  2:882401/1‑107 (MQ=255)                                                 gTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAAc                                                                                           >  4:659228/1‑149 (MQ=255)                                                 gTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAAc                                                                                           <  3:810287/149‑1 (MQ=255)                                                   ttGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACg                                                                                          >  2:550004/1‑148 (MQ=255)                                                    tGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTa                                                                                        <  3:782505/149‑1 (MQ=255)                                                           tCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACg                                                                                  <  1:541702/148‑1 (MQ=255)                                                               gCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATg                                                                             >  3:383838/1‑149 (MQ=255)                                                                cGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATg                                                                             >  2:319485/1‑148 (MQ=255)                                                                cGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATg                                                                             >  1:65937/1‑148 (MQ=255)                                                                cGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGc                                                                            <  4:20520/149‑1 (MQ=255)                                                                cGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGc                                                                            >  3:177068/1‑149 (MQ=255)                                                                cGTGTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATg                                                                             >  2:807246/1‑148 (MQ=255)                                                                  ttttGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATgg                                                                                                                                      >  4:476354/1‑89 (MQ=255)                                                                  ttttGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATgg                                                                                                                                      <  3:476354/89‑1 (MQ=255)                                                                      gTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTcg                                                                      >  3:263789/1‑149 (MQ=255)                                                                          ggggAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGACTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTcgccg                                                                   <  2:497112/148‑1 (MQ=255)                                                                           gggAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTcgccgc                                                                  >  1:594794/1‑148 (MQ=255)                                                                            ggAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCgccgcc                                                                 <  2:251332/148‑1 (MQ=255)                                                                            ggAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCgccgcc                                                                 >  2:491864/1‑148 (MQ=255)                                                                            ggAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCa                                                                <  3:489464/149‑1 (MQ=255)                                                                             gagaTGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAAc                                                                                           >  4:2638/1‑121 (MQ=255)                                                                             gagaTGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAAc                                                                                           <  3:2638/121‑1 (MQ=255)                                                                                    tttCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCgg                                                                                                                >  1:544833/1‑93 (MQ=255)                                                                                    tttCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCgg                                                                                                                <  2:544833/93‑1 (MQ=255)                                                                                    tttCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCttt                                                        <  3:327195/149‑1 (MQ=255)                                                                                     ttCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCttt                                                        >  1:787461/1‑148 (MQ=255)                                                                                           ttCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTTCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCg                                                 >  3:109674/1‑149 (MQ=255)                                                                                            tCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTTCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCg                                                 >  1:867834/1‑148 (MQ=255)                                                                                               aCGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGAtt                                             >  3:905732/1‑149 (MQ=255)                                                                                                 gcgcTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTg                                            <  2:460434/148‑1 (MQ=255)                                                                                                     tCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGAt                                        <  2:376893/148‑1 (MQ=255)                                                                                                       gCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAg                                      <  1:491864/148‑1 (MQ=255)                                                                                                         cATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTgag                                          <  2:437431/142‑1 (MQ=255)                                                                                                         cATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTgag                                          >  1:437431/1‑142 (MQ=255)                                                                                                          aTAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTg                                  >  3:140704/1‑149 (MQ=255)                                                                                                           tAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGAtt                                             >  3:301432/1‑137 (MQ=255)                                                                                                           tAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGAtt                                             <  4:301432/137‑1 (MQ=255)                                                                                                            aaTGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGc                                 >  1:523165/1‑148 (MQ=255)                                                                                                             aTGTGCGAAAAAGGCCTCGCTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCCCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCa                               >  4:822488/1‑149 (MQ=255)                                                                                                             aTGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGcc                                <  2:787461/148‑1 (MQ=255)                                                                                                               gggCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAgg                             >  3:202104/1‑149 (MQ=255)                                                                                                                ggCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAgg                             <  1:68495/148‑1 (MQ=255)                                                                                                                     aaaaGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCt                        <  2:520627/148‑1 (MQ=255)                                                                                                                       aaGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGATTGCCAGGCGGCTGAg                     <  4:383838/149‑1 (MQ=255)                                                                                                                           ccTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCg                 <  4:750689/149‑1 (MQ=255)                                                                                                                            cTCGGTTATGGCCTCGCAAACCTTCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCgg                <  4:109674/149‑1 (MQ=255)                                                                                                                            cTCGGTTATGGCCTCGCAAACCTTCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCg                 <  2:867834/148‑1 (MQ=255)                                                                                                                               ggTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCa              <  1:550004/148‑1 (MQ=255)                                                                                                                                    tGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAg         <  1:319485/148‑1 (MQ=255)                                                                                                                                     ggCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTa       >  3:680732/1‑149 (MQ=255)                                                                                                                                     ggCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGCTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTa       >  3:219372/1‑149 (MQ=255)                                                                                                                                          cGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGGCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAg  <  4:177068/149‑1 (MQ=255)                                                                                                                                          cGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCaa   >  2:759651/1‑148 (MQ=255)                                                                                                                                                 cTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCttt                                                        <  4:195004/88‑1 (MQ=255)                                                                                                                                                 cTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCttt                                                        >  3:195004/1‑88 (MQ=255)                                                                                                                                                   |                                                                                                                                            ACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTTTGGCAAGTCTGGCGTTTTGTTCGGGGAGATGATTTCCGCTTCCACGCGCTCGCCATAATGGGCGAAAAAGGCCTCGGTTATGGCCTCGCAAACCTGCTCATGGGTGATCCCCGGCAACAGCTCGGTGAGGTTGGTCACGCGGGAACGTACCGACGTAATGCCTTTCGCCGCCAGTTTCTTTTTATCCGGATTGAGATAGTTTGCCAGGCGGCTGAGGTCGGCATTGAGTAGCAAG  >  NC_000913/4623377‑4623662

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