BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,366,764	(C)_5→6	pseudogene (10/175 nt)	yhcE →	putative uncharacterized protein YhcE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,366,759	1	.	C	99.7%	1253.3 / ‑3.5	303	pseudogene (5/175 nt)	yhcE	putative uncharacterized protein YhcE
Reads supporting (aligned to +/- strand): ref base . (1/0); new base C (158/144); total (159/144)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 339 total aligned reads displayed.

GGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTT‑CCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACAATGGCTCCAGCT  >  NC_000913/3366659‑3366854

ggCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACtt                                                                                                  >  1:6133745/1‑101 (MQ=32)

ggCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACtt                                                                                                  >  2:6133749/1‑101 (MQ=32)

  cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTT‑cc                                                                                               >  1:6133753/1‑101 (MQ=35)

   gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTT‑ccc                                                                                              <  2:6133296/101‑1 (MQ=35)

     aGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTT‑ccccc                                                                                            >  1:6133757/1‑101 (MQ=38)

       gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTT‑cccccca                                                                                          <  1:6133082/101‑3 (MQ=38)

        tAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTT‑ccccccaa                                                                                         >  2:6133761/1‑98 (MQ=255)

        tAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTT‑ccccccaa                                                                                         <  2:6133548/101‑4 (MQ=255)

         aaGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTT‑ccccccaat                                                                                        <  2:6133348/101‑5 (MQ=255)

         aaGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTT‑ccccccaat                                                                                        >  1:6133765/1‑97 (MQ=255)

          aGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTT‑ccccccaata                                                                                       <  2:6133069/101‑6 (MQ=255)

           gTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCaataa                                                                                       <  1:6133405/101‑1 (MQ=255)

             tgatgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAAtt                                                                                     <  1:6132908/101‑1 (MQ=255)

               atgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTgg                                                                                   <  2:6133142/101‑1 (MQ=255)

               atgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTgg                                                                                   >  1:6133772/1‑101 (MQ=255)

               atgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTgg                                                                                   <  2:6133216/101‑1 (MQ=255)

               atgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTgg                                                                                   >  2:6133771/1‑101 (MQ=255)

                tgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGc                                                                                  >  1:6133777/1‑101 (MQ=255)

                  aCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTg                                                                                <  2:6133149/101‑1 (MQ=255)

                  aCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTg                                                                                <  2:6133186/101‑1 (MQ=255)

                  aCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTg                                                                                >  2:6133781/1‑101 (MQ=255)

                   cTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTgg                                                                               <  1:6133034/101‑1 (MQ=255)

                      atgatgAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGgtgt                                                                            >  1:6133784/1‑101 (MQ=255)

                         atgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTa                                                                         >  2:6133787/1‑101 (MQ=255)

                         atgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTa                                                                         <  2:6133092/101‑1 (MQ=255)

                          tgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTAt                                                                        <  1:6133312/101‑1 (MQ=255)

                             aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCa                                                                     >  2:6133791/1‑101 (MQ=255)

                              cccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCaa                                                                    >  1:6133794/1‑101 (MQ=255)

                               ccTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCaaa                                                                   <  1:6133328/101‑1 (MQ=255)

                                cTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCaaaa                                                                  >  2:6133798/1‑101 (MQ=255)

                                 tGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAAt                                                                 <  2:6133247/101‑1 (MQ=255)

                                 tGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAAt                                                                 >  1:6133799/1‑101 (MQ=255)

                                  gTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCa                                                                     <  2:6133285/96‑1 (MQ=255)

                                   ttCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATgg                                                               <  2:6133490/101‑1 (MQ=255)

                                    tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGt                                                              >  1:6133804/1‑101 (MQ=255)

                                     cTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTg                                                                                >  2:6133807/1‑82 (MQ=255)

                                      tATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCa                                                            <  1:6133138/101‑1 (MQ=255)

                                      tATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCa                                                            <  2:6133107/101‑1 (MQ=255)

                                       aTGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCaa                                                           <  2:6133128/101‑1 (MQ=255)

                                        tGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCaag                                                          <  2:6133230/101‑1 (MQ=255)

                                          gCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCaagaa                                                        >  1:6133812/1‑101 (MQ=255)

                                           cTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCaagaag                                                       >  2:6133813/1‑101 (MQ=255)

                                            tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGt                                                      >  2:6133816/1‑101 (MQ=255)

                                            tCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGt                                                      <  1:6133159/101‑1 (MQ=255)

                                             ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTg                                                     <  2:6133359/101‑1 (MQ=255)

                                             ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTg                                                     <  2:6133238/101‑1 (MQ=255)

                                              cAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTgg                                                    <  2:6133634/101‑1 (MQ=255)

                                                gATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTa                                                  >  2:6133824/1‑101 (MQ=255)

                                                 aTGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTaa                                                 >  1:6133826/1‑101 (MQ=255)

                                                  tGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAAt                                                <  2:6133204/101‑1 (MQ=255)

                                                   gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAAtt                                               >  2:6133829/1‑101 (MQ=255)

                                                   gACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAAtt                                               <  2:6133178/101‑1 (MQ=255)

                                                     cAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGa                                             <  2:6133221/101‑1 (MQ=255)

                                                     cAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGa                                             >  1:6133835/1‑101 (MQ=255)

                                                      aaaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGAt                                            <  1:6133353/101‑1 (MQ=255)

                                                       aaCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATa                                           >  2:6133837/1‑101 (MQ=255)

                                                        aCATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAg                                          >  1:6133840/1‑101 (MQ=255)

                                                         cATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGc                                         >  1:6133845/1‑101 (MQ=255)

                                                          aTGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCg                                        >  1:6133849/1‑101 (MQ=255)

                                                          aTGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCg                                        >  2:6133847/1‑101 (MQ=255)

                                                            gATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGt                                      <  1:6133299/101‑1 (MQ=255)

                                                              tctcATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTaa                                    <  1:6133606/101‑1 (MQ=255)

                                                              tctcATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTaa                                    >  2:6133854/1‑101 (MQ=255)

                                                               ctcATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTaaa                                   >  1:6133856/1‑101 (MQ=255)

                                                                tcATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTaaaa                                  >  1:6133860/1‑101 (MQ=255)

                                                                tcATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTaaaa                                  <  2:6133427/101‑1 (MQ=255)

                                                                 cATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAg                                 >  2:6133863/1‑101 (MQ=255)

                                                                 cATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAg                                 <  2:6133262/101‑1 (MQ=255)

                                                                 cATATCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAg                                 <  2:6133317/101‑1 (MQ=255)

                                                                  atatCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGc                                <  1:6133485/101‑1 (MQ=255)

                                                                   tatCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCg                               <  1:6133372/101‑1 (MQ=255)

                                                                    atCAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGt                              >  2:6133866/1‑101 (MQ=255)

                                                                      cAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTa                            <  1:6133222/101‑1 (MQ=255)

                                                                      cAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTa                            <  1:6133170/101‑1 (MQ=255)

                                                                      cAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTa                            >  2:6133870/1‑101 (MQ=255)

                                                                      cAGGGACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTa                            <  2:6133330/101‑1 (MQ=255)

                                                                        gggACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGc                          >  2:6133875/1‑101 (MQ=255)

                                                                         ggACTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGcc                         >  1:6133876/1‑101 (MQ=255)

                                                                            cTTGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTca                      >  2:6133881/1‑101 (MQ=255)

                                                                             ttGTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTcac                     >  1:6133885/1‑101 (MQ=255)

                                                                               gTTCGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCacac                   >  2:6133888/1‑101 (MQ=255)

                                                                                ttggCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACt                  >  1:6133892/4‑101 (MQ=255)

                                                                                  cGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGa                >  1:6133896/1‑101 (MQ=255)

                                                                                  cGCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGa                >  2:6133894/1‑101 (MQ=255)

                                                                                   gCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGAc               <  2:6133460/101‑1 (MQ=255)

                                                                                   gCACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGAc               >  2:6133901/1‑101 (MQ=255)

                                                                                    cACCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACa              <  1:6133471/101‑1 (MQ=255)

                                                                                     aCCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACaa             >  2:6133906/1‑101 (MQ=255)

                                                                                     aCCTTCCTTAGTACTTCCCCCCAATAATTGGCTGGGGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCt                        >  1:6133907/1‑90 (MQ=255)

                                                                                      ccttccttAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTaaaa                                  <  2:6133378/79‑1 (MQ=255)

                                                                                       cttccttAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACAATg           <  1:6133270/101‑1 (MQ=255)

                                                                                         tccttAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACAATGGc         <  1:6133207/101‑1 (MQ=255)

                                                                                          ccttAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACAATGGCt        >  1:6133912/1‑101 (MQ=255)

                                                                                          ccttAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACAATGGCt        >  2:6133913/1‑101 (MQ=255)

                                                                                            ttAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACAATGGCTcc      <  1:6133444/101‑1 (MQ=255)

                                                                                            ttAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACAATGGCTcc      >  2:6133919/1‑101 (MQ=255)

                                                                                             tAGTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACAATGGCTCCa     >  2:6133923/1‑101 (MQ=255)

                                                                                               gTACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACAATGGCTCCAGc   <  1:6133820/101‑1 (MQ=255)

                                                                                                tACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACAATGGCTCCAGCt  >  1:6133930/1‑101 (MQ=255)

                                                                                                tACTTCCCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACAATGGCTCCAGCt  >  1:6133926/1‑101 (MQ=255)

GGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGTACTT‑CCCCCAATAATTGGCTGGTGTTTATGCAAAATGGTCAAGAAGTGGTAATTGATAGCGGTAAAAGCGTTAGCCTCACACTGACAATGGCTCCAGCT  >  NC_000913/3366659‑3366854

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