BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,472,339	(C)_6→7	coding (95/861 nt)	ydbC →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,472,333	1	.	C	81.8%	433.4 / 101.5	170	coding (89/861 nt)	ydbC	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base . (13/18); new base C (66/73); total (79/91)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

Only 100 of 189 total aligned reads displayed.

TGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑CCCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTA  >  NC_000913/1472241‑1472420

tGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTgg                                                                                           <  2:7802402/93‑1 (MQ=255)

  cGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑cc                                                                                       <  1:4816840/93‑1 (MQ=255)

   gATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑ccc                                                                                      >  1:2698188/1‑93 (MQ=255)

   gATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGGTTTTGG‑‑cgc                                                                                      >  1:6638064/1‑91 (MQ=255)

    aTGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑cccc                                                                                     >  1:4484993/1‑93 (MQ=255)

      gAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑cccccc                                                                                   <  2:2170730/93‑1 (MQ=255)

      gAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑cccccc                                                                                   <  1:6334521/93‑1 (MQ=255)

        gcagcaATACATTTACTCTCGGTACAATATCCGTTAACCGTCTTGGTTATGGCGCGATGCAAGTGGCAGGTCCTGGAGTTTTTGG‑‑ccccccca                                                                                 <  2:2229737/93‑3 (MQ=255)

        gcagcaATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑ccccccca                                                                                 >  1:4106588/1‑91 (MQ=255)

        gcagcaATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑CCCCCCAc                                                                                 <  1:3564079/93‑1 (MQ=255)

         cagcaATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑cccccccac                                                                                >  1:4899407/1‑90 (MQ=255)

            caATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACgag                                                                              >  1:4227225/1‑93 (MQ=255)

             aaTACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑CCCCCCCCGAGAt                                                                            >  2:5707839/1‑93 (MQ=255)

              aTACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGAt                                                                            >  2:5366108/1‑93 (MQ=255)

              aTACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGAt                                                                            <  1:191860/93‑1 (MQ=255)

              aTACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGAt                                                                            <  1:6910409/93‑1 (MQ=255)

              aTACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGAt                                                                            >  1:5757633/1‑93 (MQ=255)

               tACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATc                                                                           >  1:6484233/1‑93 (MQ=255)

               tACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑CCCCCCACGAGATCg                                                                          >  1:7503155/1‑93 (MQ=255)

                aCATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCg                                                                          >  2:6116941/1‑93 (MQ=255)

                  aTTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGcc                                                                        >  1:3362270/1‑93 (MQ=255)

                   tttACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGCCCCCCCCACGAGATCGcc                                                                        <  1:6921201/93‑1 (MQ=255)

                    ttACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGccgc                                                                      <  1:3255119/93‑3 (MQ=255)

                      aCTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGt                                                                    >  2:3958667/1‑93 (MQ=255)

                       ctctCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTc                                                                   >  2:5607772/1‑93 (MQ=255)

                        tctcGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCg                                                                  >  2:7524577/1‑93 (MQ=255)

                         ctcGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGc                                                                 >  2:7422588/1‑93 (MQ=255)

                             gTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCtata                                                             >  2:2453898/1‑93 (MQ=255)

                             gTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCtata                                                             <  1:6088011/93‑1 (MQ=255)

                              tACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATaa                                                            <  2:5769249/93‑1 (MQ=255)

                               aCAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑CCCCCCACGAGATCGCCGCGTCGCTATAAcc                                                          <  2:6315594/93‑1 (MQ=255)

                                cAAAATACGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAAcc                                                          >  1:3979211/1‑93 (MQ=255)

                                 aaaaTCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCg                                                         <  2:7230028/93‑1 (MQ=255)

                                  aaaTCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGtg                                                       >  1:2137028/1‑93 (MQ=255)

                                     tCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGtgct                                                     >  1:5180758/1‑93 (MQ=255)

                                     tCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGtgct                                                     >  2:4616457/1‑93 (MQ=255)

                                     tCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGtgct                                                     <  2:4976976/93‑1 (MQ=255)

                                         ttAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGt                                                 <  1:1440776/93‑1 (MQ=255)

                                          tAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGa                                               <  1:7333595/93‑1 (MQ=255)

                                             ccGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAgg                                             >  2:3820149/1‑93 (MQ=255)

                                             ccGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAgg                                             <  2:5074717/93‑1 (MQ=255)

                                             ccGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAgg                                             >  2:1755245/1‑93 (MQ=255)

                                             ccGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGgc                                            <  1:5408598/93‑1 (MQ=255)

                                              cGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGgc                                            >  1:2342261/1‑93 (MQ=255)

                                                tCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGgcgc                                          <  1:3827075/93‑1 (MQ=255)

                                                tCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGgcgc                                          <  1:4384158/93‑1 (MQ=255)

                                                tCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCt                                         >  1:6752629/1‑93 (MQ=255)

                                                  ttGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTg                                        >  1:5250091/1‑93 (MQ=255)

                                                  ttGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTg                                        <  1:816209/93‑1 (MQ=255)

                                                   tGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTgg                                       <  1:4950239/93‑1 (MQ=255)

                                                   tGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTgg                                       <  2:3785649/93‑1 (MQ=255)

                                                   tGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTgg                                       <  2:3214764/93‑1 (MQ=255)

                                                     gTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCa                                     >  1:3236591/1‑93 (MQ=255)

                                                     gTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCa                                     >  1:81142/1‑93 (MQ=255)

                                                     gTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCAt                                    <  1:1856558/93‑1 (MQ=255)

                                                      ttATTGCGCGATGCAACCGGCAGGTCCTGGAGTTTTTGGCCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCa                                     <  1:2775664/93‑1 (MQ=255)

                                                      ttATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCAtt                                   <  2:4433020/93‑1 (MQ=255)

                                                       tATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGTGGCGCTGGCAtt                                   >  1:2540177/1‑93 (MQ=255)

                                                       tATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCAtt                                   <  1:984836/93‑1 (MQ=255)

                                                        aTGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTg                                  <  1:106021/93‑1 (MQ=255)

                                                         tGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTgg                                 <  1:870514/93‑1 (MQ=255)

                                                         tGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTgg                                 <  1:731146/93‑1 (MQ=255)

                                                         tGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTgg                                 <  2:6482001/93‑1 (MQ=255)

                                                          ggCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTggg                                <  1:3706571/93‑1 (MQ=255)

                                                           gcgcGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGc                               <  2:6777817/93‑1 (MQ=255)

                                                           gcgcGATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGc                               <  2:2872475/93‑1 (MQ=255)

                                                              cgATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTc                            <  2:5242863/93‑1 (MQ=255)

                                                               gATGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCa                           <  2:1132677/93‑1 (MQ=255)

                                                                aTGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCaa                          <  2:5794412/93‑1 (MQ=255)

                                                                aTGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCaa                          <  1:7118234/93‑1 (MQ=255)

                                                                 tGCAACTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAAt                         >  2:7014992/1‑93 (MQ=255)

                                                                    aaCTGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCat                      >  1:7523008/1‑93 (MQ=255)

                                                                    aaCTGGCAGGTCCTGGAGTTTTTGG‑‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCata                     <  1:364245/93‑1 (MQ=255)

                                                                       tGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATAtt                   <  2:7751301/93‑1 (MQ=255)

                                                                       tGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATAtt                   >  1:5953211/1‑93 (MQ=255)

                                                                       tGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATAtt                   <  2:1558666/93‑1 (MQ=255)

                                                                       tGGCAGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATAtt                   >  2:2698627/1‑93 (MQ=255)

                                                                         gCAGGTCCTTGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGa                 <  1:282597/93‑1 (MQ=255)

                                                                           aGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGAtc               >  1:2955030/1‑92 (MQ=255)

                                                                           aGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATa               <  1:6706051/93‑1 (MQ=255)

                                                                           aGGTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATa               <  2:6082280/93‑1 (MQ=255)

                                                                           aGGTCCTGGAGTGTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATa               <  2:4691500/93‑1 (MQ=255)

                                                                            ggTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATAc              <  1:7677828/93‑1 (MQ=255)

                                                                             gTCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATAcc             >  1:1707482/1‑93 (MQ=255)

                                                                             gTCCTGGAG‑TTTTGG‑‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAg           >  1:4744904/1‑93 (MQ=255)

                                                                              tCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCa            <  2:5336529/93‑1 (MQ=255)

                                                                              tCCTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCa            >  1:4833114/1‑93 (MQ=255)

                                                                               ccTGGAGTTTTTGG‑‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGc          <  1:245367/93‑1 (MQ=255)

                                                                                cTGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGc          >  2:182417/1‑93 (MQ=255)

                                                                                 tGGAGTTTTTGGCCCCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGc          >  1:6367445/1‑93 (MQ=255)

                                                                                 tGGAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCg         >  2:1234932/1‑93 (MQ=255)

                                                                                 tGGAGTTTTTGG‑‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGa        >  2:3586361/1‑93 (MQ=255)

                                                                                   gAGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGAc       >  2:2321415/1‑93 (MQ=255)

                                                                                    aGTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACt      >  1:7614158/1‑93 (MQ=255)

                                                                                     gTTTTTGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACtt     <  2:4064316/93‑1 (MQ=255)

                                                                                      tttttGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACttt    >  2:4153098/1‑93 (MQ=255)

                                                                                      tttttGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACttt    <  1:3768189/93‑1 (MQ=255)

                                                                                       ttttGG‑‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTa  <  1:591730/93‑1 (MQ=255)

                                                                                        tttGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTa  <  1:4723642/93‑1 (MQ=255)

                                                                                        tttGGC‑CCCCCCACGAGATCGCCGCGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTa  <  2:551397/93‑1 (MQ=255)

TGCGATGAGCAGCAATACATTTACTCTCGGTACAAAATCCGTTAACCGTCTTGGTTATGGCGCGATGCAACTGGCAGGTCCTGGAGTTTTTGG‑‑CCCCCCACGAGATCGCCACGTCGCTATAACCGTGCTGCGTGAGGCGCTGGCATTGGGCGTCAATCATATTGATACCAGCGACTTTTA  >  NC_000913/1472241‑1472420

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