Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,510,395	0	C	A	42.9%	77.4 / 30.1	64	S131R (AGC→AGA)	ydcR	putative DNA‑binding transcriptional regulator and putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (24/12);  new base A (0/27);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 8.32e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTAAAC  >  NC_000913/1510248‑1510545                                                                                                                                                         |                                                                                                                                                        tCGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCa                                                                                                                                                        >  1:326896/1‑149 (MQ=255)      aTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGc                                                                                                                                                >  1:616744/1‑151 (MQ=255)       tGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGcc                                                                                                                                               >  1:998503/1‑151 (MQ=255)        gAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGcc                                                                                                                                               <  2:585162/150‑1 (MQ=255)         aaGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGACCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGcc                                                                                                                                               <  2:19085/149‑1 (MQ=255)                      cAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTg                                                                                                                                >  2:92710/1‑151 (MQ=255)                        actcTTATATTTTTGATTTGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCTTTT‑CCCGC‑CCCGCGACTTTTCCCCCTCC‑AACAACAAAC‑CCGCTCGCGGGCGCAGGT‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGAt                                                                                                                              <  2:716802/148‑1 (MQ=255)                           cTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTg                                                                                                                            >  2:48415/1‑150 (MQ=255)                            ttATTTTTTTGATTTGTTGCGGGCCAGTCGCGTCCCGTCGGTCTTTCCTTTTGCCTCGCCCTTC‑CCCGC‑CCCGCGCTTTTCCCCCCCCC‑AACAACAAAC‑CCGCTCGCTGGCGCAGGT‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGAt                                                                                                                              <  2:48497/147‑1 (MQ=255)                            ttATTTTTTTGATTTGTTGCGGGCCAGTCGCGACCCGGCGGTCGTCCCTTTTCCCTCGCCCTTC‑CCCGC‑CCCGCGTTTTTCCCCCCTCC‑AACAACAAAC‑CCGCTCGTTGGCGCAGGT‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTg                                                                                                                            <  2:450678/149‑1 (MQ=255)                             tataTTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGaa                                                                                                                          >  2:268429/1‑150 (MQ=255)                              atatTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGaaa                                                                                                                         >  2:299743/1‑150 (MQ=255)                               tttttttGATTTTTTGCGGGCCAGCCGCGATCCGCGGGTCGTCCCTTTCGCCTGGGCCTTC‑CCCGC‑CCCGCGACTTTCCCCCCCCC‑AAAAACAAAC‑CCTCCCGCGGGCGCGGGA‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGAtt                                                                                                                             <  1:398537/143‑1 (MQ=255)                               tatTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAAc                                                                                                                       >  2:181734/1‑151 (MQ=255)                                 ttttttattggtgccAGCCCAGCGGCGTTCCGTCGTTCTTTCCTTTTCCCGGTGCCTTC‑CCCCC‑CCCGCGATTTCCCCCCCCCC‑AACAACAAAC‑CCGCTCGCGGGCGAGGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGaa                                                                                                                          <  1:637257/132‑1 (MQ=255)                                 tttttGATATTTGGCAGCCCGGGCGCGACCCGTCGGTCGTTCCGTTTGCCCGGGCTTCC‑CCCGA‑CCCGCGACTTTCCCCCCCCC‑AACAACAAAA‑CCGCTCGCTGGCGCAGGT‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAAc                                                                                                                       <  1:745159/149‑1 (MQ=255)                                 tttttGATATGTGGCGGGCCGGCCGCGTCCCGTCGGTCGTCCCGTTTCCCTCGGCCTTC‑CCCGC‑CCCGCGTCTTTCCCCCCCCC‑AACAACAAAC‑CCGCGCGTGGGCGCAGAA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACt                                                                                                                      <  1:804514/150‑1 (MQ=255)                                  ttttAATTTGTTGCGGGCCAGTCGCGACCCGTCGTTCTTTCCGTTTCCCTCGGCCTTC‑CCCGC‑CCCGCGACTTTCCCCCCCCC‑AACAACAAAC‑CCCCTCGCGGGCGCAGGT‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACtt                                                                                                                     <  2:284532/150‑1 (MQ=255)                                      gATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGc                                                                                                                >  2:680867/1‑151 (MQ=255)                                        ttttttGCGGGCCAGCCGCGTCCCGTCGGTCTTCCCTTTTCCCTCGCCCTTC‑CCCGC‑CCCGCGATTTTCCCCCCCCA‑AACAACTAAC‑CCGCTCGCTGGCGCAGGA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTAcc                                                                                                                  <  2:447346/143‑1 (MQ=255)                                        tatGTTGCGGCCAAGCCGCGACCCGTGGGTCTTTCCGTTTCCCCCGCCCTTC‑CCCGC‑CCCGCGCTTTTCCCCCCCCC‑AACAACTAAC‑CCGCTCGCGGGCGCAGGA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTAcc                                                                                                                  <  2:380689/147‑1 (MQ=255)                                          tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAgg                                                                                                            >  1:96404/1‑151 (MQ=255)                                            tggcgggccggccGCGACCCGTCGGTCTTCCCGTTCGCCTCGCCCTTC‑CCCGC‑CCCGCGTTTTTCCCCCCCCC‑AACAACAAAC‑CCGCCGGCGGGCGCGGGT‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCg                                                                                                                 <  2:721311/139‑1 (MQ=255)                                             tGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGaa                                                                                                          >  1:321194/1‑150 (MQ=255)                                                  ccccgccgcgccccGTGGGTCGTCCCTTTTCCCCGGGCTTTC‑CCCGC‑CCCGCGCTTTCCCCCCCCCC‑AACAATAAAC‑CCGCCCGTGGGCGGGGTA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGa                                                                                                           <  1:348752/132‑1 (MQ=255)                                                    cggccgccacccGTCGTTCTTCCCGTTTCCCCCGCCTTTC‑CCCGC‑CCCGCG‑CTTTTCCCCCCCCAAAAAACAAAC‑CCGCCCGCGGCCGCAGGT‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAg                                                                                                    <  2:470680/139‑1 (MQ=255)                                                       tCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTg                                                                                               >  2:849305/1‑151 (MQ=255)                                                         gcgATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCg                                                                                             >  2:759108/1‑151 (MQ=255)                                                           gATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTc                                                                                           >  1:758842/1‑151 (MQ=255)                                                              ccggggttcttcccGTTTCCCGCGCTCTTC‑CCCGC‑CCCGCG‑CTTTTCCCCCCCC‑AAAAACCAAACCCGCGCGGGGGGGCA‑GT‑AAGGAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAgg                                                                                        <  1:578269/139‑1 (MQ=255)                                                              ccGTCGGTC‑TTCCCTTTTCCCCCGCCTTTCCCCGC‑CCCGCGTTTTTCCCCCCCCA‑AACAACAAAC‑CCCCTCGCGGGGGCAGGA‑AAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAgg                                                                                        <  2:829341/151‑1 (MQ=255)                                                               ccgtggttgTTCCCTTTTCCCCCGCCTTTCCCCGC‑CCCGCGATTTTCCCCCCCCC‑AACAAAAAAC‑CCGCCCGCTGGCGCAGGA‑AAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAg                                                                                         <  2:27710/146‑1 (MQ=255)                                                                  cggtcgTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAt                                                                                    >  1:967081/1‑151 (MQ=255)                                                                   gttcttcccttttGCCCCGGCCTTC‑CCCGC‑CCCGCGATTTTCCCCCCCCA‑AAAAAAAAAA‑CCGCCCGCGGGCGCAGGA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAt                                                                                    <  1:586487/140‑1 (MQ=255)                                                                     ttttcccGTTTGCC‑CCGCCTTT‑CCCCC‑CCCCGGGTTTTTCCCCCCCC‑CAAAACAAAACCCCCTCGCGGGCGAAGGA‑AAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAtt                                                                                   <  2:133482/144‑1 (MQ=255)                                                                       gTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCtc                                                                               >  2:1021951/1‑151 (MQ=255)                                                                          ccttttCCCCCGCCCTTC‑CCCGC‑CCCGCG‑CTTTTCCCCCCCCAAAAAACAAAC‑CCGCTCGCTGGCGCAGAT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCtc                                                                               <  2:147541/145‑1 (MQ=255)                                                                           cGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCtcgtc                                                                            >  2:427859/1‑150 (MQ=255)                                                                           cGTTTGCCCCGGCCTTTCCCCGC‑CCCCCGTCTTTCCCCCCCCA‑AAAAAAAAAC‑CCGCTCGCGGGCGCAGTT‑AAAAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAg                                                                                         <  1:169345/138‑1 (MQ=255)                                                                                   tCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCt                                                                   >  2:841834/1‑151 (MQ=255)                                                                                      gCCTTTCCCCGC‑CCCGCGACTTTCCCCCCCCA‑AACAACAAAC‑CCGCCCGCGGGGGCAGGA‑AAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCt                                                                   <  1:649258/149‑1 (MQ=255)                                                                                      gCCTTT‑CCCCCCCCCGCGTTTTTCCCCCCCCC‑AACAACAAAC‑CCGCCCGCGGGCGCAGGA‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCt                                                                   <  1:497005/149‑1 (MQ=255)                                                                                          ttccccgc‑cccgcG‑CTTTTCCCCCCCCAAAAAACAAAC‑CCGCCCGCGGGCGCAGGA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCtt                                                                  <  2:1010550/143‑1 (MQ=255)                                                                                             cccgc‑cccgcGCTTTTCCCCCCCCC‑AACAACAACC‑CCGCGCGCGGGCGCGGAT‑AAAAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACa                                                               <  1:964857/146‑1 (MQ=255)                                                                                             cccccccccactccTTTCCCCCCCCC‑AAAAACGACC‑CCGCCCCTTGGCACAGAA‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAggg                                                            <  2:638349/137‑1 (MQ=255)                                                                                                   cccGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCa                                                       >  1:142616/1‑149 (MQ=255)                                                                                                   cccGCGACTTTCCCCCCCCA‑AACAACAAAC‑CCGCCCGCGGGCGCAGGA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCAcc                                                     <  2:478205/151‑1 (MQ=255)                                                                                                           ttttCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTcc                                             >  2:675172/1‑151 (MQ=255)                                                                                                             ttCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtg                                           >  1:800674/1‑151 (MQ=255)                                                                                                               ccccccCA‑AACAACAAAC‑CCGCCCGGTGGCGCAGGA‑AAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtg                                           <  2:273146/149‑1 (MQ=255)                                                                                                                        aacaacAAAC‑CCGCGCGCGGGCGCAGGA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGt                                 <  2:772877/151‑1 (MQ=255)                                                                                                                          caacTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTc                                >  2:50962/1‑150 (MQ=255)                                                                                                                           aacTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCAt                              >  1:893198/1‑151 (MQ=255)                                                                                                                           aacAAAC‑CCGCTCGCTGGCGCAGGA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAAtt                                   <  2:317888/146‑1 (MQ=255)                                                                                                                               aaa‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CTGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCAtt                             <  1:292826/148‑1 (MQ=255)                                                                                                                                 accccccccctGGCGCAGGA‑AACCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTGCCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCCGATGAAATTGTCATTACt                          <  1:306111/141‑1 (MQ=255)                                                                                                                                   ccGCTCGCGGGCGCAGAA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGc                        <  2:615539/150‑1 (MQ=255)                                                                                                                                      ctcgcgGGC‑CAGGTAAAACAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACt                          <  1:514053/145‑1 (MQ=255)                                                                                                                                           gggCGGAGGA‑AAAAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGc                  <  2:486090/147‑1 (MQ=255)                                                                                                                                            ggggcAGGA‑AAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTa              <  2:848041/148‑1 (MQ=255)                                                                                                                                            ggCGCAGGA‑AAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTa              <  1:953339/151‑1 (MQ=255)                                                                                                                                                cAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAgg          >  1:591287/1‑151 (MQ=255)                                                                                                                                                      aaGCAAAG‑CCGCCCCACCTATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCAt       <  1:33534/149‑1 (MQ=255)                                                                                                                                                      aaGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCAtt      <  2:305396/150‑1 (MQ=255)                                                                                                                                                         cAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTAAAc  >  2:392593/1‑151 (MQ=255)                                                                                                                                                         cAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTAAAc  >  2:380041/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                        TCGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTCGGCCTTT‑CCCGA‑CCCGCGACTTTTCCCCCTCC‑AACAACTAAA‑CCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATTAAAC  >  NC_000913/1510248‑1510545

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