BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,064,160	T→C	15.4%	L42P (CTG→CCG) ‡	yccE →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,064,160	0	T	C	15.4%	374.2 / 51.9	164	L42P (CTG→CCG)	yccE	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base T (59/77); new base C (17/8); total (78/86)
Fisher's exact test for biased strand distribution p-value = 2.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

GTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGG  >  NC_000913/1064061‑1064258

gTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACt                                                                                                    >  2:953063/1‑100 (MQ=255)

gTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACt                                                                                                    >  2:4377130/1‑100 (MQ=255)

gTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACt                                                                                                    >  1:3165345/1‑100 (MQ=255)

  tGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTgg                                                                                                  <  2:1865854/100‑1 (MQ=255)

  tGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTgg                                                                                                  >  2:4167755/1‑100 (MQ=255)

  tGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTgg                                                                                                  <  1:2508284/100‑1 (MQ=255)

      cTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGgtgtg                                                                                              <  2:2691622/100‑1 (MQ=255)

        gCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGtt                                                                                            >  2:2639386/1‑100 (MQ=255)

         cAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGttt                                                                                           <  2:4397448/100‑1 (MQ=255)

          aaaTAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGtttt                                                                                          <  2:2955892/100‑1 (MQ=255)

           aataatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTa                                                                                         <  2:2121936/100‑1 (MQ=255)

            ataatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttca                                                                                        <  2:3633363/100‑3 (MQ=255)

             taatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATg                                                                                       <  2:3649607/100‑1 (MQ=255)

             taatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATg                                                                                       >  2:4052387/1‑100 (MQ=255)

             taatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAAAACACAGGAGTGtttca                                                                                         >  1:3838745/1‑96 (MQ=255)

              aatTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttcaaa                                                                                      <  1:67445/100‑5 (MQ=255)

               attattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCg                                                                                     >  2:3458582/1‑100 (MQ=255)

               agtattGGAGACAGGACTGGGAAGATATAAGAAATGTGGACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttcaa                                                                                       <  1:3883408/96‑4 (MQ=255)

                ttattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCgg                                                                                    <  2:3831106/100‑1 (MQ=255)

                ttattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCgg                                                                                    <  2:2152712/100‑1 (MQ=255)

                ttattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCgg                                                                                    <  2:2800406/100‑1 (MQ=255)

                ttattTGAAACAGGCCTGGAACGATAGAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCgg                                                                                    >  2:2684374/1‑100 (MQ=255)

                ttattTGAAAAAAGCATCGAACAATTTAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCTCGCTTTTTGAAAACACACTGGTGTGTTTTATGTgg                                                                                    <  1:1914144/100‑1 (MQ=255)

                 tattTGAAACAGTACTGGAACAATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttcaaaaag                                                                                   <  1:893180/100‑8 (MQ=37)

                 tattTGAAACAGGCCTGGAAGGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGt                                                                                   <  1:1609205/100‑1 (MQ=255)

                 tattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGt                                                                                   >  2:3229623/1‑100 (MQ=255)

                 tattTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGt                                                                                   <  1:45773/100‑1 (MQ=255)

                 tattTGAAACAGGCCTCGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCAGTGTGttttcaaaaag                                                                                   <  1:2585228/100‑8 (MQ=38)

                  attTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGtt                                                                                  >  2:3140205/1‑100 (MQ=255)

                   tttGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTta                                                                                 >  1:3528629/1‑100 (MQ=255)

                   tttGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTta                                                                                 >  2:1342968/1‑100 (MQ=255)

                    ttGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTtat                                                                                <  2:1897813/100‑1 (MQ=255)

                      gAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATAc                                                                              <  2:860310/100‑1 (MQ=255)

                          caagCCTTGAACAATATAAAAAATGAGTACAAAAAATATTAAATGGTTGCGATCTGTCTTTTTTAAAACACATTGTTGTGTTTTATGCGGTTAGACAATg                                                                          <  1:1666993/97‑1 (MQ=255)

                           aGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGa                                                                         >  1:4358024/1‑100 (MQ=255)

                            ggCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGaa                                                                        >  1:4101904/1‑100 (MQ=255)

                            ggCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACAGGTGTGTTTTATGCGGTTATACAATGaa                                                                        >  1:2083948/1‑100 (MQ=255)

                             gCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAAc                                                                       <  2:3165363/100‑1 (MQ=255)

                               cTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGGACTc                                                                     >  2:1355114/1‑100 (MQ=255)

                                 ggAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAg                                                                   >  2:849142/1‑100 (MQ=255)

                                   aaCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGAc                                                                 <  2:1749205/100‑1 (MQ=255)

                                    aCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACg                                                                <  2:1659461/100‑1 (MQ=255)

                                    aCGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACg                                                                <  1:1392360/100‑1 (MQ=255)

                                      gATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACCGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgta                                                              <  1:2713808/100‑1 (MQ=255)

                                       atatAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaa                                                             <  2:531841/100‑1 (MQ=255)

                                         ataAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAATACACTGATGTGTATTATGCGGTTATACAATGAACTAAGACgtaaag                                                           >  2:404981/1‑100 (MQ=255)

                                         ataAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaag                                                           >  1:2379058/1‑100 (MQ=255)

                                         ataAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaag                                                           <  1:3213457/100‑1 (MQ=255)

                                          taAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaagt                                                          <  2:4298404/100‑1 (MQ=255)

                                           aaaaaaTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaagta                                                         >  2:2638326/1‑100 (MQ=255)

                                            aaaaaTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACgtaaagtaa                                                        >  2:1732752/1‑100 (MQ=255)

                                               aaTGAGTACGAAAAAAATCAAACCTATTCAATCACGGTTGTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATg                                                     <  1:3457501/100‑1 (MQ=255)

                                               aaTGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATg                                                     <  2:351508/100‑1 (MQ=255)

                                                 tGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGaa                                                   >  2:4167332/1‑100 (MQ=255)

                                                 tGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGaa                                                   <  2:596627/100‑1 (MQ=255)

                                                   aGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAaga                                                 <  2:2039345/100‑1 (MQ=255)

                                                    gTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAagag                                                <  2:403780/100‑1 (MQ=255)

                                                      aCGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGa                                              <  2:1600607/100‑1 (MQ=255)

                                                         aaaaaaaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATAc                                           >  1:179075/1‑100 (MQ=255)

                                                         aaaaaaaTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATAc                                           >  2:2324029/1‑100 (MQ=255)

                                                         aaaaaaaTCAAACATATTCAATCACGCTTTTTGAAAAAACAATGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATa                                            >  1:2462263/1‑99 (MQ=255)

                                                               aTCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATg                                     >  2:3016925/1‑100 (MQ=255)

                                                                      atatTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGaa                              <  2:619819/100‑1 (MQ=255)

                                                                       tatTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAt                             >  1:3484531/1‑100 (MQ=255)

                                                                       tatTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAt                             <  2:1270671/100‑1 (MQ=255)

                                                                         ttCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTATATGAAGAGGATACTCCATGTCTGGAAtgt                           >  1:3973320/1‑100 (MQ=255)

                                                                         ttCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgt                           >  1:10405/1‑100 (MQ=255)

                                                                         ttCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgt                           >  2:3336678/1‑100 (MQ=255)

                                                                         ttCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAAtgt                           <  1:782090/100‑1 (MQ=255)

                                                                               cACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCa                     >  2:2303521/1‑100 (MQ=255)

                                                                               cACGCCTTTTGAAAGCACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCa                     >  2:2494982/1‑100 (MQ=255)

                                                                                ataaccgcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCAc                    >  2:3342356/11‑100 (MQ=38)

                                                                                 cGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACt                   <  2:3728349/100‑1 (MQ=255)

                                                                                 cGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACt                   >  2:18361/1‑100 (MQ=255)

                                                                                 cGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACt                   <  2:2919642/100‑1 (MQ=255)

                                                                                  gCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTa                  >  1:4017798/1‑100 (MQ=255)

                                                                                  gCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTa                  <  1:131038/100‑1 (MQ=255)

                                                                                   accgcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTAAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAg                 >  1:4031838/8‑100 (MQ=38)

                                                                                    tttttGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGa                >  1:1530433/1‑100 (MQ=255)

                                                                                     cgcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaa               >  2:2217291/6‑100 (MQ=255)

                                                                                      gcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaa              >  1:992501/5‑100 (MQ=255)

                                                                                      gcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaa              >  1:1238446/5‑100 (MQ=255)

                                                                                      gcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaa              >  2:90317/5‑100 (MQ=255)

                                                                                      gcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaa               >  2:1117642/5‑99 (MQ=255)

                                                                                      gcataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGAATACTACATGTCTGGAATGTGAATAACTAGaaa              >  1:3855643/5‑100 (MQ=255)

                                                                                       ttGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa             <  1:1343845/100‑1 (MQ=255)

                                                                                       cataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa             >  1:697824/4‑100 (MQ=255)

                                                                                       cataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa             >  2:3277411/4‑100 (MQ=255)

                                                                                       cataaaaCACACCAGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa             >  1:3115222/4‑100 (MQ=255)

                                                                                       cataaaaCACACCAGTGTGTTTTAGACGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa             >  1:2168511/4‑100 (MQ=255)

                                                                                        tGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaaa            >  1:1188622/1‑100 (MQ=255)

                                                                                         gAAAAAACACTAGTGTGTTTTATGCGGTTATACAAAGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAg                 <  1:3677733/94‑1 (MQ=255)

                                                                                           aaaCACACTGGTGTGTTTTATGCGGTTATCCAATGAACTCAGACGTAAAGTACATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGaa         <  1:4357021/100‑1 (MQ=255)

                                                                                           aaaCACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGaaaa             >  1:3656170/1‑96 (MQ=255)

                                                                                             acacacGGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACCCCATGTCTGGAATGTGAATCACTAGAAAAAGAAt        <  1:2374386/99‑1 (MQ=255)

                                                                                              cacacTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACGCCATGTCTGGAATGTGAATCACTAGAAAAAGAAttt      >  2:3091268/1‑100 (MQ=255)

                                                                                                cacTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGa    >  1:2315087/1‑100 (MQ=255)

                                                                                                 acTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAg   <  1:1748515/100‑1 (MQ=255)

                                                                                                 acTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAg   <  2:607784/100‑1 (MQ=255)

                                                                                                  cTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAgg  >  2:1634583/1‑100 (MQ=255)

GTTGCACTGCAAATAATTATTTGAAACAGGCCTGGAACGATATAAAAAATGAGTACGAAAAAAATCAAACATATTCAATCACGCTTTTTGAAAACACACTGGTGTGTTTTATGCGGTTATACAATGAACTCAGACGTAAAGTAAATGAAGAGGATACTCCATGTCTGGAATGTGAATCACTAGAAAAAGAATTTGAGG  >  NC_000913/1064061‑1064258

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