BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,312,192	T→C	M19V (ATG→GTG)	yciA ←	acyl‑CoA esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,312,192	0	T	C	98.3%	405.7 / ‑4.0	117	M19V (ATG→GTG)	yciA	acyl‑CoA esterase
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (64/51); total (65/52)
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.89e-01

Only 100 of 117 total aligned reads displayed.

TGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGAT  >  NC_000913/1312095‑1312283

tGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACgg                                                                                           <  2:1274736/100‑1 (MQ=255)

 gggCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGc                                                                                          <  2:333265/100‑1 (MQ=255)

 gggCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGc                                                                                          <  2:1527071/100‑1 (MQ=255)

 gggCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGc                                                                                          >  1:1767462/1‑100 (MQ=255)

    cAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTaa                                                                                       <  2:509052/100‑1 (MQ=255)

     aaTTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTaaa                                                                                      >  1:2427009/1‑100 (MQ=255)

      aTTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAg                                                                                     >  1:273223/1‑100 (MQ=255)

        ttCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTa                                                                                   >  1:34599/1‑100 (MQ=255)

        ttCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTa                                                                                   >  2:2298675/1‑100 (MQ=255)

         tCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTAc                                                                                  >  1:263069/1‑100 (MQ=255)

         tCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTAc                                                                                  >  1:2141745/1‑100 (MQ=255)

         tCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTAc                                                                                  >  2:1076516/1‑100 (MQ=255)

           ttttGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGt                                                                                <  2:2421721/100‑1 (MQ=255)

            tttGCCAGAATAGCGCCGCCAATATCCATTTGTGACGTTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTa                                                                               >  2:2046227/1‑100 (MQ=255)

              tGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTaaa                                                                             >  1:2625067/1‑100 (MQ=255)

               gCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTaaaa                                                                            >  1:2570784/1‑100 (MQ=255)

                ccAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAAc                                                                           <  1:2692414/100‑1 (MQ=255)

                ccAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAAc                                                                           >  1:2506131/1‑100 (MQ=255)

                 cAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACa                                                                          >  1:1882562/1‑100 (MQ=255)

                     aTAGCGCCGCCAGTATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATTGGCGGGCACGGCTAAAGTACGTAAAACAAGAt                                                                      <  1:2671550/100‑1 (MQ=255)

                      tAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATc                                                                     <  2:752848/100‑1 (MQ=255)

                      tAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATc                                                                     <  1:1025574/100‑1 (MQ=255)

                           ccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATcacac                                                                 >  2:554116/1‑95 (MQ=255)

                            cgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTg                                                               <  2:180098/100‑1 (MQ=255)

                            cgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCt                                                                >  1:2031294/1‑99 (MQ=255)

                             gccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGa                                                              >  2:2062910/1‑100 (MQ=255)

                             gccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGa                                                              <  1:394238/100‑1 (MQ=255)

                             gccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGa                                                              >  2:1458396/1‑100 (MQ=255)

                               cAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAgg                                                            <  2:686191/100‑1 (MQ=255)

                                    tCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGt                                                       <  1:1021770/100‑1 (MQ=255)

                                     ccATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGtt                                                      >  1:830885/1‑100 (MQ=255)

                                     ccATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGtt                                                      <  2:2334258/100‑1 (MQ=255)

                                     ccATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGtt                                                      <  2:1067258/100‑1 (MQ=255)

                                      cATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGCGGGACGTTa                                                     >  1:2473770/1‑100 (MQ=255)

                                      cATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTa                                                     >  1:990268/1‑100 (MQ=255)

                                          tgtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTAtgtg                                                 <  2:194655/100‑1 (MQ=255)

                                          tgtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAATACGTAAAACAAGATCGCCCTGACGGAAGTTAtgtg                                                 >  1:2302875/1‑100 (MQ=255)

                                           gtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATgtgt                                                <  1:429540/100‑1 (MQ=255)

                                            tgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgtt                                               >  1:2081463/1‑100 (MQ=255)

                                              aCATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgttgt                                             <  2:822243/100‑1 (MQ=255)

                                              aCATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgttgt                                             >  1:2556092/1‑100 (MQ=255)

                                              aCATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgatgt                                             >  1:170457/1‑100 (MQ=255)

                                               cATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTa                                            >  2:2444151/1‑100 (MQ=255)

                                               cATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTa                                            >  1:1708930/1‑100 (MQ=255)

                                               cATTAACCAACCACCAAAGATGTCACCAATGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTa                                            <  1:887244/100‑1 (MQ=255)

                                                aTTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAg                                           >  2:820948/1‑100 (MQ=255)

                                                aTTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAg                                           <  1:2148388/100‑1 (MQ=255)

                                                 ttgaccaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGa                                          <  1:869826/97‑1 (MQ=255)

                                                  tAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGAc                                         >  2:1968200/1‑100 (MQ=255)

                                                   aaccaaccACCAAAGATGTCTCCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACCAGATCGCCCTGAGGGACCTTATGTGTTGTAg                                           <  2:1954592/97‑1 (MQ=255)

                                                    accaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACAt                                       >  2:793859/1‑100 (MQ=255)

                                                    accaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACAt                                       >  1:2292962/1‑100 (MQ=255)

                                                    accaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACAt                                       >  1:2293293/1‑100 (MQ=255)

                                                     ccaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATg                                      <  2:2556125/100‑1 (MQ=255)

                                                     ccaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATg                                      <  2:887334/100‑1 (MQ=255)

                                                     ccaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGACTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATg                                      >  2:208889/1‑100 (MQ=255)

                                                      caaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGACTAAAGTACGTAAAACAAGAACGCCCTGAGGGACGTTATGTGTTGTAGACATgg                                     >  1:997927/1‑100 (MQ=255)

                                                       aaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGt                                    >  1:2242293/1‑100 (MQ=255)

                                                       aaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGt                                    <  1:2440152/100‑1 (MQ=255)

                                                         ccaccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTaa                                  <  2:1353888/100‑1 (MQ=255)

                                                           accaAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTaaaa                                <  2:1646309/100‑1 (MQ=255)

                                                            ccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAAc                               >  1:1722278/1‑100 (MQ=255)

                                                            ccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAAc                               <  1:751909/100‑1 (MQ=255)

                                                             caAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAAcc                              >  2:1660360/1‑100 (MQ=255)

                                                             caAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAAcc                              >  2:166522/1‑100 (MQ=255)

                                                              aaaGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCg                             >  2:1792221/1‑100 (MQ=255)

                                                              aaaGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCg                             >  2:2505007/1‑100 (MQ=255)

                                                               aaGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGa                            <  2:1881980/100‑1 (MQ=255)

                                                               aaGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGa                            <  2:1882588/100‑1 (MQ=255)

                                                                aGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGAc                           <  2:28563/100‑1 (MQ=255)

                                                                 gATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACt                          <  1:2457007/100‑1 (MQ=255)

                                                                  aTGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACtt                         <  2:2847869/100‑1 (MQ=255)

                                                                   tGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTa                        <  1:1947453/100‑1 (MQ=255)

                                                                   tGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTa                        <  2:648574/100‑1 (MQ=255)

                                                                    gTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTaa                       >  2:1656996/1‑100 (MQ=255)

                                                                    gTCACCATNGGCATTGGTATCGGCGGNCACGNCTNAANTACGNANAACAAGGTCGCCCTGCGGGACGTTATGTGTTGTAGACATGGTaaagtcnacgta                        >  2:2257559/1‑90 (MQ=255)

                                                                      cACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTaaaa                     <  1:2073968/100‑1 (MQ=255)

                                                                           ttggcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGaaa                >  1:1209515/1‑100 (MQ=255)

                                                                           ttggcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGaaa                <  1:1357560/100‑1 (MQ=255)

                                                                           ttggcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGaaa                >  2:724787/1‑100 (MQ=255)

                                                                              gcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAg                    <  1:2504975/93‑1 (MQ=255)

                                                                              gcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATc             >  2:1665372/1‑100 (MQ=255)

                                                                              gcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATc             >  2:4039/1‑100 (MQ=255)

                                                                              gcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTAGGTGTTGTAGACATGGTAATACCGACTTAAAAGGAAAATc             <  1:1769481/100‑1 (MQ=255)

                                                                                 ttggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATcaca          <  2:34600/100‑1 (MQ=255)

                                                                                  tggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAg         <  2:2625100/100‑1 (MQ=255)

                                                                                   ggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAgg        <  2:1319659/100‑1 (MQ=255)

                                                                                   ggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAgg        >  1:1090995/1‑100 (MQ=255)

                                                                                   ggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAgg        >  2:697561/1‑100 (MQ=255)

                                                                                   ggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGAAAAGGTAAAACCGACTTAAAAGGAAAATCACAgg        >  1:1001225/1‑100 (MQ=255)

                                                                                    gTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGgc       <  1:924501/100‑1 (MQ=255)

                                                                                    gTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGgc       >  1:584173/1‑100 (MQ=255)

                                                                                      aTCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGgcgc     <  2:394742/100‑1 (MQ=255)

                                                                                      aTCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGgcgc     >  1:1955651/1‑100 (MQ=255)

                                                                                        cggcggGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGa   >  1:96832/1‑100 (MQ=255)

                                                                                        cggcggGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGa   >  1:1708543/1‑100 (MQ=255)

                                                                                        cggcggGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGa   >  1:2822447/1‑100 (MQ=255)

                                                                                         ggcggGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGAt  <  2:2724425/100‑1 (MQ=255)

                                                                                         ggcggGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGAt  >  2:1338468/1‑100 (MQ=255)

                                                                                           cggGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGTGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGa   >  1:435986/1‑97 (MQ=255)

TGGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGAT  >  NC_000913/1312095‑1312283

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