BRESEQ :: Evidence

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,781,152	0	A	G	21.1%	259.8 / 20.7	136	intergenic (+135/‑63)	lldD/trmL	L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand): ref base A (71/34); new base G (0/28); total (71/65)
Fisher's exact test for biased strand distribution p-value = 7.75e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 141 total aligned reads displayed.

GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCG  >  NC_000913/3781054‑3781250

ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACACTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑Ga                                                                                                     >  1:150279/1‑99 (MQ=255)

 gagaGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAg                                                                                                    >  2:174071/1‑99 (MQ=255)

  agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAgg                                                                                                   <  1:4207791/99‑1 (MQ=255)

  agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAgg                                                                                                   <  1:3394010/99‑1 (MQ=255)

  agagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTC‑TTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAggg                                                                                                  >  1:3797356/1‑99 (MQ=255)

   gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAggg                                                                                                  <  2:2508457/99‑1 (MQ=255)

   gagGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAggg                                                                                                  <  2:3156264/99‑1 (MQ=255)

        gACCGTTCGGCGCTGTATGTTCTCCCTCTCTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGa‑aaa                                                                                            <  1:1236933/99‑1 (MQ=255)

        gACCGTTCGGCGCTGTATGTACTCTCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGa‑aaa                                                                                            <  1:629028/99‑1 (MQ=255)

        gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑TGGGGCTGGCGT‑GAGGGGa‑aaa                                                                                            >  2:388322/1‑99 (MQ=255)

        gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGa‑aaa                                                                                            <  1:3936959/99‑1 (MQ=255)

        gACCGTTCGGCGCTGTATAGACTCCCCAACTCTGAAACGACACCGCACTCTTTGTTACTC‑CCTCGCCCCACCGGGGA‑GAGGCCCGGGGC‑GAAGGCa‑aaa                                                                                            >  2:5161518/1‑99 (MQ=255)

         aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAg                                                                                           >  2:4593734/1‑99 (MQ=255)

          ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAgg                                                                                          <  1:2045110/99‑1 (MQ=255)

           cGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGc                                                                                         >  1:2266378/1‑99 (MQ=255)

             ttCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCg                                                                                       >  2:3422580/1‑99 (MQ=255)

               cGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGC‑AAAGGCCGCa                                                                                     >  2:77684/1‑99 (MQ=255)

                ggCGCTGTATGTACTCCCTCACTCTGAAACGACCCCGCACTCTTTTTATCAC‑TCTCGCCCCACCGGGGA‑GAGGACCGGGGT‑GAGGGGA‑AAAGGCCGCaa                                                                                    >  1:4105465/1‑98 (MQ=255)

                ggCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCAc                                                                                    >  2:1878597/1‑99 (MQ=255)

                ggCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCAc                                                                                    >  2:1985991/1‑99 (MQ=255)

                 gcgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACt                                                                                   >  2:4088229/1‑99 (MQ=255)

                 gcgcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACt                                                                                   >  1:3240388/1‑99 (MQ=255)

                   gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGc                                                                                 >  1:3277549/1‑99 (MQ=255)

                       tatgtaCCCCCTCACCCTGAAACGACCCCCCCCTTTTTTTTTCTC‑CCCCGCCCCGCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAGGGCCGCACTGCTCAc                                                                             <  2:2650455/99‑1 (MQ=25)

                           taCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAAt                                                                         >  2:524527/1‑99 (MQ=255)

                            atccccTCATCCTAAAAAGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCCCGGGGGA‑GAGGCCCGGGGT‑GAGGGAA‑TAAGGCCGCACTGCTCACCAAtt                                                                        <  2:339155/96‑1 (MQ=255)

                            aCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAAtt                                                                        >  1:2735316/1‑99 (MQ=255)

                              tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAAGGCCGGGGA‑GAGGGGA‑AAAGGCCGCACTGCTCACCAAtttt                                                                      >  2:3819720/1‑99 (MQ=255)

                               ccctccctCTGAAAGGACCCCGCACTTTTTTTTTCTC‑CCTCGCCCCCCCGGGGG‑GGGGGCGGGGGT‑GGGGGGA‑AAAGGCCGCACTGCTCACTAATTTTc                                                                     <  2:4960728/99‑1 (MQ=255)

                                ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCt                                                                    >  1:1841788/1‑99 (MQ=255)

                                ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGA‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCt                                                                    >  1:1689491/1‑99 (MQ=255)

                                ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCACCTCCGGGGA‑GAGGGCCGGGGT‑GAGGTGA‑AAAGGCCGGTCTGCTCACCAATTTTCt                                                                    >  2:3124861/1‑99 (MQ=38)

                                 ccccctctGAACCGATTCCGCTTTTTTTTTTTTCC‑CTTCGCCCTTCGGGGGA‑GGGGGCGGGGGT‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTg                                                                   <  1:3705709/95‑1 (MQ=255)

                                 cccACTCGGAAACGACACCGCACTTTTTTTTTTTC‑CCTCGCCCCCCGGGGGA‑GAGGGCCGGGGG‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTg                                                                   <  2:1302262/97‑1 (MQ=255)

                                 cTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTg                                                                   >  1:1152098/1‑99 (MQ=255)

                                 cTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTg                                                                   >  2:119342/1‑99 (MQ=255)

                                 cTCACACTGAACCGACCCCGCCTTTTTTTTTTCTC‑CCCCGCCCCTCCGGGGG‑GAGGGCGGGGGT‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTg                                                                   <  2:1805896/99‑1 (MQ=255)

                                  tCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGA‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGc                                                                  >  1:39024/1‑99 (MQ=255)

                                  cccctctGAAACGACCCCCCCCTC‑TTTTTTCTC‑CCTCGCCCCCCCGGGGA‑GAGGGCCGGGGG‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCt                                                                 <  2:3740556/96‑1 (MQ=255)

                                     ctctGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTAt                                                               >  1:4748871/1‑99 (MQ=255)

                                      tctGAAACGACCCCTCTTTTTTTTTTTTTC‑CCTCCCCCTTCGGGGGA‑GGGGGCGGGGGG‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                              <  1:947847/99‑1 (MQ=255)

                                          aaaCGACACCGCACTCTTTTTTTCTT‑TTTTTCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                                          >  2:1276478/1‑99 (MQ=38)

                                          aaaCGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                                          >  2:4994760/1‑99 (MQ=255)

                                           aaCGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcc                                                         >  2:908388/1‑99 (MQ=255)

                                            aCGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccc                                                        >  1:3337177/1‑99 (MQ=255)

                                             cGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑GAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGcccc                                                       >  1:4944582/1‑99 (MQ=255)

                                                 cccGCACTCTTTTTTCCTC‑CCGCCCCCCTCCGGGGA‑GAGGGCCGGGGT‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCa                                                   <  2:3348241/98‑1 (MQ=255)

                                                     cACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCACCCGCACTaa                                               >  1:1975508/1‑99 (MQ=255)

                                                        tttttttttCTC‑CTCCCCCCCTCGGGGAG‑GAGGGCGGGGGT‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                            <  1:5192431/97‑1 (MQ=255)

                                                        tctTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                            >  2:407096/1‑99 (MQ=255)

                                                         tttttcttCTC‑CCTCGCCCCCCCGGGGG‑GAGGGCCGGGGT‑GGGGGCA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                           <  2:1308383/98‑1 (MQ=255)

                                                          tttttttCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                          >  1:2316884/1‑99 (MQ=255)

                                                          tttttttCCCTCCCCCCCCCTCGGGGTA‑GAGGGCGTGGGT‑GAGGCGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                           <  2:1846296/99‑1 (MQ=255)

                                                          tttttttC‑C‑CCGAGGCCCCCCGGGGAGGAGGGCGGGGGG‑GGGGGAA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                          <  2:683832/99‑1 (MQ=255)

                                                          tttttttC‑C‑CCCCCCCCCTTCGGGGG‑GGGGGGCGGGGTGGGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                          <  1:5058558/99‑1 (MQ=255)

                                                            tttttCTC‑CTTCCCCCTTCGGGGGG‑GGGGGCGGGGGT‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCa                                        <  1:2578030/99‑1 (MQ=255)

                                                            tttttCTC‑CCTCGCCCCTCCGGGGA‑GGGGGCGGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCa                                        <  1:3368859/99‑1 (MQ=255)

                                                            tttttCTC‑CCTCCCCCCTCGGGGGA‑GGGGGCCGGGGT‑GAGGGTA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCa                                        <  1:151544/99‑1 (MQ=39)

                                                                tccc‑ccGCGCCCCTCGGGGGG‑GGGGGCGGGGGT‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                    <  1:4010487/96‑1 (MQ=255)

                                                                cccc‑cccCGCCCCTCGGGGGG‑GAGGGCGGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAt                                    <  2:2075427/96‑1 (MQ=38)

                                                                     ccTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                                >  2:826821/1‑99 (MQ=255)

                                                                     ccTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                                >  2:2419280/1‑99 (MQ=255)

                                                                     ccTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTa                                >  1:4494148/1‑99 (MQ=255)

                                                                      cTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTAACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaa                               >  1:638843/1‑99 (MQ=255)

                                                                       tCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                              >  2:1439822/1‑99 (MQ=255)

                                                                       tCGCCCCTCCGGGGA‑GAGGGCCGAGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                              <  2:2539806/99‑1 (MQ=255)

                                                                       gcGCCCCCCCGGGGG‑GGGGGTGGGGGT‑GTGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                              <  1:4107487/98‑1 (MQ=255)

                                                                        cGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAAc                             >  1:5246193/1‑99 (MQ=255)

                                                                           cccTCGGGGGG‑GGGGGCTGTGGT‑GTGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                          <  1:4411406/99‑1 (MQ=255)

                                                                           cccTCGGGGGG‑GAGGGCCGGGGG‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                          <  2:4646869/99‑1 (MQ=255)

                                                                           cccTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                          >  1:2523856/1‑99 (MQ=255)

                                                                           cccTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                          >  2:3654095/1‑99 (MQ=255)

                                                                           cccTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATc                          >  1:1010097/1‑99 (MQ=255)

                                                                            ccTCGGGGGG‑GGGGGCCGGGTT‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                         <  1:3503031/99‑1 (MQ=255)

                                                                             cTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                        >  2:4404136/1‑99 (MQ=255)

                                                                             cTCCGGGGA‑GAGGACCTAGGT‑GAGGGGA‑AAAGGTCGCACTGCTAACCAATTTTCTGCTATCCTGCCCCCGTACTTAGGGGGCAGTATGCTAAACATCGt                        >  2:4460250/1‑99 (MQ=255)

                                                                             atcgggggg‑gTGGGCGGGGGG‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGt                        <  1:3825563/95‑1 (MQ=255)

                                                                              tCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa                       >  1:2073755/1‑99 (MQ=255)

                                                                              tCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa                       >  1:1531700/1‑99 (MQ=255)

                                                                               tctgggA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTAc                      <  1:3632453/96‑1 (MQ=255)

                                                                                  gggg‑GAGGGCCGGGGT‑GGGTGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt                   <  2:4615619/99‑1 (MQ=255)

                                                                                  gggA‑GGGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt                   <  1:2973338/99‑1 (MQ=255)

                                                                                   gggcgagGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt                   <  1:972199/95‑1 (MQ=255)

                                                                                   ggA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa                  >  1:2554020/1‑99 (MQ=255)

                                                                                      ggggggCGGGGGG‑GGGGGCA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTGCTTTACg                <  2:1654989/96‑1 (MQ=255)

                                                                                       gggggCGGGGTT‑GTGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGa               <  1:1786429/97‑1 (MQ=255)

                                                                                        agGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGATATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGaa              >  1:95836/1‑99 (MQ=255)

                                                                                         gggCCGGTCT‑GGGGGGA‑GAAGGCCGTACTGCTCACCAATTTGCTGCTATCCGGCCCCCGCACTACGGGGGCAGTATGCTAAACATCGTACTTTACGAAc             <  2:309034/99‑1 (MQ=255)

                                                                                          gtgcGGGGG‑GGGGGGATAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAc             <  1:1596306/96‑1 (MQ=255)

                                                                                          ggCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc            >  2:292017/1‑99 (MQ=255)

                                                                                          ggCCGGGGT‑GAGGGAA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAcc            <  1:5011943/99‑1 (MQ=255)

                                                                                           gcggggtT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCa           <  2:1124313/96‑1 (MQ=255)

                                                                                           gCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCa           >  1:3241334/1‑99 (MQ=255)

                                                                                           ccgggggT‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCa           <  1:4695117/96‑1 (MQ=255)

                                                                                            cggggtT‑GGGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg          <  2:3685203/97‑1 (MQ=255)

                                                                                            ccGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg          >  2:3070739/1‑99 (MQ=255)

                                                                                              ggggT‑GAGGGGA‑GAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGaa        <  2:2501665/99‑1 (MQ=255)

                                                                                                    aaGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTcc   <  2:3225736/98‑1 (MQ=255)

                                                                                                     aGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCg  <  2:5138714/99‑1 (MQ=255)

                                                                                                     aGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCg  >  1:3131521/1‑99 (MQ=255)

GGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTC‑CCTCGCCCCTCCGGGGA‑GAGGGCCGGGGT‑GAGGGGA‑AAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCG  >  NC_000913/3781054‑3781250

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