BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,033,458	T→G	51.7%	I9I (ATA→ATC)	yedR ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,033,458	0	T	G	51.7%	22.7 / 232.1	156	I9I (ATA→ATC)	yedR	inner membrane protein
Reads supporting (aligned to +/- strand): ref base T (42/33); new base G (37/43); total (79/77)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

Only 100 of 163 total aligned reads displayed.

CAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGAATACA  >  NC_000913/2033362‑2033557

cAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTAtg                                                                                                  >  2:2147749/1‑100 (MQ=255)

 aGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTAtgt                                                                                                 <  2:1839040/100‑1 (MQ=255)

   tttGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGCAAATTCAGACTTAGTACTATGATGTGa                                                                                               >  1:2531329/1‑100 (MQ=255)

      gccgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATaa                                                                                            <  2:883836/100‑1 (MQ=255)

       ccgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATaaa                                                                                           >  2:1174788/1‑100 (MQ=255)

        cgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATaaaa                                                                                          >  1:1529878/1‑100 (MQ=255)

        cgccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATaaaa                                                                                          >  2:2460143/1‑100 (MQ=255)

         gccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAg                                                                                         >  2:3453254/1‑100 (MQ=255)

         gccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAg                                                                                         >  2:3887640/1‑100 (MQ=255)

         gccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACACGGAAAATTCAGACTTAGTACTATGATGTGATAAAAg                                                                                         <  1:923274/100‑1 (MQ=255)

          ccAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGt                                                                                        >  1:3226811/1‑100 (MQ=255)

           cAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGt                                                                                        >  1:3902870/1‑99 (MQ=255)

               gCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAACATTCAGACTTAGTACTATGATGTGATAAAAGTcaca                                                                                    >  1:196616/1‑99 (MQ=255)

               gCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACAt                                                                                   >  1:282211/1‑100 (MQ=255)

                   cGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACAttttt                                                                               >  1:1992682/1‑100 (MQ=255)

                   cGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACAttttt                                                                               >  2:2507759/1‑100 (MQ=255)

                    gagagaGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGCCTTAGTACTATTATGTGATAGAAGTCACAtttt                                                                                >  1:2663003/1‑98 (MQ=255)

                    gagagaGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTc                                                                              <  2:26256/100‑1 (MQ=255)

                    gagagaGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTc                                                                              >  1:2187198/1‑100 (MQ=255)

                    gagagaGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTc                                                                              <  1:354709/100‑1 (MQ=255)

                    gagagaGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTc                                                                              <  1:4104694/100‑1 (MQ=255)

                    gagagaGTCTCTCTTCCATCTTTGTCGAACCGTATACCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTc                                                                              >  2:2993714/1‑100 (MQ=255)

                     agagagTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTcc                                                                             <  2:2260459/100‑1 (MQ=255)

                      gagagTCTCTCTTCCATCTTTGTCGAACCGTATTCGATTTTTTGATAACCTGGAATATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCca                                                                            >  1:3999349/1‑100 (MQ=255)

                             tctctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGa                                                                     <  2:29086/100‑1 (MQ=255)

                             tctctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGa                                                                     >  2:3227116/1‑100 (MQ=255)

                             tctctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTAAAAAAAAAAAAAAAAAGTCACATTTTTCCACATTg                                                                       <  1:3807271/98‑1 (MQ=255)

                             tctctTCCATCTTTGTCGAACCGTATTCCATCTTGAGTTAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGa                                                                     <  1:2873491/100‑1 (MQ=255)

                              ctctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGAt                                                                    >  1:689019/1‑100 (MQ=255)

                                ctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGa                                                                  >  2:897714/1‑100 (MQ=255)

                                ctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGa                                                                  >  2:1494187/1‑100 (MQ=255)

                                ctTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATTTGATAAAAGTCACATTTTTCCACATTGGAt                                                                    >  1:1869597/1‑98 (MQ=255)

                                 ttCTATCATTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGGGATAAAAGTCACATTTTTCCACATTGGATGAg                                                                 <  1:4009985/100‑1 (MQ=255)

                                 ttCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAg                                                                 <  2:2080479/100‑1 (MQ=255)

                                 ttCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAg                                                                 <  2:796921/100‑1 (MQ=255)

                                   ccATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGcc                                                               <  2:945382/100‑1 (MQ=255)

                                   ccATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGcc                                                               <  1:444330/100‑1 (MQ=255)

                                   ccATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGAACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGcc                                                               <  1:424770/100‑1 (MQ=255)

                                      tCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGt                                                            >  1:2894940/1‑100 (MQ=255)

                                      tCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGt                                                            >  2:1121167/1‑100 (MQ=255)

                                         ttGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCaa                                                         >  2:2746562/1‑100 (MQ=255)

                                         tggTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCaa                                                         <  2:330338/98‑1 (MQ=255)

                                          tGTCGAACCGTATTCCATTTTGAGATAACCCGGGAAATTCAGATTTAGTACTATGATGTGATAAAAATAACATGTTTCCACATTGGATGaaatagtc                                                           >  1:1548568/1‑90 (MQ=255)

                                           gTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGa                                                       >  1:887684/1‑100 (MQ=255)

                                            tCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGAc                                                      >  2:2941646/1‑100 (MQ=255)

                                             cGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACt                                                     >  2:3706214/1‑100 (MQ=255)

                                             cGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACt                                                     >  2:4098847/1‑100 (MQ=255)

                                              gAACCGTATTCCATTTTGAGATAACCTGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTAAAATTTTTCCACATTGGATGAGCTAGTCaa                                                         >  1:2452009/1‑95 (MQ=255)

                                              gAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTc                                                    >  2:4113603/1‑100 (MQ=255)

                                              gAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTc                                                    <  2:2497756/100‑1 (MQ=255)

                                              gAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTc                                                    <  1:1601474/100‑1 (MQ=255)

                                               aaCCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCa                                                   <  2:2810452/100‑1 (MQ=255)

                                               aaCCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCa                                                   >  2:2776947/1‑100 (MQ=255)

                                               aaCCGTATTCCATTTTGAGATAACCCGGAAAATTCACATTTAGTACTATTATGTGATAAAAGTCACATTTTACCACATTGGATGAGCCAGTCAAGACTCa                                                   <  1:2046429/100‑1 (MQ=255)

                                                aCCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCaa                                                  >  2:2182752/1‑100 (MQ=255)

                                                 ccGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCTCATTGGATGAGCCAGTCAAGACTCAAt                                                 <  2:45980/100‑1 (MQ=255)

                                                 ccGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAAt                                                 <  2:605644/100‑1 (MQ=255)

                                                     aTTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtt                                             >  2:272179/1‑100 (MQ=255)

                                                      ttCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAttt                                            >  2:3548866/1‑100 (MQ=255)

                                                      ttCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAttt                                            <  2:1143401/100‑1 (MQ=255)

                                                       tCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtttt                                           >  1:2130481/1‑100 (MQ=255)

                                                       tCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtttt                                           <  2:3045846/100‑1 (MQ=255)

                                                        ccATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAttttt                                          >  2:3541899/1‑100 (MQ=255)

                                                         cATTTTGAGATAACCCGAAAAATTCAGACTGAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCAtttttt                                         <  1:3414628/100‑1 (MQ=255)

                                                            tttGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGaa                                      <  2:3822269/100‑1 (MQ=255)

                                                             ttGAGATAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAAc                                     <  1:2649724/100‑1 (MQ=255)

                                                              tGAGATAACCCTGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACa                                    <  2:2018556/100‑1 (MQ=255)

                                                                 gaTAACCCGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACaaan                                 <  2:745760/100‑2 (MQ=255)

                                                                     aCCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGa                              >  1:3779241/1‑99 (MQ=255)

                                                                     aCCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGAt                             <  2:682265/100‑1 (MQ=255)

                                                                     aCCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGAt                             >  2:174990/1‑100 (MQ=255)

                                                                      cccGGAAAATTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATa                            >  2:2127958/1‑100 (MQ=255)

                                                                       ccGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATaa                           <  2:1016654/100‑1 (MQ=255)

                                                                       ccGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATaa                           <  2:1869440/100‑1 (MQ=255)

                                                                       ccGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATaa                           >  2:1630593/1‑100 (MQ=255)

                                                                         ggAAAATTCAGTGTTAGTACTAGTATGGGGTAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAAt                         <  2:3835539/100‑1 (MQ=255)

                                                                           aaaaTTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTACCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAAttt                       <  1:4179065/100‑1 (MQ=255)

                                                                             aaTTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTg                     >  2:2778800/1‑100 (MQ=255)

                                                                              aTTCAGACTTAGTATTATTATGTGATAAAAGTCACATTTTTCCACAGTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgc                    <  1:340784/100‑1 (MQ=255)

                                                                              aTTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCATGACTCAATCATTTTTTGAACAAATAGATAAATTTTgc                    <  1:2278037/100‑1 (MQ=255)

                                                                              aTTCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCGATCATTTTTTGAACAAATAGATAAAtttttc                    <  2:2879303/100‑3 (MQ=255)

                                                                               ttCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgcg                   >  1:1297645/1‑100 (MQ=255)

                                                                               ttCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgcg                   <  1:826733/100‑1 (MQ=255)

                                                                               ntCAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgcg                   >  1:2700439/2‑100 (MQ=255)

                                                                                tCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgcgc                  >  1:2391085/1‑100 (MQ=255)

                                                                                tCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgcg                   <  1:1690089/99‑1 (MQ=255)

                                                                                 cAGACTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTgcgc                  <  2:2366473/99‑1 (MQ=255)

                                                                                  aGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCaa                >  2:2992990/1‑100 (MQ=255)

                                                                                     cTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACg              >  1:623877/1‑99 (MQ=255)

                                                                                     cTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACgg             <  2:1025928/100‑1 (MQ=255)

                                                                                     cTTAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGCATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACgg             <  1:1669586/100‑1 (MQ=255)

                                                                                      ttAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGc            <  1:383179/100‑1 (MQ=255)

                                                                                      ttAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGc            <  2:3283422/100‑1 (MQ=255)

                                                                                       tAGTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCa           >  2:929138/1‑100 (MQ=255)

                                                                                         gTACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGt         >  1:115019/1‑100 (MQ=255)

                                                                                          tACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTg        <  2:3525644/100‑1 (MQ=255)

                                                                                          tACTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTg        >  1:3829246/1‑100 (MQ=255)

                                                                                           aCTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGa       >  1:4259456/1‑100 (MQ=255)

                                                                                           aCTATGATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTg        <  1:845341/99‑1 (MQ=255)

                                                                                               ttatGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGAATAc   <  1:4222898/100‑1 (MQ=255)

CAGTTTGCCGCCAGGGCTGCGAGAGAGTCTCTCTTCCATCTTTGTCGAACCGTATTCCATTTTGAGATAACCCGGAAAATTCAGACTTAGTACTATTATGTGATAAAAGTCACATTTTTCCACATTGGATGAGCCAGTCAAGACTCAATCATTTTTTGAACAAATAGATAAATTTTGCGCAACGGCAGTGAATACA  >  NC_000913/2033362‑2033557

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: `‑`