BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,613,804	C→T	S6S (AGC→AGT)	b1527 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,613,804	0	C	T	99.5%	744.8 / ‑5.6	208	S6S (AGC→AGT)	b1527	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/1); new base T (107/100); total (107/101)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 211 total aligned reads displayed.

TTTTTGATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGCATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGATGC  >  NC_000913/1613714‑1613893

 ttttGATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATc                                                                                        <  2:2520425/93‑1 (MQ=255)

 ttttGATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATc                                                                                        <  2:2147333/93‑1 (MQ=255)

 ttttGATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATc                                                                                        <  1:6883648/93‑1 (MQ=255)

   ttGATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAg                                                                                      <  2:1028558/93‑1 (MQ=255)

     gATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCa                                                                                    >  1:6635141/1‑93 (MQ=255)

     gATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCa                                                                                    >  1:2395089/1‑93 (MQ=255)

     gATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCa                                                                                    >  2:6970785/1‑93 (MQ=255)

     gATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCa                                                                                    >  1:7431965/1‑93 (MQ=255)

      aTAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCaa                                                                                   <  1:2905923/93‑1 (MQ=255)

      aTAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCaa                                                                                   >  2:5354918/1‑93 (MQ=255)

       tAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAAt                                                                                  >  1:1242556/1‑93 (MQ=255)

          ttCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTAc                                                                               >  2:2392134/1‑93 (MQ=255)

            cTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACat                                                                             <  2:4145193/93‑1 (MQ=255)

            cTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACat                                                                             <  1:1714054/93‑1 (MQ=255)

             ttGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACata                                                                            <  2:4409695/93‑1 (MQ=255)

             ttGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACata                                                                            <  1:2481010/93‑1 (MQ=255)

             ttGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACata                                                                            <  2:183995/93‑1 (MQ=255)

             ttGCGGTTAATATGCTCTATAGAGTGATGTTCCGATGACTTATTACTATATGGGGCAAATATGGTTACGCCAGTAAGCATCAGCAATTACata                                                                            <  1:5044284/93‑1 (MQ=255)

                cGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATAtct                                                                         <  2:4571463/93‑1 (MQ=255)

                cGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATAtct                                                                         <  1:3659768/93‑1 (MQ=255)

                 ggTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATAtctc                                                                        <  2:2541402/93‑1 (MQ=255)

                   ttAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCtt                                                                      <  1:946032/93‑1 (MQ=255)

                   ttAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCtt                                                                      >  1:7035504/1‑93 (MQ=255)

                   ttAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCtt                                                                      >  1:2614776/1‑93 (MQ=255)

                   ttAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCtt                                                                      >  1:5091875/1‑93 (MQ=255)

                   atAATATGCTCTATATAGTGATGTTCCGATTACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCtt                                                                      <  1:4003859/92‑1 (MQ=255)

                     aaTATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTcc                                                                    <  1:4289399/93‑1 (MQ=255)

                      atatGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCt                                                                   <  1:1964609/93‑1 (MQ=255)

                        atGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCtga                                                                 <  1:4560340/93‑1 (MQ=255)

                         tGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCtgat                                                                >  1:3945322/1‑93 (MQ=255)

                              tatataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGAttt                                                           <  2:4473154/93‑1 (MQ=255)

                              tatataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGAttt                                                           <  1:128071/93‑1 (MQ=255)

                              tatataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGAttt                                                           >  2:4065221/1‑93 (MQ=255)

                              tatataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGAttt                                                           >  1:6193984/1‑93 (MQ=255)

                               atataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGAtttt                                                          >  1:3186569/1‑93 (MQ=255)

                                tataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTc                                                         >  1:5397442/1‑93 (MQ=255)

                                 ataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTcc                                                        >  1:799758/1‑93 (MQ=255)

                                 ataGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTcc                                                        <  1:2527977/93‑1 (MQ=255)

                                   aGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCcg                                                      >  1:5838964/1‑93 (MQ=255)

                                    gTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCcgt                                                     >  2:7920960/1‑93 (MQ=255)

                                    gTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCcgt                                                     >  2:7438959/1‑93 (MQ=255)

                                     tGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCcgta                                                    <  2:1663757/93‑1 (MQ=255)

                                     tGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCcgta                                                    >  1:723272/1‑93 (MQ=255)

                                          ttCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTaa                                               <  2:4242818/93‑1 (MQ=255)

                                          ttCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTaa                                               >  2:7497727/1‑93 (MQ=255)

                                          ttCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTaa                                               >  1:4841630/1‑93 (MQ=255)

                                          ttCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTaa                                               <  1:477551/93‑1 (MQ=255)

                                            ccGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAata                                             <  2:510116/93‑1 (MQ=255)

                                             cGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAatat                                            >  1:4654553/1‑93 (MQ=255)

                                             cGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAatat                                            <  1:1529468/93‑1 (MQ=255)

                                              gATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATAtt                                           <  2:6242316/93‑1 (MQ=255)

                                              gATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATAtt                                           >  2:5361559/1‑93 (MQ=255)

                                              gATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATAtt                                           <  1:4157502/93‑1 (MQ=255)

                                               aTGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTg                                          <  2:3618132/93‑1 (MQ=255)

                                               aTGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTg                                          >  2:1409633/1‑93 (MQ=255)

                                               aTGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTg                                          >  1:2377398/1‑93 (MQ=255)

                                                tGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGc                                         <  1:1657815/93‑1 (MQ=255)

                                                tGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGc                                         >  1:4280907/1‑93 (MQ=255)

                                                tGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGc                                         >  1:4088220/1‑93 (MQ=255)

                                                 gACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCa                                        >  2:3502542/1‑93 (MQ=255)

                                                 gACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCa                                        <  2:5011394/93‑1 (MQ=255)

                                                  aCTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCAc                                       <  2:6478598/93‑1 (MQ=255)

                                                    ttATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCg                                     <  2:5409227/93‑1 (MQ=255)

                                                     tATGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTGCATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGc                                    >  1:5264061/1‑93 (MQ=255)

                                                       tGACTATATGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAg                                  <  2:3128821/93‑1 (MQ=255)

                                                           tataTGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTaa                              >  2:2595392/1‑93 (MQ=255)

                                                           tataTGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTaa                              >  1:7757178/1‑93 (MQ=255)

                                                           tataTGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTaa                              >  1:3652808/1‑93 (MQ=255)

                                                           tataTGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGAACCGCAGGTaa                              >  2:1709543/1‑93 (MQ=255)

                                                             tatGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTaaaa                            >  1:1368444/1‑93 (MQ=255)

                                                             tatGGGGCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTaaaa                            <  1:1111924/93‑1 (MQ=255)

                                                                ggggCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAaagaa                         >  1:2725148/1‑93 (MQ=255)

                                                                ggggCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAaagaa                         <  2:1624915/93‑1 (MQ=255)

                                                                ggggCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAaagaa                         <  2:5944575/93‑1 (MQ=255)

                                                                 gggCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAaagaag                        <  2:1778916/93‑1 (MQ=255)

                                                                 gggCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAaagaag                        <  2:5976158/93‑1 (MQ=255)

                                                                 gggCAAATATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAaagaag                        >  1:5882024/1‑93 (MQ=255)

                                                                     aaaTATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGtttt                    >  1:5755540/1‑93 (MQ=255)

                                                                      aaTATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTa                   >  2:4150166/1‑93 (MQ=255)

                                                                      aaTATGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTa                   >  1:668490/1‑93 (MQ=255)

                                                                          tGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAgg               >  1:7807533/1‑93 (MQ=255)

                                                                          tGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAgg               >  2:6802137/1‑93 (MQ=255)

                                                                          tGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAgg               <  2:3894961/93‑1 (MQ=255)

                                                                          tGGTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAgg               <  2:2954931/93‑1 (MQ=255)

                                                                           ggTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGa              <  1:4251323/93‑1 (MQ=255)

                                                                            gTTACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGAt             >  1:4929052/1‑93 (MQ=255)

                                                                             ttACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGAtt            <  2:813655/93‑1 (MQ=255)

                                                                              tACGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGAttt           <  1:1444386/93‑1 (MQ=255)

                                                                                cGCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTa         <  2:477020/93‑1 (MQ=255)

                                                                                 gCCAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTAt        <  1:871478/93‑1 (MQ=255)

                                                                                  ccAGTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTAtt       <  1:6236951/93‑1 (MQ=255)

                                                                                    agtaagtaTCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGa     <  2:2851385/93‑1 (MQ=255)

                                                                                     gTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGAt    <  1:6012035/93‑1 (MQ=255)

                                                                                     gTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGAt    >  1:216416/1‑93 (MQ=255)

                                                                                     gTAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGAt    >  2:7821521/1‑93 (MQ=255)

                                                                                      tAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGATg   >  1:2004030/1‑93 (MQ=255)

                                                                                      tAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGATg   >  1:7272385/1‑93 (MQ=255)

                                                                                      tAAGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGATg   >  2:294519/1‑93 (MQ=255)

                                                                                       aaGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGATGc  <  2:691928/93‑1 (MQ=255)

                                                                                       aaGTATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGATGc  <  2:1225653/93‑1 (MQ=255)

TTTTTGATAGTTCTTGCGGTTAATATGCTCTATATAGTGATGTTCCGATGACTTATGACTATATGGGGCAAATATGGTTACGCCAGTAAGCATCAGCAATTACATATCTCTTCCTGATGATTTTCCCGTACGTAATATTGCACCGCAGGTAAAAGAAGTTTTAAAGGATTTTATTGATGC  >  NC_000913/1613714‑1613893

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