BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,008,480	+A	coding (637/654 nt)	yedK →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,008,476	1	.	A	99.2%	914.4 / ‑1.3	248	coding (633/654 nt)	yedK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/2); new base A (129/117); total (129/119)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 277 total aligned reads displayed.

GCGGTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑AAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTATTA  >  NC_000913/2008385‑2008563

gCGGTAAAGAAGCCCCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑a                                                                                         <  2:8543640/93‑1 (MQ=255)

  ggTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑aaa                                                                                       <  1:9736226/93‑1 (MQ=255)

  ggTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑aaa                                                                                       <  1:4883635/93‑1 (MQ=255)

   gTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑aaaa                                                                                      >  2:4069394/1‑93 (MQ=255)

   gTAAAGAAGCCGCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑aaaa                                                                                      <  1:3467497/93‑1 (MQ=255)

     aaaGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑aaaaac                                                                                    >  1:1481207/1‑91 (MQ=255)

     aaaGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGGGGGTAATGTT‑aaaaac                                                                                    >  1:6242307/1‑91 (MQ=255)

      aagaagCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑aaaaacc                                                                                   <  2:8597537/93‑4 (MQ=255)

      aagaagCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑aaaaacc                                                                                   <  1:6208924/93‑4 (MQ=255)

      aagaagCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑aaaaacc                                                                                   >  1:7110707/1‑90 (MQ=255)

      aagaagCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑aaaaacc                                                                                   >  1:3654088/1‑90 (MQ=255)

      aagaagCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑aaaaacc                                                                                   >  2:4556655/1‑90 (MQ=255)

       agaagCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑aaaaacca                                                                                  <  1:4453255/93‑5 (MQ=255)

       agaagCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑aaaaacca                                                                                  <  2:1008201/93‑5 (MQ=255)

          agCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAggg                                                                               <  2:650992/93‑1 (MQ=255)

                 gAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAg                                                                        >  1:1202200/1‑93 (MQ=255)

                 gAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAg                                                                        >  1:1116605/1‑93 (MQ=255)

                 gAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAg                                                                        >  2:3783141/1‑93 (MQ=255)

                    aTCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTa                                                                     <  2:4745749/93‑1 (MQ=255)

                    aTCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTa                                                                     >  1:4308372/1‑93 (MQ=255)

                     tCGCGGCTAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTaa                                                                    <  2:8356750/93‑1 (MQ=255)

                     aggcgGCCAGTTTCTGTGTTCCAGCAAACCATTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑AAAAGCCAGGGCGCGGAGTTaa                                                                    <  2:7949061/91‑1 (MQ=255)

                       gcgGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAAtt                                                                  >  1:6660136/1‑93 (MQ=255)

                         ggCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCa                                                                >  1:9499667/1‑93 (MQ=255)

                          gCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCaa                                                               <  1:816560/93‑1 (MQ=255)

                            cAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAAcc                                                             <  1:5861355/93‑1 (MQ=255)

                            cAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAAcc                                                             >  1:4033421/1‑93 (MQ=255)

                               tGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGt                                                          >  1:7984159/1‑93 (MQ=255)

                               tGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGt                                                          >  1:7495003/1‑93 (MQ=255)

                               tGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGt                                                          >  1:1152975/1‑93 (MQ=255)

                                    gtgtTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGAt                                                     >  2:2439372/1‑93 (MQ=255)

                                     tgtTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATa                                                    <  1:8822606/93‑1 (MQ=255)

                                      gtTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATaa                                                   <  2:2748912/93‑1 (MQ=255)

                                       ttCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAg                                                  >  1:1807316/1‑93 (MQ=255)

                                        tCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGc                                                 <  1:8666073/93‑1 (MQ=255)

                                         ccAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGcc                                                >  2:1902108/1‑93 (MQ=255)

                                         ccAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGcc                                                <  2:5965849/93‑1 (MQ=255)

                                         ccAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGcc                                                <  1:7064238/93‑1 (MQ=255)

                                            gCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTgg                                             <  1:5088381/93‑1 (MQ=255)

                                             cAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGga                                            <  2:6278119/93‑1 (MQ=255)

                                             cAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGga                                            <  1:92002/93‑1 (MQ=255)

                                             cAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGga                                            <  1:3031296/93‑1 (MQ=255)

                                               aaCCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGgaga                                          >  2:7916713/1‑93 (MQ=255)

                                                aCCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGat                                         >  2:353562/1‑93 (MQ=255)

                                                aCCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGat                                         <  2:7061001/93‑1 (MQ=255)

                                                 ccAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGatt                                        <  2:3269259/93‑1 (MQ=255)

                                                 ccAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGatt                                        >  1:830344/1‑93 (MQ=255)

                                                  cAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGatta                                       <  1:5253075/93‑1 (MQ=255)

                                                   aGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGattat                                      >  1:7597944/1‑93 (MQ=255)

                                                    ggtctcCTCGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGattatt                                     <  2:1590209/91‑1 (MQ=255)

                                                    gTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGattatt                                     >  1:2738688/1‑93 (MQ=255)

                                                     ttCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTg                                    >  1:9115923/1‑93 (MQ=255)

                                                     ttCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTg                                    <  1:2191877/93‑1 (MQ=255)

                                                     gtctcCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTg                                    <  2:6777306/92‑1 (MQ=255)

                                                      tctcCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGTTTATTGa                                   >  2:9262242/1‑93 (MQ=255)

                                                       ctcCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGAt                                  <  2:9052739/93‑1 (MQ=255)

                                                          cTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGA‑AAGCCTGGAGATTATTGATcaaa                              >  2:7173842/1‑92 (MQ=255)

                                                          cTGGCACCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGATAATTCAACCTGTTTGATAAGCCTAGAGATTATTAATcaa                               >  1:8458226/1‑93 (MQ=255)

                                                               aCCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATa                          <  2:6674001/93‑1 (MQ=255)

                                                               aCCCGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATa                          >  1:1961698/1‑93 (MQ=255)

                                                                cccGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATAc                         >  1:325577/1‑93 (MQ=255)

                                                                 ccGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACt                        >  2:8157534/1‑93 (MQ=255)

                                                                  cGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTg                       <  1:9347295/93‑1 (MQ=255)

                                                                  cGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTg                       >  1:5518681/1‑93 (MQ=255)

                                                                  cGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTg                       >  2:3503187/1‑93 (MQ=255)

                                                                  cGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTg                       >  2:3503184/1‑93 (MQ=255)

                                                                  cGGTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTg                       >  1:4586318/1‑93 (MQ=255)

                                                                   ggTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGc                      >  1:1335409/1‑93 (MQ=255)

                                                                   ggTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGc                      >  1:5608292/1‑93 (MQ=255)

                                                                    gTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCg                     >  2:2293191/1‑93 (MQ=255)

                                                                    gTATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCg                     >  1:3595626/1‑93 (MQ=255)

                                                                     tATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGt                    <  1:8586890/93‑1 (MQ=255)

                                                                     tATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGt                    >  1:2843638/1‑93 (MQ=255)

                                                                     tATCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGt                    <  1:2342658/93‑1 (MQ=255)

                                                                      aTCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTc                   <  1:1663694/93‑1 (MQ=255)

                                                                       tCGCGCGCGGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCa                  >  2:6943175/1‑93 (MQ=255)

                                                                        cgcgcgcgGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCAt                 <  1:1370120/93‑1 (MQ=255)

                                                                         gcgcgcgGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATa                >  1:6893278/1‑93 (MQ=255)

                                                                             gcgGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGaa            <  2:8905307/93‑1 (MQ=255)

                                                                              cgGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGaaa           >  2:9900893/1‑93 (MQ=255)

                                                                              cgGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGaaa           <  2:4468420/93‑1 (MQ=255)

                                                                              cgGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGaaa           <  2:9507162/93‑1 (MQ=255)

                                                                              cgGTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGaaa           >  2:2175175/1‑93 (MQ=255)

                                                                               ggTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAAt          <  1:5581267/93‑1 (MQ=255)

                                                                               ggTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAAt          >  2:8775591/1‑93 (MQ=255)

                                                                                gTGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAAtc         <  2:9020251/93‑1 (MQ=255)

                                                                                 tGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAAtct        >  2:5434228/1‑93 (MQ=255)

                                                                                 tGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAAtct        >  1:657602/1‑93 (MQ=255)

                                                                                 tGGGTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAAtct        <  2:1537224/93‑1 (MQ=255)

                                                                                  gtgTAATCTT‑AAAACCCAGGGCGCGTAGTTAATTCAACCTGTTTGATAAGCCTGGATATTATTGATCAACAATACTGCGTCATAAGAAAtctc       <  2:4941615/91‑1 (MQ=255)

                                                                                    gTAATGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTa     >  2:8065217/1‑93 (MQ=255)

                                                                                      aaTGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAGCCTGTTTGATAAGCCTGGAGATTATTGATCAACAATCCTGCGTCATAAGAAATCTCTAtt   <  1:3038141/93‑1 (MQ=255)

                                                                                      aaTGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTAtt   <  2:1675636/93‑1 (MQ=255)

                                                                                      aaTGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTAtt   <  2:8450805/93‑1 (MQ=255)

                                                                                       aTGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGCAATCTCTATTa  >  1:7155315/1‑93 (MQ=255)

                                                                                       aTGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTATTa  <  2:7399170/93‑1 (MQ=255)

                                                                                       aTGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTATTa  <  2:4104843/93‑1 (MQ=255)

                                                                                       aTGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTATTa  >  2:201940/1‑93 (MQ=255)

                                                                                       aTGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTATTa  >  1:4663359/1‑93 (MQ=255)

                                                                                       aTGTTAAAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTATTa  <  1:6058749/93‑1 (MQ=255)

GCGGTAAAGAAGCCTCAGAAATCGCGGCCAGTGGCTGTGTTCCAGCAAACCAGTTCTCCTGGCACCCGGTATCGCGCGCGGTGGGTAATGTT‑AAAA‑CCAGGGCGCGGAGTTAATTCAACCTGTTTGATAAGCCTGGAGATTATTGATCAACAATACTGCGTCATAAGAAATCTCTATTA  >  NC_000913/2008385‑2008563

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