BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	916,439	T→C	D248D (GAT→GAC)	ybjD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	916,439	0	T	C	99.5%	684.0 / ‑5.7	218	D248D (GAT→GAC)	ybjD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (105/112); total (105/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 226 total aligned reads displayed.

TTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCAAAAGGC  >  NC_000913/916349‑916529

ttattTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGacat                                                                                          >  2:6653032/1‑93 (MQ=255)

ttattTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGacat                                                                                          >  2:7307191/1‑93 (MQ=255)

       tctGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAAc                                                                                   >  1:1421495/1‑93 (MQ=255)

       tctGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAAc                                                                                   >  2:7541303/1‑93 (MQ=255)

        ctGAGCAGGGGGCCGGACAGGCGCGATATCTTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAAcc                                                                                  <  2:1521268/93‑1 (MQ=255)

         tGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCg                                                                                 >  2:5347731/1‑93 (MQ=255)

           aGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGga                                                                               >  1:7682084/1‑93 (MQ=255)

              aGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGgatga                                                                            <  1:7117718/93‑1 (MQ=255)

              aGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGgatga                                                                            <  2:5090540/93‑1 (MQ=255)

              aGGGGGCCGAACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGgatga                                                                            <  1:174026/93‑1 (MQ=255)

               gggggCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGgatgat                                                                           >  1:5639073/1‑93 (MQ=255)

               gggggCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGgatgat                                                                           >  2:7090608/1‑93 (MQ=255)

                ggggCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGAtt                                                                          <  2:6845164/93‑1 (MQ=255)

                ggggCCGAACAGGCGCGAGATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGAtt                                                                          <  1:2270592/93‑1 (MQ=255)

                 gggCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTg                                                                         <  2:2827220/93‑1 (MQ=255)

                   gCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTgac                                                                       <  1:4746186/93‑1 (MQ=255)

                   gCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTgac                                                                       <  1:1325786/93‑1 (MQ=255)

                    ccGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTgacc                                                                      <  2:166203/93‑1 (MQ=255)

                      ggACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTgaccga                                                                    >  1:798727/1‑93 (MQ=255)

                       gACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTgaccgac                                                                   >  2:315937/1‑93 (MQ=255)

                       gACAGGCGCGATATCGTTTAATGCGCCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTgaccgac                                                                   <  2:814189/93‑1 (MQ=255)

                        aCAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGTATGATTgaccgacc                                                                  >  1:2419050/1‑93 (MQ=255)

                        aCAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTgaccgacc                                                                  <  2:5675870/93‑1 (MQ=255)

                        aCAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTgaccgacc                                                                  <  1:2653013/93‑1 (MQ=255)

                         cAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCt                                                                 <  1:1042439/93‑1 (MQ=255)

                           ggCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTAACCGACCtgg                                                               >  2:5955887/1‑93 (MQ=255)

                            gcgcGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCtggt                                                              >  2:1586995/1‑93 (MQ=255)

                             cgcgATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCtggtg                                                             <  1:4537141/93‑1 (MQ=255)

                             cgcgATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCtggtg                                                             <  1:5138509/93‑1 (MQ=255)

                             cgcgATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCtggtg                                                             <  2:1251073/93‑1 (MQ=255)

                                   atCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCt                                                       >  1:2835957/1‑93 (MQ=255)

                                   atCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCt                                                       >  2:7335265/1‑93 (MQ=255)

                                   atCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCt                                                       <  1:6326465/93‑1 (MQ=255)

                                      gTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGt                                                    >  2:3784918/1‑93 (MQ=255)

                                      gTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGt                                                    >  1:2552087/1‑93 (MQ=255)

                                      gTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGt                                                    >  1:7913056/1‑93 (MQ=255)

                                      gTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGt                                                    >  1:7767677/1‑93 (MQ=255)

                                       tttAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTa                                                   >  1:2053971/1‑93 (MQ=255)

                                        ttAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTAt                                                  <  1:1112950/93‑1 (MQ=255)

                                         tAATGCGGCGGCGAGCCAGCAATGGGCAACGAAGCTGGCGCTATCTGGACATCATCACCCGGATGATTGACCGACCTGGTGGGCGCTCGTATc                                                 <  2:2479790/93‑1 (MQ=255)

                                         tAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATc                                                 >  1:7398236/1‑93 (MQ=255)

                                         tAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATc                                                 <  1:3359029/93‑1 (MQ=255)

                                           aTGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCgg                                               >  1:7554248/1‑93 (MQ=255)

                                               ggcggcgAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTa                                           <  2:3750091/93‑1 (MQ=255)

                                               ggcggcgAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTa                                           <  2:4506978/93‑1 (MQ=255)

                                               ggcggcgAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTa                                           <  2:5740155/93‑1 (MQ=255)

                                                 cggcgAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTAtt                                         <  1:2381277/93‑1 (MQ=255)

                                                 cggcgAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTAtt                                         <  2:2758832/93‑1 (MQ=255)

                                                  ggcgAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTAttt                                        <  2:6124756/93‑1 (MQ=255)

                                                  ggcgAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTAttt                                        >  2:2728824/1‑93 (MQ=255)

                                                   gcgAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTAtttt                                       >  2:6437344/1‑93 (MQ=255)

                                                   gcgAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTAtttt                                       <  2:3626644/93‑1 (MQ=255)

                                                   gcgAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTAtttt                                       <  2:669999/93‑1 (MQ=255)

                                                   gcgAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTAtttt                                       <  1:2450289/93‑1 (MQ=255)

                                                   gcgAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTAtttt                                       <  2:7657674/93‑1 (MQ=255)

                                                      aGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCt                                    <  1:6979265/93‑1 (MQ=255)

                                                       gCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTc                                   <  1:5308486/93‑1 (MQ=255)

                                                        ccAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCg                                  <  2:1503040/93‑1 (MQ=255)

                                                        ccAGCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCg                                  <  1:4576053/93‑1 (MQ=255)

                                                           gCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGcc                               >  1:2680716/1‑93 (MQ=255)

                                                           gCAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGcc                               >  2:538700/1‑93 (MQ=255)

                                                            cAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCt                              >  1:835357/1‑93 (MQ=255)

                                                            cAATGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCt                              >  1:3719729/1‑93 (MQ=255)

                                                             aaTGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTa                             <  2:5078220/93‑1 (MQ=255)

                                                              aTGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTAt                            >  1:1771513/1‑93 (MQ=255)

                                                              aTGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTAt                            <  2:328168/93‑1 (MQ=255)

                                                              aTGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCAACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTAt                            <  2:1187433/93‑1 (MQ=255)

                                                               tGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTAtt                           <  2:217552/93‑1 (MQ=255)

                                                               tGAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTAtt                           >  2:2263885/1‑93 (MQ=255)

                                                                gAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTAttt                          <  2:2319208/93‑1 (MQ=255)

                                                                gAGCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTAttt                          <  2:4544000/93‑1 (MQ=255)

                                                                  gCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTgt                        >  2:6174305/1‑92 (MQ=255)

                                                                  gCAACGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGc                        >  2:4114955/1‑93 (MQ=255)

                                                                   cAACGAAGCTGGCGCTATCTGGACATCATCAACCGAATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCt                       >  2:7799817/1‑93 (MQ=255)

                                                                    aaCGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCtc                      >  1:4241664/1‑92 (MQ=255)

                                                                    aaCGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTa                      >  2:2923244/1‑93 (MQ=255)

                                                                     aCGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTAc                     >  2:1570394/1‑93 (MQ=255)

                                                                     aCGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTAc                     >  1:582292/1‑93 (MQ=255)

                                                                      cGAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACt                    <  1:5940681/93‑1 (MQ=255)

                                                                       gAAGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACtt                   >  1:311192/1‑93 (MQ=255)

                                                                        aaGCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACttt                  <  2:4219229/93‑1 (MQ=255)

                                                                          gCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTttgt                <  1:1829172/93‑1 (MQ=255)

                                                                          gCTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTttgt                <  2:3467090/93‑1 (MQ=255)

                                                                           cTGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTttgtt               <  1:6815481/93‑1 (MQ=255)

                                                                            tGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGTCCTATTTGCTACTttgttg              >  1:5697095/1‑93 (MQ=255)

                                                                            tGGCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTttgttg              >  2:1915726/1‑93 (MQ=255)

                                                                             ggCGCTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGc             <  1:370269/93‑1 (MQ=255)

                                                                              gcgcTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCa            <  1:6291030/93‑1 (MQ=255)

                                                                              gcgcTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCa            <  1:956826/93‑1 (MQ=255)

                                                                                gcTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAgg          <  1:2885777/93‑1 (MQ=255)

                                                                                 cTATCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGc         <  2:583322/93‑1 (MQ=255)

                                                                                   aTCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCaa       >  1:6533751/1‑93 (MQ=255)

                                                                                    tCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCaaa      >  2:5194754/1‑93 (MQ=255)

                                                                                    tCTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCaaa      <  1:258613/93‑1 (MQ=255)

                                                                                     cTGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCaaaa     <  1:5103189/93‑1 (MQ=255)

                                                                                      tGGATATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCAAAAg    <  1:3806939/93‑1 (MQ=255)

                                                                                      tGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCAAAAg    >  1:6087388/1‑93 (MQ=255)

                                                                                      tGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCAAAAg    >  2:5039313/1‑93 (MQ=255)

                                                                                      tGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCAAAAg    >  2:4525890/1‑93 (MQ=255)

                                                                                      tGGACATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCAAAAg    >  1:5367894/1‑93 (MQ=255)

TTATTTCTCTGAGCAGGGGGCCGGACAGGCGCGATATCGTTTAATGCGGCGGCGAGCCAGCAATGAGCAACGAAGCTGGCGCTATCTGGATATCATCAACCGGATGATTGACCGACCTGGTGGGCGCTCGTATCGGGTTATTTTGCTCGGCCTATTTGCTACTTTGTTGCAGGCAAAAGGC  >  NC_000913/916349‑916529

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