BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	11,127	C→G	15.8%	S77T (AGT→ACT)	yaaW ←	UPF0174 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	11,127	0	C	G	15.8%	333.1 / 52.9	126	S77T (AGT→ACT)	yaaW	UPF0174 family protein
Reads supporting (aligned to +/- strand): ref base C (49/56); new base G (8/12); total (58/68)
Fisher's exact test for biased strand distribution p-value = 6.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-01

Only 100 of 131 total aligned reads displayed.

GTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCA  >  NC_000913/11028‑11226

gTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATAc                                                                                                     <  1:1784279/100‑1 (MQ=255)

 tAGACATCTCTTTGTCGTCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTCGGTGGCGATAc                                                                                                     <  1:825136/99‑1 (MQ=255)

  aGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTa                                                                                                   <  1:1259791/100‑1 (MQ=255)

    aCATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATc                                                                                                 <  2:662380/100‑1 (MQ=255)

       tctcTTTGTCGGCTTTCAGCTTCACTCGCTTTGAAACATCGAGCAAAATGGCCAGATACAATTTACGGTGTGCGCGCAGTTTGTTGGCGATACTATCGcc                                                                                              >  2:1095391/1‑100 (MQ=255)

       tctcTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGcc                                                                                              <  2:345353/100‑1 (MQ=255)

       tctcTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCAcc                                                                                              <  1:1244023/100‑1 (MQ=255)

        ttctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGcca                                                                                             <  2:619697/99‑1 (MQ=255)

        ctctTTGTCGGCTTTCAGCTTCAATCGCTTGGAGACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGAGACTATCGcca                                                                                             >  1:1699703/1‑100 (MQ=255)

        ccgatcttCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGcca                                                                                             <  1:817951/93‑1 (MQ=255)

         tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATATTATCGccac                                                                                            >  1:145342/1‑100 (MQ=255)

         tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGccac                                                                                            >  2:1132253/1‑100 (MQ=255)

         tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGccac                                                                                            <  1:2177919/100‑1 (MQ=255)

         tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGccac                                                                                            <  2:1405996/100‑1 (MQ=255)

         tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGccac                                                                                            >  2:2108820/1‑100 (MQ=255)

         tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGccac                                                                                            >  2:1382609/1‑100 (MQ=255)

         tctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGccaa                                                                                            >  1:619694/1‑99 (MQ=255)

          ctTTGTCGGCTTTCATCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGccacc                                                                                           <  2:573278/100‑1 (MQ=255)

          ctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATAGTATCGccaac                                                                                           <  2:1018915/100‑3 (MQ=255)

          ctTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATAGTATCGccaac                                                                                           <  2:555478/100‑3 (MQ=255)

           tttGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATATAATTTACCGTGTCCGCGCATTTTGGTGGTGATAGTATCGccaaca                                                                                          <  1:11618/100‑1 (MQ=38)

           tttGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGccacca                                                                                          <  1:972016/100‑1 (MQ=255)

             tGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCCAAATGGCCCGATACAATTTACCGTATCCGCGCAGTTTGTTGACGATACTATCGccaccta                                                                                         >  1:84990/1‑97 (MQ=255)

              gTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTCACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCaa                                                                                         >  1:1653338/1‑98 (MQ=255)

               tCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGccaagatcgg                                                                                      >  2:817954/1‑93 (MQ=255)

               tCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAAt                                                                                      <  2:406165/100‑1 (MQ=255)

               tCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAAt                                                                                      >  2:223737/1‑100 (MQ=255)

                ctgCTTACAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATAGTATCGCCAACaaac                                                                                       <  1:2187899/96‑2 (MQ=38)

                cGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGGTGGCGATAGTATCGCCAACAAACTg                                                                                     <  2:1830902/100‑1 (MQ=38)

                cGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGGTGGCGATAGTATCGCCAACAAACTg                                                                                     <  1:610274/100‑1 (MQ=38)

                 ggCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGc                                                                                    <  2:973671/100‑1 (MQ=255)

                 ggCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGc                                                                                    <  2:2289164/100‑1 (MQ=255)

                 ggCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGc                                                                                    <  2:344984/100‑1 (MQ=255)

                 ggCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGGTGGCGATAGTATCGCCAACAAACTGc                                                                                    <  1:1200527/100‑1 (MQ=38)

                  gCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCATTTTGGTGGCGATAGTATCGCCAACAAACTGcg                                                                                   <  2:2129233/100‑2 (MQ=37)

                  gCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCt                                                                                   >  1:2432504/1‑100 (MQ=255)

                   cTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCATTTTGGTGGCGATAGTATCGCCAACAAACTGcgc                                                                                  <  2:586897/100‑3 (MQ=25)

                    tttCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGt                                                                                 <  1:2210628/100‑1 (MQ=255)

                      tCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACTGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAActgctg                                                                                  >  1:1609372/1‑97 (MQ=255)

                      tCAGCTTCAATCACTTTGAAACATCGAGCAAAATGGTCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGt                                                                                 <  1:1279029/98‑1 (MQ=255)

                       cAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAAt                                                                              >  2:11676/1‑100 (MQ=255)

                          cTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATtct                                                                           >  2:2099899/1‑100 (MQ=255)

                             cAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCg                                                                        >  2:1048025/1‑100 (MQ=255)

                              ttccgatctGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCgg                                                                       <  2:2371250/92‑1 (MQ=255)

                              aaTCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCgg                                                                       <  1:628632/100‑1 (MQ=255)

                                  gCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAAt                                                                   >  2:1076993/1‑100 (MQ=255)

                                   cTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATc                                                                  >  2:2323564/1‑100 (MQ=255)

                                       gAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGACACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCt                                                              <  1:1412307/100‑1 (MQ=255)

                                         aaCATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGc                                                            >  2:2374622/1‑100 (MQ=255)

                                         aaCATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGc                                                            <  1:1074573/100‑1 (MQ=255)

                                          aCATCGAGCAAAAGGTCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGGAATGCTCCGGCAATCAGCTGcc                                                           <  1:247204/100‑1 (MQ=255)

                                             tCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTTTTGGCGATACTGTCGCCACCAAAAGGCTGTAATTCTCCGGCAATCATCCGCCAg                                                         <  1:1541110/99‑1 (MQ=255)

                                             tCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGt                                                        <  1:353830/100‑1 (MQ=255)

                                               gAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTg                                                      >  2:2461858/1‑100 (MQ=255)

                                                   aaaaTGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTgcggc                                                  <  2:783776/100‑1 (MQ=255)

                                                   aaaaTGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTgcggc                                                  >  2:1636659/1‑100 (MQ=255)

                                                    aaaTGGCCCGATACAATTTACTGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCACCTGCCAGTTgcggcg                                                 >  1:1339589/1‑100 (MQ=255)

                                                     aaTGGCCCCATACAATTTACCGTGCCCGCGCATTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGc                                                               <  1:2544012/85‑1 (MQ=255)

                                                      aTGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGAt                                               <  1:22982/100‑1 (MQ=255)

                                                       tGGCGCGATACCATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATg                                              <  1:975538/100‑1 (MQ=255)

                                                        ggCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGt                                             >  1:314422/1‑100 (MQ=255)

                                                         gCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGtt                                            >  1:2400569/1‑100 (MQ=255)

                                                                 cAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGa                                    >  2:80520/1‑100 (MQ=255)

                                                                 cAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGa                                    <  1:674270/100‑1 (MQ=255)

                                                                   aTTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATg                                  >  1:1711290/1‑100 (MQ=255)

                                                                    tttCCCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATg                                  <  1:1675141/99‑1 (MQ=255)

                                                                    tttACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGc                                 >  1:1767240/1‑100 (MQ=255)

                                                                    tttACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGc                                 <  1:300630/100‑1 (MQ=255)

                                                                    tttACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGc                                 >  2:1922165/1‑100 (MQ=255)

                                                                     ttACCGTGTCCGCGCAGTTTGTTGGCGATACTGTCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGcc                                <  1:2418665/100‑1 (MQ=255)

                                                                     ttACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGcc                                <  1:2076794/100‑1 (MQ=255)

                                                                     ttACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGcc                                <  2:473840/100‑1 (MQ=255)

                                                                      tACCGTGTCCGCGCAGTTTGTTGGCGATACGATCGCCACCAAAATGCTGTAATGCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGc                                 <  1:1775782/98‑1 (MQ=255)

                                                                       aCCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCt                              <  1:1326185/100‑1 (MQ=255)

                                                                       aCCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCt                              <  1:2094878/100‑1 (MQ=255)

                                                                         cGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTc                            >  1:750260/1‑100 (MQ=255)

                                                                         cGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTc                            >  1:1989585/1‑100 (MQ=255)

                                                                            gtCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCAt                         >  2:1413779/1‑100 (MQ=255)

                                                                             ttacaacattttGGTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATc                        <  2:15308/91‑1 (MQ=255)

                                                                                gcgcAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGAt                     <  2:1733848/100‑1 (MQ=255)

                                                                                gcgcAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGAt                     <  2:742220/100‑1 (MQ=255)

                                                                                cagcattttGGTGGCGATAGTATCGCCAACAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGAt                     >  2:1462200/7‑100 (MQ=37)

                                                                                  gcAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGAttt                   <  2:1909302/100‑1 (MQ=255)

                                                                                    attttGGTGGCGATAGTATCGCCAACAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTaa                 >  2:2504119/3‑100 (MQ=38)

                                                                                    aGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTaa                 <  2:2400584/100‑1 (MQ=255)

                                                                                    aGTTTGTTGGCGATACTATCGCCACCAAAATGCTGAAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTaa                 <  1:803265/100‑1 (MQ=255)

                                                                                        tgttgGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAg             >  1:739283/1‑100 (MQ=255)

                                                                                        tgttgGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAg             >  2:1493856/1‑100 (MQ=255)

                                                                                        tgttgGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAg             <  1:1090375/100‑1 (MQ=255)

                                                                                         gttgGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGt            >  2:1761989/1‑100 (MQ=255)

                                                                                         ggtgGCGATAGTATCGCCAACAAACTGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGt            >  2:2312072/3‑100 (MQ=38)

                                                                                         ggtgGCGATAGTATCGCCAACAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGGTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGt            >  2:118433/3‑100 (MQ=255)

                                                                                          ttgGCGATACTATCGCCACCAAAATGCTGTAATTCGCCGGCAATCAGCTGCCAGTTGCTGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGtt           <  1:264970/100‑1 (MQ=255)

                                                                                           tgGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTc          >  2:346396/1‑100 (MQ=255)

                                                                                            ggCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCg         >  1:1077662/1‑100 (MQ=255)

                                                                                             gCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGt        >  2:1143662/1‑100 (MQ=255)

                                                                                                 tACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTgcg    <  1:1493847/100‑1 (MQ=255)

                                                                                                   cTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCa  <  2:2561217/100‑1 (MQ=255)

                                                                                                   cTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCa  <  2:548537/100‑1 (MQ=255)

                                                                                                   cTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCa  <  1:1219915/100‑1 (MQ=255)

GTAGACATCTCTTTGTCGGCTTTCAGCTTCAATCGCTTTGAAACATCGAGCAAAATGGCCCGATACAATTTACCGTGTCCGCGCAGTTTGTTGGCGATACTATCGCCACCAAAATGCTGTAATTCTCCGGCAATCAGCTGCCAGTTGCGGCGATGTTGCTCGGGATGCCCTTCCATCGATTTAAACAGTTCGTTGCGCA  >  NC_000913/11028‑11226

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