BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	836,115	A→G	Y181C (TAC→TGC)	ybiC →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	836,115	0	A	G	87.5%	576.9 / 67.3	231	Y181C (TAC→TGC)	ybiC	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (14/15); new base G (123/79); total (137/94)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

Only 100 of 234 total aligned reads displayed.

GCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGTTAACG  >  NC_000913/836024‑836204

gcgcCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTAc                                                                                          >  1:4206914/1‑93 (MQ=255)

  gcCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTgcgc                                                                                        <  2:5827866/93‑1 (MQ=255)

  gcCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTgcgc                                                                                        >  2:4884978/1‑93 (MQ=255)

   ccGTTCCGCGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGcc                                                                                       <  1:6763942/93‑1 (MQ=255)

   ccGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGcc                                                                                       >  1:786868/1‑93 (MQ=255)

   ccGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGcc                                                                                       >  1:6345271/1‑93 (MQ=255)

     gTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGccac                                                                                     >  1:1425072/1‑93 (MQ=255)

      ttCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGccacc                                                                                    >  1:1516071/1‑93 (MQ=255)

      ttCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGccacc                                                                                    <  1:4359206/93‑1 (MQ=255)

      ttCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGccacc                                                                                    >  1:461721/1‑93 (MQ=255)

      ttCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGccacc                                                                                    >  2:7260706/1‑93 (MQ=255)

      ttCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGccacc                                                                                    >  2:6072991/1‑93 (MQ=255)

       tCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGccacca                                                                                   <  2:5920369/93‑1 (MQ=255)

       gccACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGccacca                                                                                   <  1:3450949/92‑1 (MQ=255)

        ccACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAg                                                                                  >  1:7128509/1‑93 (MQ=255)

         cACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAgc                                                                                 <  1:5534744/93‑1 (MQ=255)

           cGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAgcgc                                                                               <  2:3215689/93‑1 (MQ=255)

              tCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCAt                                                                            <  1:3788604/93‑1 (MQ=255)

               cgcgACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCAtt                                                                           <  1:6878778/93‑1 (MQ=255)

                  gACAGGCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCa                                                                        <  2:535023/93‑1 (MQ=255)

                   aCAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCAt                                                                       <  2:6392498/93‑1 (MQ=255)

                      gCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTg                                                                    >  1:2213531/1‑93 (MQ=255)

                      gCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTg                                                                    <  2:1656028/93‑1 (MQ=255)

                      gCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTg                                                                    >  2:4282603/1‑93 (MQ=255)

                      gCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTg                                                                    <  1:6043579/93‑1 (MQ=255)

                         gCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCa                                                                 <  2:421903/93‑1 (MQ=255)

                         gCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCa                                                                 >  1:3207440/1‑93 (MQ=255)

                           tttGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCaaa                                                               >  1:234219/1‑93 (MQ=255)

                           tttGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCaaa                                                               >  2:6251277/1‑93 (MQ=255)

                            ttGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCaaaa                                                              <  1:1349257/93‑1 (MQ=255)

                            ttGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCaaaa                                                              <  2:2023320/93‑1 (MQ=255)

                            ttGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCaaaa                                                              <  1:3095478/93‑1 (MQ=255)

                              ggCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAAcc                                                            <  2:1226789/93‑1 (MQ=255)

                               gCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAAccc                                                           <  2:2772273/93‑1 (MQ=255)

                                cACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCcg                                                          >  1:5452056/1‑93 (MQ=255)

                                 aCCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCcgc                                                         <  2:4675730/93‑1 (MQ=255)

                                 aCCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCcgc                                                         <  2:365128/93‑1 (MQ=255)

                                 aCCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCcgc                                                         >  2:6529892/1‑93 (MQ=255)

                                  ccAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATATGGCAAAACCcgcg                                                        >  2:412637/1‑93 (MQ=255)

                                    aaTCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTc                                                      <  2:3429222/93‑1 (MQ=255)

                                    aaTCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTc                                                      <  1:2088648/93‑1 (MQ=255)

                                     aTCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCg                                                     >  2:3840297/1‑93 (MQ=255)

                                      tCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGc                                                    >  1:7038304/1‑93 (MQ=255)

                                         gTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTg                                                 <  2:1063445/93‑1 (MQ=255)

                                         gTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTg                                                 <  1:4624355/93‑1 (MQ=255)

                                          ttCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTgg                                                >  2:2162523/1‑93 (MQ=255)

                                           tCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGc                                               <  2:4233792/93‑1 (MQ=255)

                                           tCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGc                                               >  2:5874245/1‑93 (MQ=255)

                                             tgtgtgGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCAt                                             >  1:5225482/1‑93 (MQ=255)

                                             tgtgtgGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCAt                                             >  2:7116619/1‑93 (MQ=255)

                                              gtgtgGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATa                                            >  1:5899369/1‑93 (MQ=255)

                                              gtgtgGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATa                                            >  1:51998/1‑93 (MQ=255)

                                              gtgtgGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATa                                            >  1:866561/1‑93 (MQ=255)

                                               tgtgGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGTTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATaa                                           >  2:890134/1‑93 (MQ=255)

                                               tgtgGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATaa                                           >  2:3515228/1‑93 (MQ=255)

                                                gtgGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATaaa                                          <  1:2474116/93‑1 (MQ=255)

                                                gtgGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATaaa                                          >  2:1427611/1‑93 (MQ=255)

                                                gtgGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATaaa                                          >  1:5861227/1‑93 (MQ=255)

                                                 tgGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAg                                         >  2:4781482/1‑93 (MQ=255)

                                                 tgGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAg                                         <  1:278620/93‑1 (MQ=255)

                                                  ggTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAgg                                        <  2:7380251/93‑1 (MQ=255)

                                                   gTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGc                                       <  1:5091139/93‑1 (MQ=255)

                                                   gTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGc                                       >  1:6694829/1‑93 (MQ=255)

                                                    ttttCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCg                                      <  1:4775393/93‑1 (MQ=255)

                                                        cccTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTccc                                  <  2:6773333/93‑1 (MQ=255)

                                                        cccTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTccc                                  <  1:6228892/93‑1 (MQ=255)

                                                        cccTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGACAAAACCCGCGTCGCCTGGCATAAAGGCGTccc                                  <  1:1910716/93‑1 (MQ=255)

                                                         ccTCGTAAAGATAATTTCGCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTcccc                                 >  1:4002329/1‑93 (MQ=255)

                                                         ccTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTcccc                                 >  1:1053861/1‑93 (MQ=255)

                                                          cTCGTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCg                                >  2:6084242/1‑93 (MQ=255)

                                                          cTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCg                                <  1:1398176/93‑1 (MQ=255)

                                                             gTAAAGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTgc                             <  1:6718169/93‑1 (MQ=255)

                                                               aaaGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTgccg                           >  2:766601/1‑93 (MQ=255)

                                                               aaaGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTgccg                           >  1:720559/1‑93 (MQ=255)

                                                                aaGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTgccgc                          >  1:5790279/1‑93 (MQ=255)

                                                                 aGATAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCAGCGTCGCCTGGCATAAAGGCGTCCCCGTgccgcc                         <  2:6999327/93‑1 (MQ=255)

                                                                   aTAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAg                       <  2:7478515/93‑1 (MQ=255)

                                                                   aTAATTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAg                       >  1:7289033/1‑93 (MQ=255)

                                                                      aTTTCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGtt                    <  1:3256254/93‑1 (MQ=255)

                                                                      aTTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGtt                    >  2:3004628/1‑93 (MQ=255)

                                                                       tttCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTg                   <  1:6784918/93‑1 (MQ=255)

                                                                        ttCCCGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGc                  <  2:942258/93‑1 (MQ=255)

                                                                          cccGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCt                >  1:1022060/1‑93 (MQ=255)

                                                                          cccGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCt                >  2:1553918/1‑93 (MQ=255)

                                                                           ccGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTg               >  1:7365672/1‑93 (MQ=255)

                                                                            cGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGa              <  1:6602202/93‑1 (MQ=255)

                                                                            cGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGa              <  1:916812/93‑1 (MQ=255)

                                                                            cGCTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGa              >  2:1284977/1‑93 (MQ=255)

                                                                             gCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGAt             <  1:4470221/93‑1 (MQ=255)

                                                                              cTGTTGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGAtt            >  2:4053442/1‑93 (MQ=255)

                                                                                 ttGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGAc         >  2:4533241/1‑93 (MQ=255)

                                                                                  tGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACg        >  2:3891572/1‑93 (MQ=255)

                                                                                  tGCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACg        <  2:2330818/93‑1 (MQ=255)

                                                                                   gCTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGt       >  2:4915339/1‑93 (MQ=255)

                                                                                    cTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGtt      >  2:5619787/1‑93 (MQ=255)

                                                                                    cTTGATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGtt      >  2:5708402/1‑93 (MQ=255)

                                                                                    cTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGtt      >  1:2440842/1‑93 (MQ=255)

                                                                                       gATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGTTAAc   >  1:288061/1‑93 (MQ=255)

                                                                                       gATTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGTTAAc   <  1:3011799/93‑1 (MQ=255)

                                                                                        aTTGCGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGTTAACg  >  1:6450818/1‑93 (MQ=255)

GCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTGATTACGCCACCAGCGCCATTGCATTTGGCAAAACCCGCGTCGCCTGGCATAAAGGCGTCCCCGTGCCGCCAGGTTGCCTGATTGACGTTAACG  >  NC_000913/836024‑836204

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