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.2%	890.7 / ‑4.6	245	R28R (CGG→CGA)	yiePs	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (2/0); new base T (143/100); total (145/100)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

Only 100 of 249 total aligned reads displayed.

GTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAA  >  NC_000913/3938633‑3938812

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

  gATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAg                                                                                       <  2:2016502/93‑1 (MQ=255)

      gATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAATTTCCAGTTCGATTGTGGCATGACTCGAGTAcc                                                                                   >  1:137271/1‑93 (MQ=255)

      gATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTAcc                                                                                   >  1:3181744/1‑93 (MQ=255)

      gATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTAcc                                                                                   >  1:2090958/1‑93 (MQ=255)

       aTGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCa                                                                                  >  2:4586878/1‑93 (MQ=255)

        tGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCaa                                                                                 >  1:1952877/1‑93 (MQ=255)

         gACCTGATGAAAGTTCTCTTCTGTCATCAACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAAt                                                                                <  1:4509010/93‑1 (MQ=255)

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

           ccTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtc                                                                              >  2:1908878/1‑93 (MQ=255)

             tgatgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgc                                                                            >  2:3183356/1‑93 (MQ=255)

             tgatgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgc                                                                            >  1:4658257/1‑93 (MQ=255)

             tgatgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgc                                                                            >  1:4591533/1‑93 (MQ=255)

              gatgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcg                                                                           >  1:2716957/1‑93 (MQ=255)

              gatgaAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACAAATtcgcg                                                                           >  1:3588068/1‑93 (MQ=255)

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

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

                   aaGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcggtcgc                                                                      <  1:1928741/93‑1 (MQ=255)

                    aGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcggtcgcg                                                                     >  2:4382434/1‑93 (MQ=255)

                     gTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATtcgcggtcgcgg                                                                    <  2:2967839/93‑1 (MQ=255)

                      ttCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTcgcggtcgcggt                                                                   >  1:577653/1‑93 (MQ=255)

                       tctcTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGATTCTAGTACCAATTCGCGGTCGCGGTa                                                                  >  2:725202/1‑93 (MQ=255)

                       tctcTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTa                                                                  >  1:669302/1‑93 (MQ=255)

                        ctctTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTaa                                                                 <  1:2929120/93‑1 (MQ=255)

                        ctctTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTaa                                                                 <  1:927981/93‑1 (MQ=255)

                         tcttctGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTaaa                                                                >  2:2511359/1‑93 (MQ=255)

                           ttctGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAAc                                                              <  1:796346/93‑1 (MQ=255)

                            tctGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAAcc                                                             >  1:1702028/1‑93 (MQ=255)

                              tGTCATCCACCAGGTAAGCAATTCCTAATAAAGAAAATTCCAGTTCGATTGTGGCATGCCTCGAGTACCAATTCGCGGTCGCGGTAAAACCAt                                                           <  2:1878620/93‑1 (MQ=255)

                               gTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATc                                                          <  1:4196217/93‑1 (MQ=255)

                               gTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATc                                                          >  1:3999677/1‑93 (MQ=255)

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

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

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

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

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

                                      accaGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTg                                                   <  2:3807271/93‑1 (MQ=255)

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

                                        caGGTAAGCAATTCCTGATCAAGAAAATTCCCGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGcg                                                 >  1:4460395/1‑92 (MQ=255)

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

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

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

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

                                          ggTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGt                                               >  2:1515247/1‑93 (MQ=255)

                                          ggTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGt                                               >  1:2760092/1‑93 (MQ=255)

                                          ggTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGt                                               >  2:1412099/1‑93 (MQ=255)

                                           gTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGtt                                              <  2:2880195/93‑1 (MQ=255)

                                            tAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTa                                             <  1:3413327/93‑1 (MQ=255)

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

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

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

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

                                               gCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGAGGTAAAACCATCCCTTTTGCCGTTAACg                                          >  1:4233407/1‑93 (MQ=255)

                                                cAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGt                                         >  2:2437011/1‑93 (MQ=255)

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

                                                  aTTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGttt                                       <  2:3049649/93‑1 (MQ=255)

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

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

                                                   ttCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGtttt                                      <  1:949890/93‑1 (MQ=255)

                                                    tCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTg                                     <  1:284666/93‑1 (MQ=255)

                                                    tCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTg                                     <  2:2309104/93‑1 (MQ=255)

                                                     ccTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGa                                    <  1:4328109/93‑1 (MQ=255)

                                                     ccTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGa                                    >  1:464381/1‑93 (MQ=255)

                                                      cTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGAc                                   >  1:1351024/1‑93 (MQ=255)

                                                      cTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTCAc                                   >  1:3909566/1‑93 (MQ=255)

                                                       tGATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGaac                                  >  2:1226681/1‑91 (MQ=255)

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

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

                                                        gATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCg                                 >  2:3261878/1‑93 (MQ=255)

                                                        gATCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCg                                 >  1:3361507/1‑93 (MQ=255)

                                                          tCAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCt                               >  2:3095445/1‑93 (MQ=255)

                                                           cAAGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCtt                              <  2:3903858/93‑1 (MQ=255)

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

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

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

                                                             aGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTcg                            <  2:2365727/93‑1 (MQ=255)

                                                             aGAAAATTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTcg                            <  2:129218/93‑1 (MQ=255)

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

                                                               aaaaTTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTcgcg                          <  1:3727889/93‑1 (MQ=255)

                                                                aaaTTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGt                         >  2:4664274/1‑93 (MQ=255)

                                                                aaaTTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGt                         >  2:2572463/1‑93 (MQ=255)

                                                                 aaTTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTa                        >  2:1108503/1‑93 (MQ=255)

                                                                  aTTCCAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTAc                       >  2:2118426/1‑93 (MQ=255)

                                                                      cAGTTCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCt                   <  1:3798186/93‑1 (MQ=255)

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

                                                                         ttCGATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTa                <  1:1522196/93‑1 (MQ=255)

                                                                            gATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACgc             >  2:340367/1‑92 (MQ=255)

                                                                            gATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGa             >  2:3556070/1‑93 (MQ=255)

                                                                            gATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGa             >  2:4501032/1‑93 (MQ=255)

                                                                            gATTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGa             >  1:2226806/1‑93 (MQ=255)

                                                                             aTTGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGac            <  2:2245420/93‑1 (MQ=255)

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

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

                                                                              ttGTGGCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGact           >  2:2662860/1‑93 (MQ=255)

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

                                                                                  ggCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGactcact       <  1:1593312/93‑1 (MQ=255)

                                                                                  ggCATGACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGactcact       >  2:1677702/1‑93 (MQ=255)

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

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

                                                                                       gACTCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCaaa  <  1:1551745/93‑1 (MQ=255)

GTGATCGATGACCTGATGAAAGTTCTCTTCTGTCATCCACCAGGTAAGCAATTCCTGATCAAGAAAATTCCAGTTCGATTGTGGCATGACCCGAGTACCAATTCGCGGTCGCGGTAAAACCATCCCTTTTGCCGTTAACGTTTTGACCGCTTCGCGTACCGCTGTACGACTCACTCCAAA  >  NC_000913/3938633‑3938812

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