Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	758,466	G→T	6.1%	intergenic (+61/‑240)	sdhB → / → sucA	succinate dehydrogenase, FeS subunit/2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	758,466	0	G	T	6.1%	229.3 / 5.0	66	intergenic (+61/‑240)	sdhB/sucA	succinate dehydrogenase, FeS subunit/2‑oxoglutarate decarboxylase, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base G (27/35);  new base T (2/2);  total (29/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 = 5.31e-01

AACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACTACGTAGA  >  NC_000913/758319‑758612                                                                                                                                                    |                                                                                                                                                   aaCTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTg                                                                                                                                                    <  2:427332/148‑1 (MQ=255)   cTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGaa                                                                                                                                                  <  1:126064/148‑1 (MQ=255)   cTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAAc                                                                                                                                                 <  4:246953/149‑1 (MQ=255)    tGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAAc                                                                                                                                                 <  1:901287/148‑1 (MQ=255)          aGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTa                                                                                                                                           >  2:84998/1‑148 (MQ=255)           gtgtATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                         >  4:394793/1‑149 (MQ=255)           gtgtATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATc                                                                                                                                         >  3:782073/1‑149 (MQ=255)             gtATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                       >  4:509374/1‑149 (MQ=255)                tGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                    <  4:25234/149‑1 (MQ=255)                 gTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                    <  1:87204/148‑1 (MQ=255)                 gTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCt                                                                                                                                   >  3:93280/1‑149 (MQ=255)                  tCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTa                                                                                                                                  >  4:44254/1‑149 (MQ=255)                       aGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAg                                                                                                                              <  1:614387/148‑1 (MQ=255)                             tGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTAc                                                                                                                       <  3:694012/149‑1 (MQ=255)                               aaCCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTAccc                                                                                                                     <  4:782073/149‑1 (MQ=255)                                   cGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTa                                                                                                                  >  2:309896/1‑148 (MQ=255)                                         gcgcCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGa                                                                                                           <  3:820566/149‑1 (MQ=255)                                                 ggCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAagcgca                                                                                                   <  3:509374/149‑1 (MQ=255)                                                      tatCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATAAGGCGTTCACGCCGCATCCGGcactg                                                                                                                                 >  3:819475/1‑110 (MQ=25)                                                      tatCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATAAGGCGTTCACGCCGCATCCGGcactg                                                                                                                                 <  4:819475/114‑5 (MQ=25)                                                            gTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACTAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAgtc                                                                                        <  4:650199/149‑3 (MQ=255)                                                            gTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                        <  3:855295/149‑1 (MQ=255)                                                             tCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                                                        >  2:362188/1‑148 (MQ=255)                                                                 tgttgtTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCAtt                                                                                                                         <  2:765914/111‑1 (MQ=255)                                                                 tgttgtTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCAtt                                                                                                                         >  1:765914/1‑111 (MQ=255)                                                                  gttgttGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAc                                                                                   >  2:735664/1‑148 (MQ=255)                                                                  gttgttGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAAc                                                                                   >  1:408040/1‑148 (MQ=255)                                                                      ttgCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGa                                                                              <  4:93280/149‑1 (MQ=255)                                                                       tgCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGaa                                                                             <  3:44254/149‑1 (MQ=255)                                                                            cGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGTGTAACAAAGAAATGc                                                                         <  1:309896/148‑1 (MQ=255)                                                                               aaTGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAgg                                                                      >  1:768788/1‑148 (MQ=255)                                                                                   cGTAAACCGTAGGCCTGATATGACGCGCAAGCGTCGCATCAGGCAACCTGTGCCGGATGCGGCGTGATCGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAAt                                                                  <  1:702745/148‑1 (MQ=255)                                                                                    gTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATc                                                                 >  1:833560/1‑148 (MQ=255)                                                                                     tAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCCGGCGTAACAAAGAAATGCAGGAAATCtt                                                               >  4:697959/1‑149 (MQ=255)                                                                                     tAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCtt                                                               >  4:507124/1‑149 (MQ=255)                                                                                     tAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCtt                                                               >  3:237995/1‑149 (MQ=255)                                                                                      aaaCCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCtt                                                               <  1:84998/148‑1 (MQ=255)                                                                                          cGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTaaa                                                           <  1:446840/148‑1 (MQ=255)                                                                                            tAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAAc                                                        <  3:804416/149‑1 (MQ=255)                                                                                                   gATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGcccc                                                  <  1:735664/148‑1 (MQ=255)                                                                                                    aTAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCt                                                 <  2:768788/148‑1 (MQ=255)                                                                                                     tAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTg                                                <  2:408040/148‑1 (MQ=255)                                                                                                     tAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTg                                                <  1:32240/148‑1 (MQ=255)                                                                                                     tAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGa                                               >  4:319683/1‑149 (MQ=255)                                                                                                        gACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGacac                                            <  4:237995/149‑1 (MQ=255)                                                                                                        gACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGaca                                             <  1:901682/148‑1 (MQ=255)                                                                                                                 gCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAg                                    >  2:67453/1‑148 (MQ=255)                                                                                                                  cGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGtt                                  <  3:195027/149‑1 (MQ=255)                                                                                                                    tCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACa                                                                                  <  3:447274/99‑1 (MQ=255)                                                                                                                    tCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACa                                                                                  >  4:447274/1‑99 (MQ=255)                                                                                                                         tCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAg                           <  3:394793/149‑1 (MQ=255)                                                                                                                         tCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAg                           >  4:430733/1‑149 (MQ=255)                                                                                                                           aGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGtt                                  >  3:71591/1‑140 (MQ=255)                                                                                                                           aGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGtt                                  <  4:71591/140‑1 (MQ=255)                                                                                                                           aGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGt                         >  3:337807/1‑149 (MQ=255)                                                                                                                             gCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGtt                        >  1:667858/1‑148 (MQ=255)                                                                                                                             gCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGtt                        <  1:292739/148‑1 (MQ=255)                                                                                                                             gCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTc                       >  4:703111/1‑149 (MQ=255)                                                                                                                                 ccAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCtt                    <  2:81945/148‑1 (MQ=255)                                                                                                                                       ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGa                                                                     <  1:149782/93‑1 (MQ=25)                                                                                                                                       ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGa                                                                     >  2:149782/1‑93 (MQ=25)                                                                                                                                            tGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACt         <  2:833560/148‑1 (MQ=255)                                                                                                                                             gcggcgTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCt                                                 >  4:103113/1‑107 (MQ=255)                                                                                                                                             gcggcgTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCt                                                 <  3:103113/107‑1 (MQ=255)                                                                                                                                              cggcgTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGcccc                                                  >  2:837483/1‑105 (MQ=255)                                                                                                                                              cggcgTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGcccc                                                  <  1:837483/105‑1 (MQ=255)                                                                                                                                                   tGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACTACGTAGa  >  2:716743/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                   AACTGCGTCAGTGTATGTCCGAAGGGGCTGAACCCGACGCGCGCCATCGGCCATATCAAGTCGATGTTGTTGCAACGTAATGCGTAAACCGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCAACCAGTGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGTCATTACCCGTAGGCCTGATAAGCGCAGCGCATCAGGCGTAACAAAGAAATGCAGGAAATCTTTAAAAACTGCCCCTGACACTAAGACAGTTTTTAAAGGTTCCTTCGCGAGCCACTACGTAGA  >  NC_000913/758319‑758612

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