BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,207,804	G→A	A180V (GCG→GTG)	ycfA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,207,804	0	G	A	87.2%	606.3 / 66.3	226	A180V (GCG→GTG)	ycfA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (15/14); new base A (90/107); total (105/121)
Fisher's exact test for biased strand distribution p-value = 5.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-03

Only 100 of 239 total aligned reads displayed.

ATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCAGACTAATG  >  NC_000913/1207714‑1207896

aTCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcacg                                                                                            <  1:1085194/93‑1 (MQ=21)

aTCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACcacg                                                                                            <  1:2960656/93‑1 (MQ=21)

 tCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGt                                                                                           <  1:3865351/93‑1 (MQ=18)

    aGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCaa                                                                                        >  2:5351708/1‑93 (MQ=34)

     gTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAAt                                                                                       <  2:6635318/93‑1 (MQ=33)

     gTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAAt                                                                                       >  2:3527172/1‑93 (MQ=21)

     gTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAAt                                                                                       <  2:252225/93‑1 (MQ=18)

      tGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATg                                                                                      <  2:3266421/93‑1 (MQ=18)

         gcCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTaa                                                                                   >  1:271724/1‑93 (MQ=21)

          ccAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAAt                                                                                  >  2:944009/1‑93 (MQ=21)

          ccAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAAt                                                                                  >  1:2539600/1‑93 (MQ=21)

           cAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATc                                                                                 <  1:2005212/93‑1 (MQ=18)

           cAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATc                                                                                 <  1:4363744/93‑1 (MQ=18)

            aGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATcc                                                                                <  1:5243532/93‑1 (MQ=18)

             gATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCa                                                                               >  2:4517664/1‑93 (MQ=34)

             gATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCa                                                                               >  1:2648892/1‑93 (MQ=21)

              atatTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAg                                                                              >  1:3729024/1‑93 (MQ=34)

               tatTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGc                                                                             <  2:1494983/93‑1 (MQ=18)

               tatTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGc                                                                             >  2:5696648/1‑93 (MQ=21)

                 ttGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGcac                                                                           >  2:3376940/1‑93 (MQ=21)

                  tGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGcaca                                                                          <  2:6258797/93‑1 (MQ=18)

                  tGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGcaca                                                                          >  2:5188870/1‑93 (MQ=21)

                   gAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAg                                                                         <  2:5513572/93‑1 (MQ=18)

                     aTGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCg                                                                       >  2:55594/1‑93 (MQ=37)

                      tGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGt                                                                      <  1:3352375/93‑1 (MQ=255)

                      tGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGt                                                                      <  2:6225502/93‑1 (MQ=37)

                       ggCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGtt                                                                     <  2:1146243/93‑1 (MQ=255)

                       ggCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGtt                                                                     <  1:2743581/93‑1 (MQ=255)

                        gCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTa                                                                    <  1:5737261/93‑1 (MQ=37)

                        gCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTa                                                                    >  1:2279055/1‑93 (MQ=255)

                          cTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAg                                                                  <  1:596205/93‑1 (MQ=37)

                           tACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGt                                                                 >  2:673484/1‑93 (MQ=37)

                           tACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGt                                                                 >  1:3621887/1‑93 (MQ=37)

                            aCGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTc                                                                <  1:1312137/93‑1 (MQ=37)

                            aCGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTc                                                                <  2:6573343/93‑1 (MQ=255)

                            aCGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTc                                                                >  2:6305348/1‑93 (MQ=37)

                            aCGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTc                                                                <  1:605575/93‑1 (MQ=255)

                             cGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCt                                                               >  2:2669569/1‑93 (MQ=37)

                             cGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCt                                                               <  1:3971601/93‑1 (MQ=37)

                              gCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTg                                                              >  2:5118414/1‑93 (MQ=37)

                                cTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGgt                                                            <  2:246856/93‑1 (MQ=255)

                                 tCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGgtt                                                           >  2:6488962/1‑93 (MQ=255)

                                 tCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGgtt                                                           >  2:2990031/1‑93 (MQ=255)

                                  ccGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGgttg                                                          <  2:2309352/93‑1 (MQ=255)

                                    ggCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGgttgtt                                                        <  2:4087990/93‑1 (MQ=255)

                                     gCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGttt                                                       >  2:2479516/1‑93 (MQ=255)

                                     gCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGttt                                                       >  2:4387017/1‑93 (MQ=255)

                                      cAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTc                                                      <  2:4883638/93‑1 (MQ=255)

                                      cAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTc                                                      <  2:527554/93‑1 (MQ=255)

                                      cAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTc                                                      <  1:4048567/93‑1 (MQ=255)

                                       aGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTct                                                     >  2:1665965/1‑93 (MQ=255)

                                        gTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctg                                                    >  1:4572220/1‑93 (MQ=255)

                                         ttCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgc                                                   >  2:6019018/1‑93 (MQ=255)

                                         ttCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgc                                                   >  2:2098729/1‑93 (MQ=255)

                                          tCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                                                  <  2:3589463/93‑1 (MQ=255)

                                          tCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcc                                                  <  1:4815510/93‑1 (MQ=255)

                                           cAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcct                                                 <  1:791899/93‑1 (MQ=255)

                                             ggCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTctgcctgc                                               <  2:6332279/93‑1 (MQ=255)

                                              gCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCg                                              >  1:1410085/1‑93 (MQ=255)

                                                cAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTc                                            >  1:1164644/1‑93 (MQ=255)

                                                 aGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCa                                           <  1:1259773/93‑1 (MQ=255)

                                                   aTGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGt                                         >  2:796597/1‑93 (MQ=255)

                                                   aTGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGt                                         >  1:1689418/1‑93 (MQ=255)

                                                   aTGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGt                                         >  2:3072005/1‑93 (MQ=255)

                                                     gACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGttt                                       >  2:1353426/1‑93 (MQ=255)

                                                     gACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGttt                                       <  2:1077220/93‑1 (MQ=255)

                                                      aCATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTa                                      <  2:849203/93‑1 (MQ=255)

                                                       cATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTAc                                     >  2:3150849/1‑93 (MQ=255)

                                                        aTCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACg                                    >  2:5946592/1‑93 (MQ=255)

                                                         tCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGt                                   >  2:5832573/1‑93 (MQ=255)

                                                          ccGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTc                                  <  2:618839/93‑1 (MQ=255)

                                                           cGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTcc                                 <  1:3261783/93‑1 (MQ=255)

                                                           cGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTcc                                 <  2:3399168/93‑1 (MQ=255)

                                                           cGGCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTcc                                 <  1:6224960/93‑1 (MQ=255)

                                                            ggCGCGGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCg                                >  1:2753684/1‑93 (MQ=255)

                                                                cgGTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCt                            >  2:3326936/1‑93 (MQ=255)

                                                                  gTGCTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGc                          >  1:6570300/1‑93 (MQ=255)

                                                                     cTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAAt                       <  1:4421629/93‑1 (MQ=255)

                                                                     cTGGTATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAAt                       <  2:4961242/93‑1 (MQ=255)

                                                                      tGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAAtt                      >  1:3524442/1‑93 (MQ=255)

                                                                        gTATCTGTTGCCGTCACCACGTCAATGTAGTCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTc                    >  1:4525191/1‑93 (MQ=255)

                                                                         tATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCa                   <  2:3542270/93‑1 (MQ=255)

                                                                         tATCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCa                   <  1:3047734/93‑1 (MQ=255)

                                                                          aTCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAg                  <  1:725358/93‑1 (MQ=255)

                                                                          aTCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAg                  <  1:2203905/93‑1 (MQ=255)

                                                                           tCTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGt                 >  1:1805223/1‑93 (MQ=255)

                                                                            cTGTTGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGtt                <  1:4099025/93‑1 (MQ=255)

                                                                              gTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGa              >  1:477369/1‑93 (MQ=255)

                                                                               ttGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGaa             >  2:1375021/1‑93 (MQ=255)

                                                                                tGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAAt            >  2:3729793/1‑93 (MQ=255)

                                                                                tGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAAt            <  1:544851/93‑1 (MQ=255)

                                                                                tGCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAAt            >  1:4568817/1‑93 (MQ=255)

                                                                                 gCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATc           <  1:3027338/93‑1 (MQ=255)

                                                                                 gCCGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATc           <  2:2596089/93‑1 (MQ=255)

                                                                                   cGTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCAg         >  1:5339334/1‑93 (MQ=255)

                                                                                    gTCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCTTGCAATTTCAGTTGAATCAGa        <  2:5475863/93‑1 (MQ=255)

                                                                                     tCACCACGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCAGAc       >  2:4262928/1‑93 (MQ=255)

                                                                                      caccacGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCAGACt      <  1:4170738/93‑1 (MQ=255)

                                                                                      caccacGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCAGACt      <  1:62338/93‑1 (MQ=255)

                                                                                      caccacGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCAGACt      <  1:1075049/93‑1 (MQ=255)

ATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGTAATCCAGCACAGCGTTAAGTCTGGTTGTTTCTGCCTGCGTCAGTTTACGTCCGGCCTGCAATTTCAGTTGAATCAGACTAATG  >  NC_000913/1207714‑1207896

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