BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,938,723	C→T	R28R (CGG→CGA)	yiePs ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,938,723	0	C	T	99.5%	1520.8 / ‑4.8	417	R28R (CGG→CGA)	yiePs	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (2/0); new base T (239/176); total (241/176)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

Only 100 of 425 total aligned reads displayed.

GTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAA  >  NC_000913/3938633‑3938812

 tGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGa                                                                                        <  2:4813230/93‑1 (MQ=255)

  gATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGCAAGCAATTCCTGATCAAGAACATTCCAGTTCGATTGTGGCATGACTCGAg                                                                                       >  1:80290/1‑93 (MQ=255)

    tCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTa                                                                                     >  2:1933531/1‑93 (MQ=255)

    tCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCATGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTa                                                                                     <  2:4536329/93‑1 (MQ=255)

     cGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTAc                                                                                    <  2:5979508/93‑1 (MQ=255)

      gATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGCCTCGAGGAcc                                                                                   >  1:379827/1‑93 (MQ=255)

      gATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTAcc                                                                                   <  1:7553344/93‑1 (MQ=255)

        tGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCaa                                                                                 <  1:430231/93‑1 (MQ=255)

        tGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCaa                                                                                 <  2:226539/93‑1 (MQ=255)

          aCCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAAtt                                                                               <  2:7219899/93‑1 (MQ=255)

          aCCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAAtt                                                                               >  1:4941150/1‑93 (MQ=255)

           ccTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtc                                                                              <  1:1204578/93‑1 (MQ=255)

            cTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcg                                                                             >  1:4826226/1‑93 (MQ=255)

             tgatgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgc                                                                            <  1:7241112/93‑1 (MQ=255)

             tgatgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgc                                                                            <  2:2852299/93‑1 (MQ=255)

             tgatgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACCCGAGTACCAATtcgc                                                                            >  1:5550472/1‑93 (MQ=255)

               atgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcgg                                                                          >  2:3607725/1‑93 (MQ=255)

               atgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcgg                                                                          >  1:295921/1‑93 (MQ=255)

               atgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcgg                                                                          >  2:5790827/1‑93 (MQ=255)

                tgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcggt                                                                         >  1:5885669/1‑93 (MQ=255)

                tgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcggt                                                                         <  1:1034092/93‑1 (MQ=255)

                tgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCACGACTCGAGTACCAATtcgcggt                                                                         <  2:4977973/93‑1 (MQ=255)

                  aaaGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcggtcg                                                                       >  1:7760405/1‑93 (MQ=255)

                    aGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcggtcgcg                                                                     <  2:1688238/93‑1 (MQ=255)

                      ttCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTcgcggtcgcggt                                                                   <  2:7436996/93‑1 (MQ=255)

                          cttctGTCATCCACCAGTTAAGCAATTCCGGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTaaaa                                                               <  1:6005977/93‑1 (MQ=255)

                          cttctGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTaaaa                                                               >  1:6317922/1‑93 (MQ=255)

                          cttctGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTaaaa                                                               <  2:5416833/93‑1 (MQ=255)

                           ttctGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAAc                                                              >  2:3501805/1‑93 (MQ=255)

                           ttctGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAAc                                                              >  2:5665634/1‑93 (MQ=255)

                             ctGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCa                                                            >  2:616650/1‑93 (MQ=255)

                              tGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCAt                                                           <  2:6073847/93‑1 (MQ=255)

                              tGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCAt                                                           >  2:1796144/1‑93 (MQ=255)

                              tGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCAt                                                           <  1:4554864/93‑1 (MQ=255)

                               gTCATCCACCAGGTACGCACTTCCTGATCAAGATAATTCCAGTTCGATTGTGGCATGACTCGAATCCCAATTCGCGGTCGCGGTAAAACCATc                                                          <  1:7041370/93‑1 (MQ=255)

                                tCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATcc                                                         >  1:2619430/1‑93 (MQ=255)

                                 cATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATccc                                                        <  1:4372725/93‑1 (MQ=255)

                                 cATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATccc                                                        <  2:6664886/93‑1 (MQ=255)

                                 cATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATccc                                                        >  2:2722864/1‑93 (MQ=255)

                                  aTCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCt                                                       >  1:3334340/1‑93 (MQ=255)

                                  aTCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCt                                                       >  2:3747823/1‑93 (MQ=255)

                                     caccaGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCtttt                                                    <  2:3059854/93‑1 (MQ=255)

                                     caccaGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCtttt                                                    >  1:2212664/1‑93 (MQ=255)

                                       ccaGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGc                                                  >  1:943243/1‑93 (MQ=255)

                                       ccaGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGc                                                  >  1:5114237/1‑93 (MQ=255)

                                        caGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGcc                                                 >  2:4941492/1‑93 (MQ=255)

                                        caGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGcc                                                 <  1:5933830/93‑1 (MQ=255)

                                        caGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGcc                                                 >  2:251284/1‑93 (MQ=255)

                                         aGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCg                                                >  2:3240046/1‑93 (MQ=255)

                                             aaGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTaa                                            <  1:4047806/93‑1 (MQ=255)

                                              aGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTacc                                           >  1:6829390/1‑91 (MQ=255)

                                              aGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAAc                                           >  1:1719226/1‑93 (MQ=255)

                                              aGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAAc                                           >  1:2055875/1‑93 (MQ=255)

                                              aGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAAc                                           >  2:2051615/1‑93 (MQ=255)

                                               gCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACg                                          >  1:1857929/1‑93 (MQ=255)

                                               gCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACg                                          <  2:7338954/93‑1 (MQ=255)

                                                 aaTTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGtt                                        >  2:3388918/1‑93 (MQ=255)

                                                 aaTTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGtt                                        >  2:7133569/1‑93 (MQ=255)

                                                 aaTTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGtt                                        >  2:5239803/1‑93 (MQ=255)

                                                  aTTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGGGAACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGCTAACGttt                                       >  2:64773/1‑93 (MQ=255)

                                                  aTTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGttt                                       >  2:2947928/1‑93 (MQ=255)

                                                   ttCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGtttt                                      >  1:6187864/1‑93 (MQ=255)

                                                    tCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTg                                     >  1:2799300/1‑93 (MQ=255)

                                                       tGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGAcc                                  <  2:7598437/93‑1 (MQ=255)

                                                       tGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGAcc                                  >  2:5514280/1‑93 (MQ=255)

                                                       tGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGAcc                                  >  1:1116275/1‑93 (MQ=255)

                                                       tGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGAcc                                  >  1:7819233/1‑93 (MQ=255)

                                                          tCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCt                               <  1:3474111/93‑1 (MQ=255)

                                                           cAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCtt                              >  1:2397605/1‑93 (MQ=255)

                                                            aaGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTc                             <  1:4859796/93‑1 (MQ=255)

                                                            aaGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTc                             >  2:4253151/1‑93 (MQ=255)

                                                            aaGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTc                             >  1:6560370/1‑93 (MQ=255)

                                                             aGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTcg                            >  2:1258410/1‑93 (MQ=255)

                                                              gAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTcgc                           <  2:2416986/93‑1 (MQ=255)

                                                              gAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTcgc                           <  1:5632417/93‑1 (MQ=255)

                                                                 aaTTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTa                        <  1:5350661/93‑1 (MQ=255)

                                                                   ttCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTAcc                      <  1:3925509/93‑1 (MQ=255)

                                                                   ttCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTAcc                      <  2:5371011/93‑1 (MQ=255)

                                                                   ttCCAGTTCGATGGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTAcc                      <  2:1039070/93‑1 (MQ=255)

                                                                     ccAGTTCGCTTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGc                    <  1:3085566/93‑1 (MQ=255)

                                                                     ccAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGc                    <  2:4095155/93‑1 (MQ=255)

                                                                     ccAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGc                    <  1:1427502/93‑1 (MQ=255)

                                                                     ccAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGc                    >  2:2192356/1‑93 (MQ=255)

                                                                        gTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGt                 >  1:710879/1‑93 (MQ=255)

                                                                        gTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGt                 <  2:3972657/93‑1 (MQ=255)

                                                                        gTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGt                 <  1:3511808/93‑1 (MQ=255)

                                                                         ttCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTa                >  1:3573229/1‑93 (MQ=255)

                                                                          tCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTAc               <  1:5284155/93‑1 (MQ=255)

                                                                           cGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACg              >  2:6940962/1‑93 (MQ=255)

                                                                             aTTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGac            >  1:7379974/1‑93 (MQ=255)

                                                                              ttGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTACCGTTTTGACCGCTTCGCGTACCGCTGTCGGact           >  1:6680668/1‑93 (MQ=255)

                                                                                gtgGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGactca         <  2:1438071/93‑1 (MQ=255)

                                                                                 tgGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGactcac        <  1:4743758/93‑1 (MQ=255)

                                                                                   gCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTCCGactcactc      >  1:2653357/1‑93 (MQ=255)

                                                                                   gCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGactcactc      >  2:6349526/1‑93 (MQ=255)

                                                                                   gCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGactcactc      >  2:6987732/1‑93 (MQ=255)

                                                                                   gCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGactcactc      >  1:1254409/1‑93 (MQ=255)

                                                                                     aTGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCa    >  1:5084334/1‑93 (MQ=255)

                                                                                      tGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCaa   >  1:519353/1‑93 (MQ=255)

                                                                                       gACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCaaa  >  2:7783475/1‑93 (MQ=255)

GTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAA  >  NC_000913/3938633‑3938812

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