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	97.5%	408.9 / ‑1.8	121	M19V (ATG→GTG)	yciA	acyl‑CoA esterase
Reads supporting (aligned to +/- strand): ref base T (2/1); new base C (55/63); total (57/64)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 121 total aligned reads displayed.

GGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATG  >  NC_000913/1312096‑1312284

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

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

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

 ggCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCAACATTGGCATTGGTATCGGCGGGCACGGc                                                                                           >  2:1522602/1‑99 (MQ=255)

  gCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTa                                                                                         >  2:1310385/1‑100 (MQ=255)

  gCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTa                                                                                         <  1:442479/100‑1 (MQ=255)

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

       ttCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTa                                                                                    <  2:1282639/100‑1 (MQ=255)

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

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

         cTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACg                                                                                  <  2:2525383/100‑1 (MQ=255)

         cTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACg                                                                                  >  1:302816/1‑100 (MQ=255)

          ttttGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGt                                                                                 >  2:1072387/1‑100 (MQ=255)

             tGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTaaa                                                                              <  2:2608082/100‑1 (MQ=255)

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

              gCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTaaaa                                                                             <  1:2494749/100‑1 (MQ=255)

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

               ccAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCACAGATGTAACCGTTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTaaacc                                                                            <  2:1831462/100‑3 (MQ=255)

               ccAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTAACCATTGGAATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAAc                                                                            <  2:1988429/100‑1 (MQ=255)

                 aGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACaa                                                                          >  2:2253976/1‑100 (MQ=255)

                  gAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAg                                                                         <  2:180675/100‑1 (MQ=255)

                   aaTAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAg                                                                         <  2:368648/99‑1 (MQ=255)

                          ccgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCt                                                                 >  1:184557/1‑100 (MQ=255)

                           cgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTg                                                                <  2:627790/100‑1 (MQ=255)

                           cgccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTg                                                                <  1:2517340/100‑1 (MQ=255)

                             ccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAg                                                              <  1:2989302/100‑1 (MQ=255)

                             ccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAg                                                              <  1:1092275/100‑1 (MQ=255)

                             ccAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAg                                                              <  1:2871070/100‑1 (MQ=255)

                                 tatCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCCGAGGGAc                                                          >  2:640433/1‑100 (MQ=255)

                                 tatCCATTTGTGACATGAACCAACCACCAAAGATGTCACGATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAGACAAGAGCGCCCTTAGGGAc                                                          >  1:2676368/1‑100 (MQ=255)

                                  atCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACg                                                         <  2:87105/100‑1 (MQ=255)

                                   tCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGt                                                        >  1:626892/1‑100 (MQ=255)

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

                                     cATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTa                                                      <  2:2547969/100‑1 (MQ=255)

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

                                      aTTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTAt                                                     >  1:2124340/1‑100 (MQ=255)

                                       tttGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTAtg                                                    <  2:256960/100‑1 (MQ=255)

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

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

                                          gtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATgtgt                                                 <  2:707604/100‑1 (MQ=255)

                                          gtgACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATgtgt                                                 >  1:1493034/1‑100 (MQ=255)

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

                                            gACATTAACCATCCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgttg                                               <  1:30944/100‑1 (MQ=255)

                                            gACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGtgttg                                               >  2:1676118/1‑100 (MQ=255)

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

                                              cATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTa                                             <  2:2528329/100‑1 (MQ=255)

                                              cATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTa                                             <  2:3324140/100‑1 (MQ=255)

                                                ttAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGa                                           <  2:2466988/100‑1 (MQ=255)

                                                  aaccaaccACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACa                                         >  1:204814/1‑100 (MQ=255)

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

                                                   accaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACAt                                        <  2:932478/100‑1 (MQ=255)

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

                                                    ccaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAATATCGCCCTGAGGGACGTTATGTGTTGTAGACATg                                       >  2:1774989/1‑100 (MQ=255)

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

                                                    ccaaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATg                                       >  2:3339880/1‑100 (MQ=255)

                                                     caaccaCCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATgg                                      <  1:2885102/100‑1 (MQ=255)

                                                     caaccaCCAAAGATGTCACCATTGGCATTGGTAACGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCGCTGAGGGACGTTATGTGTTGTAGACATgg                                      >  2:2018346/1‑100 (MQ=255)

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

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

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

                                                       accaccAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTa                                    <  1:2203281/100‑1 (MQ=255)

                                                       accaccAAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTa                                    >  2:2116023/1‑100 (MQ=255)

                                                        ccaccaAATATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAGGATCGCCCTGAGGGACGTTAGGTGTTGTAGACATGGt                                     <  2:417959/98‑1 (MQ=255)

                                                        ccaccaAAGATGTCACCATTGGCATTTGTATCGGCGGGCACGGCGAAAGTACGTAAATCAATATCGCCCTGACGGACGTTATGTGTTGTAGACATGGTa                                    <  2:3281730/99‑1 (MQ=255)

                                                        ccaccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGt                                     <  1:3163081/98‑1 (MQ=255)

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

                                                         caccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTaaa                                  >  1:1627636/1‑100 (MQ=255)

                                                         caccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTaaa                                  >  1:1592667/1‑100 (MQ=255)

                                                         caccaAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTaaa                                  >  2:1680312/1‑100 (MQ=255)

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

                                                          accaAAGATGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTaaaa                                 >  1:3211658/1‑100 (MQ=255)

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

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

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

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

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

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

                                                                  tGTCACCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTa                         >  1:2442965/1‑100 (MQ=255)

                                                                      aCCATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAg                     <  1:1100259/100‑1 (MQ=255)

                                                                       ccATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAgg                    <  1:2541541/100‑1 (MQ=255)

                                                                       ccATTGGCATTGGTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAgg                    >  2:574067/1‑100 (MQ=255)

                                                                            ggcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAAt               >  1:943403/1‑100 (MQ=255)

                                                                            ggcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAAt               >  1:1494635/1‑100 (MQ=255)

                                                                            ggcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAAt               <  2:1930655/100‑1 (MQ=255)

                                                                             gcattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATc              <  1:1178986/100‑1 (MQ=255)

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

                                                                              cattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATca             >  2:2902207/1‑100 (MQ=255)

                                                                              cattggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATca             >  2:876547/1‑100 (MQ=255)

                                                                               attggTATCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATcac            <  1:2462344/100‑1 (MQ=255)

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

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

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

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

                                                                                   gTATCGGCGGGCACGGCGAAAGCCCGTGAAACCATGCCGCCCTGAGGGACGTTGTGTGTTGTCGCCCCGGTCATAGCGACATAAAATGAAAATtac            >  2:368141/1‑96 (MQ=255)

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

                                                                                      tCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGcgcg     <  1:1676122/100‑1 (MQ=255)

                                                                                      tCGGCGGGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGcgcg     >  2:2390389/1‑100 (MQ=255)

                                                                                        ggcggGCACGGCTAAAGGACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATcatagtagcatc   <  2:2670988/100‑11 (MQ=255)

                                                                                         gcggGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATg  >  1:1983216/1‑100 (MQ=255)

                                                                                          cggGCACGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATg  <  2:3152180/99‑1 (MQ=255)

GGGCAATTTCTTTTGCCAGAATAGCGCCGCCAATATCCATTTGTGACATTAACCAACCACCAAAGATGTCACCATTGGCATTGGTATCGGCGGGCATGGCTAAAGTACGTAAAACAAGATCGCCCTGAGGGACGTTATGTGTTGTAGACATGGTAAAACCGACTTAAAAGGAAAATCACAGGCGCGATG  >  NC_000913/1312096‑1312284

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