BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,980,297	(C)_8→9	noncoding (79/87 nt)	leuT →	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,980,289	1	.	C	76.6%	633.8 / 196.6	257	noncoding (71/87 nt)	leuT	tRNA‑Leu
Reads supporting (aligned to +/- strand): ref base . (22/38); new base C (87/109); total (110/147)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 325 total aligned reads displayed.

CATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGT  >  NC_000913/3980197‑3980380

cATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑t                                                                                                >  1:1836051/1‑93 (MQ=255)

cATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGCTCAAG‑‑t                                                                                                >  2:4350850/1‑93 (MQ=255)

 attattAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑c                                                                                            >  1:3370559/1‑93 (MQ=255)

  ttattaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑cc                                                                                           >  1:1692088/1‑93 (MQ=255)

  ttattaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑cc                                                                                           >  2:7374123/1‑93 (MQ=255)

  ttattaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑cc                                                                                           >  1:4107854/1‑93 (MQ=255)

   tattaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑ccc                                                                                          <  1:2421794/93‑1 (MQ=255)

    attaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑cccc                                                                                         >  2:953981/1‑93 (MQ=255)

    attaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑cccc                                                                                         >  1:7193073/1‑93 (MQ=255)

    attaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑cccc                                                                                         <  1:6682191/93‑1 (MQ=255)

    attaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑cccc                                                                                         <  1:891915/93‑1 (MQ=255)

    attaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑cccc                                                                                         <  2:6584650/93‑1 (MQ=255)

    attaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑cccc                                                                                         <  2:7420165/93‑1 (MQ=255)

     ttaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑ccccc                                                                                        <  2:6956144/93‑1 (MQ=255)

     ttaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑ccccc                                                                                        >  1:5961716/1‑93 (MQ=255)

     ttaGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAG‑‑T‑‑‑ccccc                                                                                        >  1:3757186/1‑93 (MQ=255)

       aGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑ccccccc                                                                                      >  2:3983642/1‑93 (MQ=255)

        gAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑cccccccc                                                                                     <  1:4789007/93‑1 (MQ=255)

          aGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑ccccccccct                                                                                   >  2:785473/1‑91 (MQ=255)

          aGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑ccccccccct                                                                                   <  2:6678037/93‑3 (MQ=255)

          aGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑cccccccccc                                                                                   <  2:1863938/93‑3 (MQ=255)

           gTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑cccccccccct                                                                                  <  1:4596069/93‑4 (MQ=255)

           gTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑ccccccccccc                                                                                  >  1:4303829/1‑90 (MQ=255)

            ttGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTTAAG‑‑T‑‑‑cccccccccccc                                                                                 >  1:1309969/1‑89 (MQ=255)

            ttGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑cccccccccctc                                                                                 <  1:2829554/93‑5 (MQ=255)

               tgACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGcacc                                                                              >  1:6445580/1‑93 (MQ=35)

                gACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGcacc                                                                              >  1:275398/1‑93 (MQ=255)

                  cAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCC‑CCCCCCCCTCGcacca                                                                             >  1:4418601/1‑93 (MQ=11)

                     tGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGAc                                                                         >  1:2197375/1‑93 (MQ=255)

                     tGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGAc                                                                         <  1:7442142/93‑1 (MQ=255)

                      gCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACt                                                                        <  2:6174157/93‑1 (MQ=255)

                      gCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACt                                                                        <  2:6848277/93‑1 (MQ=255)

                      gCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACt                                                                        <  2:7077211/93‑1 (MQ=255)

                        gAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTa                                                                     <  2:719655/93‑1 (MQ=39)

                         aaGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTa                                                                     >  2:6985161/1‑93 (MQ=255)

                          aGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTaa                                                                    >  2:4072023/1‑93 (MQ=255)

                          aGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTaa                                                                    <  2:2403936/93‑1 (MQ=255)

                          aGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTaaa                                                                   <  1:3325714/93‑1 (MQ=255)

                          aGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACCTTaaa                                                                   >  2:5339336/1‑93 (MQ=255)

                           ggTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTaaa                                                                   <  1:7215585/93‑1 (MQ=255)

                           ggTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAg                                                                  <  2:1578265/93‑1 (MQ=255)

                            gTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAg                                                                  <  2:1249403/93‑1 (MQ=255)

                            gTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAg                                                                  <  2:1668564/93‑1 (MQ=255)

                             tGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGa                                                                 <  2:6777725/93‑1 (MQ=255)

                              ggCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGaa                                                                <  2:1355331/93‑1 (MQ=255)

                              ggCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAAt                                                               <  1:1520974/93‑1 (MQ=255)

                              ggCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAAt                                                               <  1:1918578/93‑1 (MQ=255)

                                 ggAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTg                                                             >  1:6086605/1‑93 (MQ=255)

                                   aaTTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGaa                                                           <  2:7014383/93‑1 (MQ=255)

                                     ttGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACt                                                         >  1:5394210/1‑93 (MQ=255)

                                      tGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTa                                                        <  2:3420358/93‑1 (MQ=255)

                                      tGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAAAATTGAACTa                                                        >  1:5100208/1‑93 (MQ=255)

                                        gTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTaaa                                                      <  1:4932157/93‑1 (MQ=255)

                                         tAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTaaaa                                                     >  2:2805470/1‑93 (MQ=255)

                                         tAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTaaaa                                                     >  2:4606060/1‑93 (MQ=255)

                                          aGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTaaaaa                                                    <  1:1102885/93‑1 (MQ=255)

                                          aGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTaaaaa                                                    <  1:3556344/93‑1 (MQ=255)

                                          aGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAG‑‑TC‑‑CCCCCCCCCCCCACCACCACTTTTAAAAAATGAAATaaaaa                                                    >  1:6340492/1‑93 (MQ=14)

                                           gACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTTAAAAATTGAAACaaaaaa                                                   >  1:1755516/1‑92 (MQ=14)

                                           gACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAAt                                                   >  2:5826580/1‑93 (MQ=255)

                                           gACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAAt                                                   <  2:5945641/93‑1 (MQ=255)

                                           gACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAAt                                                   >  2:3106129/1‑93 (MQ=255)

                                           gACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAAtt                                                  >  1:2273795/1‑93 (MQ=255)

                                            aCGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTc                                                 <  2:6392256/93‑1 (MQ=255)

                                              gcgcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCa                                                >  1:5583138/1‑93 (MQ=255)

                                               cgcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaa                                              <  2:1146355/93‑1 (MQ=255)

                                                gcTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGGTCAAG‑‑TC‑‑CCCCCCCCTCGGACCACCACTTTAAAGAATTGAACTAAAAATTCaaa                                              >  1:4913364/1‑93 (MQ=21)

                                                  tAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaaaaa                                            <  1:6514355/93‑1 (MQ=255)

                                                  tAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAg                                           <  2:2358504/93‑1 (MQ=255)

                                                    gCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGc                                          >  1:3514956/1‑93 (MQ=255)

                                                    gCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGc                                          <  2:4877051/93‑1 (MQ=255)

                                                     cTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCa                                         <  1:2966216/93‑1 (MQ=255)

                                                        cAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTa                                      <  2:898044/93‑1 (MQ=255)

                                                         aGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAt                                     <  1:7354494/93‑1 (MQ=255)

                                                         aGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAt                                     >  2:5763890/1‑93 (MQ=255)

                                                         aGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAtt                                    >  2:5514918/1‑93 (MQ=255)

                                                          ggTGTTAGTGTCCTTACGGACGGGGGGGTTCAAGGCC‑‑‑CCCCCCCCTCGCTCTATGGCTTTAAAGAATGGAACTAAAAATTCAAAAAGCAGTAt                                     >  2:7294071/1‑93 (MQ=18)

                                                           gtgtTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTAttt                                   <  2:6344945/93‑1 (MQ=255)

                                                             gtTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCg                                 <  1:786089/93‑1 (MQ=255)

                                                              ttAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCgg                                <  2:6510106/93‑1 (MQ=255)

                                                                aGTGTCCTTACGGACGTGGGGGTTCAAG‑‑TCCCCCCCCCCCCCCCACCCCCCCTTTAAAAAATTTAAATAAAAATTCAAAAAACAATATTTCgg                                >  2:3212568/1‑93 (MQ=9)

                                                                aGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGa                             <  2:5130452/93‑1 (MQ=255)

                                                                 gtgtCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAg                            >  1:1033030/1‑93 (MQ=255)

                                                                  tgtCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAg                            <  1:6206080/93‑1 (MQ=255)

                                                                    tCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTa                          <  2:4921318/93‑1 (MQ=255)

                                                                    tCCTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTa                          <  2:3344077/93‑1 (MQ=255)

                                                                     ccTTACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAg                         <  1:2569437/93‑1 (MQ=255)

                                                                       ttACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgcg                       >  2:2663039/1‑93 (MQ=255)

                                                                        tACGGACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAgcgc                      <  2:4826273/93‑1 (MQ=255)

                                                                           ggACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGc                   <  2:4208290/93‑1 (MQ=255)

                                                                           ggACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGc                   >  2:5711438/1‑93 (MQ=255)

                                                                            gACGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCt                  <  2:2477029/93‑1 (MQ=255)

                                                                             aCGTGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCtt                 <  1:5683965/93‑1 (MQ=255)

                                                                                tGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGt              >  2:3744630/1‑93 (MQ=255)

                                                                                tGGGGGTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGt              <  1:986367/93‑1 (MQ=255)

                                                                                 gggggTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTa             >  1:6892962/1‑93 (MQ=255)

                                                                                 gggggTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAg            >  2:2894553/1‑93 (MQ=255)

                                                                                     gTTCAAG‑‑TC‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAgcgc         <  2:5013024/93‑1 (MQ=255)

                                                                                      ttCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCaa       >  2:209874/1‑93 (MQ=255)

                                                                                          aG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTgg   >  1:3090398/1‑93 (MQ=255)

CATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAG‑‑T‑‑‑CCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCGGCGAGTAGCGCAGCTTGGTAGCGCAACTGGT  >  NC_000913/3980197‑3980380

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