Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	298,960	0	A	C	40.9%	115.5 / 22.9	66	Y659D (TAT→GAT)	paoC	PaoABC aldehyde oxidoreductase, Moco‑containing subunit
Reads supporting (aligned to +/- strand):  ref base A (36/3);  new base C (0/27);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.66e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC  >  NC_000913/298811‑299107                                                                                                                                                         |                                                                                                                                                    gTTATACACCGCGTTGGCGATAGCCGCGCTCGCGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATAc                                                                                                                                                    <  2:120884/151‑1 (MQ=255) gTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATAc                                                                                                                                                    >  1:205201/1‑151 (MQ=255)  ttATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATAcc                                                                                                                                                   >  1:257471/1‑151 (MQ=255)   tataCACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATAccc                                                                                                                                                  >  1:619052/1‑151 (MQ=255)       cacCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGcc                                                                                                                                              <  1:312503/151‑1 (MQ=255)             gTTGGCGATGGCCGCGCTCCCGCCGCACGGGCCCAGCTCACCGACCCCTTGGGCCTTCATCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCTCCTGTTTGGGG‑ATATCCGAATGACCCGG‑CC‑CCTAATCCCCCGCCata                                                                                                                                           <  2:326550/148‑1 (MQ=255)               tGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTg                                                                                                                                      >  2:101134/1‑151 (MQ=255)                    aTAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTa                                                                                                                                  >  2:548379/1‑150 (MQ=255)                       gCCGCGCTCACGCCGCACAGCCCCAGCCCCCCGACACTTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTATCATCCAGAAAAATCCCCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CC‑CCTCATCCCCCGCCATATCGTGATTAACg                                                                                                                               <  1:534401/150‑1 (MQ=255)                        ccGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGa                                                                                                                              >  1:705088/1‑150 (MQ=255)                         cgcgCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGa                                                                                                                              >  1:595491/1‑149 (MQ=255)                            gcTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGt                                                                                                                           >  2:295137/1‑149 (MQ=255)                             cTCCCGCCGCACGGCCCCAGCTCCCCGCCCCCTTTGGCCTTCATCGGGGAGGATAGGGGGTCGGTACCACCCAGAAAAACCCCCTCCTTTTTGGGG‑ATA‑CCCAAGGGACCGGCCC‑CCTCATCCCCCGCCATATCGTGATTAACGAAg                                                                                                                            <  2:177758/147‑1 (MQ=255)                             cTCACGCCGCACAGGCCCAGCTCACCGACCCCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCTCCTGTTTGGGA‑ATACCCGCATGAACCGG‑CA‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAg                                                                                                                         <  2:865609/150‑1 (MQ=255)                                      cacaGGCCCAGCTCCCCGACCCTTTGGGCCTTCATGGGGGAGGATATGGGGTGGGTATCATCCAGAAAAATCCCCTCCTGTTTTGGA‑ATATCCGAATGAACCGG‑CC‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGt                                                                                                               <  1:213428/151‑1 (MQ=255)                                        caGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTc                                                                                                              >  1:289603/1‑150 (MQ=255)                                             ccAGCTCCCCGCCCCCTTTGGCCTTCACCGGGGAGGATATGGGGTCGGTATCACCCAGAAAAACCCCCTCCTTTTTGGGG‑ATACCCGAATGACCCGC‑CC‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCAc                                                                                                        <  2:218720/151‑1 (MQ=255)                                                   ccccGACCCCTTGGGCCTTCAGCGGGGAGGATATGGGGTCGGTACCATCCAGAAAAATCCCCCCCTGTTTTGGG‑ATATCCGCATGACCCGG‑CA‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGcc                                                                                                    <  1:626544/147‑1 (MQ=255)                                                         cacCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcctc                                                                                            >  2:406374/1‑151 (MQ=255)                                                             ttgggCCTTCACCGGGGAGGATAGGGGGTGGGTATCACCCAGAAAAACCCCCTCCTTTTTTGGG‑ATATCCGCATGACCCGG‑CC‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGctcctc                                                                                            <  1:298578/144‑1 (MQ=255)                                                                ggCCTCCAGGGGGGAGAATAGGGGGTCGGTATCACCCAGAAAAACCCCCCCCTGTTTTGGA‑ATTCCCGAATGAACCGC‑CC‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAg                                                                                      <  1:581992/150‑1 (MQ=255)                                                                  ccTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGc                                                                                   >  2:317533/1‑151 (MQ=255)                                                                      cATCGGGGAGGATAGGGGGTGGTTATCCCCCAGAAAACCCCCCTCCTGTTTTGGG‑ATACCCGCATGACCCGC‑CC‑CCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcgc                                                                                <  1:570771/150‑1 (MQ=255)                                                                          ggggAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgccc                                                                            >  2:581480/1‑150 (MQ=255)                                                                          ggggAGGATAGGGGGTCGTTACCACCAAGAAAAACCCCCCCCTGTTTTGGG‑ATACCCGAATAACCCGC‑CC‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCcg                                                                                 <  1:132152/145‑1 (MQ=255)                                                                           gggagattttGGGGTCGGTATCACCCAGGAAAATCCCCCCTTTTTTGGGG‑ATACCCGAATGAACCGG‑CC‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                          <  1:318981/142‑1 (MQ=255)                                                                           gggAGGATATGGGGGCGTTACCACCAAGGAAAACCCCCCCCGGTTTGGGA‑ATACCCGCGTGACCCGA‑CT‑CCTCACCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccat                                                                          <  2:670174/151‑1 (MQ=255)                                                                              aggaTATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcc                                                                       >  2:866835/1‑151 (MQ=255)                                                                                 atatGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgcccat                                                                    >  1:270203/1‑151 (MQ=255)                                                                                     ggggTCGTTACCACCCAGAAAACCCCCCTCTTGTTTTGGG‑ATACCCGCAGGACCCGG‑CC‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCgcccatgccca                                                                     <  2:873986/146‑1 (MQ=255)                                                                                     ggggTCGGTACCACCCAGAAAAACCCCCTCTTGTTTGGGA‑ATACCCGCAGGAACCGC‑CC‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATa                                                                   <  2:690880/148‑1 (MQ=255)                                                                                      gggTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTc                                                                >  1:758750/1‑150 (MQ=255)                                                                                           gttATCAACAAGAAAAACCCCCCCCTTTTTGGGG‑ATACCCGCAGGACCCGC‑CC‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCAtt                                                             <  1:577858/146‑1 (MQ=255)                                                                                           ggTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcg                                                          >  1:85753/1‑151 (MQ=255)                                                                                           ggTATCACCCAGAAAATTCCCCCCTTGTTTTGGA‑ATACCCGAATGACCCGC‑CC‑CCTAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTg                                                            <  1:927240/149‑1 (MQ=255)                                                                                              atcatcCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc                                                        >  1:502715/1‑150 (MQ=255)                                                                                                catcCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt                                                     >  2:33342/1‑151 (MQ=255)                                                                                                 acccAGAAAAACCCCCCCCTGTTTGGGG‑ATACCCGCATGACCCGC‑CC‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGcc                                                        <  1:417339/145‑1 (MQ=255)                                                                                                   ccAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGAc                                                  >  2:379651/1‑151 (MQ=255)                                                                                                   ccAGAAAAATCCCCTCTTTTTTGGGA‑ATATCCGAATGACCCGC‑CC‑CCCAATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAAt                                                     <  2:654983/148‑1 (MQ=255)                                                                                                    cgagAAAAGCCCCCCCCGTTTTGGG‑GAACCCCAAAGAACCCG‑CCCCCTCCTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTgcg                                                          <  2:339128/140‑1 (MQ=255)                                                                                                     cagAAAACCCCCCCCCGTTTTGGGAATCCCCGCATGACCCGG‑CC‑CCCCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGa                                                   <  2:747736/147‑1 (MQ=255)                                                                                                     aGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCt                                                >  2:519950/1‑151 (MQ=255)                                                                                                      ggAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTg                                               >  1:810697/1‑151 (MQ=255)                                                                                                       gAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTgg                                              >  2:835822/1‑151 (MQ=255)                                                                                                        aaaaTCCCCCCTTTTTTTGGA‑TTACCCCCATAACCCGC‑CC‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTg                                               <  1:943347/149‑1 (MQ=255)                                                                                                         aaacccccccctTTTTGGGG‑ATCCCCGCAAACACCGG‑CC‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCt                                            <  1:137896/142‑1 (MQ=255)                                                                                                            tCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcg                                         >  2:560948/1‑151 (MQ=255)                                                                                                            ccccctcctTTTTGGGG‑ATACCCGCGGGACCGGC‑CC‑CCCCCTCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTg                                               <  2:37160/142‑1 (MQ=255)                                                                                                               cctcttTTTTTGGG‑ATACCCCCAAGACCCGG‑CC‑CCTCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTgcgcgc                                      <  1:742667/151‑1 (MQ=255)                                                                                                                     gTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt                                >  1:14643/1‑151 (MQ=255)                                                                                                                      tttgGGA‑ATACCCGAATGACCCGC‑CC‑CCTCCCCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGtttt                                <  2:706339/150‑1 (MQ=255)                                                                                                                           gg‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGAtt                          >  1:884529/1‑151 (MQ=255)                                                                                                                                  ccGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGAt                    >  1:556809/1‑151 (MQ=255)                                                                                                                                  ccGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGAt                    >  1:396828/1‑151 (MQ=255)                                                                                                                                      aTGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGt                >  2:993576/1‑151 (MQ=255)                                                                                                                                          aCCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTg            >  2:842496/1‑151 (MQ=255)                                                                                                                                           ccGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc           >  2:763024/1‑151 (MQ=255)                                                                                                                                            cGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgc           >  2:893787/1‑150 (MQ=255)                                                                                                                                              g‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcagc        >  2:529371/1‑151 (MQ=255)                                                                                                                                                 cccttCATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTgcagca       <  2:595479/147‑1 (MQ=255)                                                                                                                                                    cTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacac    >  2:451711/1‑151 (MQ=255)                                                                                                                                                    cTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacac    >  1:402951/1‑151 (MQ=255)                                                                                                                                                     tCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacac    >  1:567727/1‑150 (MQ=255)                                                                                                                                                      cATCCCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATAGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacacac  <  2:86378/151‑1 (MQ=255)                                                                                                                                                       aTACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGcacacac  >  2:242571/1‑150 (MQ=255)                                                                                                                                                         |                                                                                                                                                    GTTATACACCGCGTTGGCGATAGCCGCGCTCACGCCGCACAGGCCCAGCTCACCGACACCTTTGGCCTTCATCGGGGAGGATATGGGGTCGGTATCATCCAGGAAAATCACCTCCTGTTTTGGG‑ATATCCGCATGAACCGG‑CA‑CCTCATACCCCGCCATATCGTGATTAACGAAGTAGCCCAAACGGTCATCCACCGCCAGCTCCTCCATCAGTGCCGCGCCCATGCCCATAGTCATTGCGCCAATGACCTGGCTGCGCGCAGTTTTCGGATTCAGGATGCGTCCTGCAGCACACAC  >  NC_000913/298811‑299107

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