BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,132,550	T→C	G135G (GGT→GGC)	yghT →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,132,550	0	T	C	99.4%	1073.5 / ‑3.2	322	G135G (GGT→GGC)	yghT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (2/0); new base C (161/159); total (163/159)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

Only 100 of 367 total aligned reads displayed.

TCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAA  >  NC_000913/3132460‑3132639

tcgCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAggcgg                                                                                         >  1:579807/1‑93 (MQ=255)

 cgcgGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAggcggt                                                                                        >  1:4498646/1‑93 (MQ=255)

 cgcgGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAggcggt                                                                                        <  1:6867975/93‑1 (MQ=255)

  gcgGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTc                                                                                       <  1:6956411/93‑1 (MQ=255)

  gcgGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTc                                                                                       <  1:1825027/93‑1 (MQ=255)

    gggTTAATTGTTCTAACCGACCGTTATCCTCAGGACAAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGt                                                                                     >  2:4372980/1‑93 (MQ=255)

    gggTTAATTGTTCTAACCGACCGTTACCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGt                                                                                     <  2:6606128/93‑1 (MQ=255)

     ggTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGtt                                                                                    >  1:5994535/1‑93 (MQ=255)

      gTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGttt                                                                                   <  1:2261947/93‑1 (MQ=255)

      gTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGAGCTTACGATGGTACGGTGTTCCCCCCTAACGTTGAAGGCGGTCGttt                                                                                   <  1:3502020/93‑1 (MQ=255)

       ttAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGtttt                                                                                  <  1:1035381/93‑1 (MQ=255)

       ttAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGtttt                                                                                  <  1:509032/93‑1 (MQ=255)

         aaTTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGt                                                                                >  1:3130497/1‑93 (MQ=255)

         aaTTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGt                                                                                <  2:1102527/93‑1 (MQ=255)

          aTTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGtc                                                                               >  2:2981458/1‑93 (MQ=255)

          aTTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGtc                                                                               <  2:2857613/93‑1 (MQ=255)

          aTTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGtc                                                                               >  2:2594271/1‑93 (MQ=255)

          aTTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGtc                                                                               >  2:6471539/1‑93 (MQ=255)

           ttGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGtct                                                                              <  2:831625/93‑1 (MQ=255)

           ttGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGtct                                                                              <  1:7987420/93‑1 (MQ=255)

             gTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCa                                                                            <  2:4846986/93‑1 (MQ=255)

              ttCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCAt                                                                           <  1:4710529/93‑1 (MQ=255)

               tCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCAtg                                                                          <  2:8015946/93‑1 (MQ=255)

               tCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTGGTCTCAtg                                                                          >  1:4940373/1‑93 (MQ=255)

                cTAACCGACCGTTATCCTCAGGACCACATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCAtgg                                                                         <  1:7394731/93‑1 (MQ=255)

                cTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCAtgg                                                                         <  1:1096793/93‑1 (MQ=255)

                 tAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCAtggc                                                                        >  1:6708800/1‑93 (MQ=255)

                 tAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCAtggc                                                                        >  2:3460436/1‑93 (MQ=255)

                  aaCCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCAtggct                                                                       >  1:6286662/1‑93 (MQ=255)

                  aaCCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCAtggct                                                                       <  2:5666594/93‑1 (MQ=255)

                      gaccgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTGTCGTCTCATGGCTGGCa                                                                   >  2:7810604/1‑93 (MQ=255)

                       accgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCaa                                                                  <  2:4661105/93‑1 (MQ=255)

                       accgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCaa                                                                  <  1:8460904/93‑1 (MQ=255)

                       accgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCaa                                                                  >  2:867950/1‑93 (MQ=255)

                        ccgTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAg                                                                 <  1:8670153/93‑1 (MQ=255)

                          gTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGcc                                                               <  2:5032687/93‑1 (MQ=255)

                          gTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGcc                                                               >  2:2295794/1‑93 (MQ=255)

                          gTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGcc                                                               >  1:6538286/1‑93 (MQ=255)

                          gTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTGCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAggc                                                               >  1:411210/1‑91 (MQ=255)

                           ttATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCa                                                              <  2:7522718/93‑1 (MQ=255)

                           ttATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCa                                                              <  1:7078241/93‑1 (MQ=255)

                           ttATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCa                                                              >  1:1309691/1‑93 (MQ=255)

                            tATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAg                                                             <  1:2945398/93‑1 (MQ=255)

                            tATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAg                                                             <  2:5893938/93‑1 (MQ=255)

                             ttCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAgg                                                            <  1:3882294/92‑1 (MQ=255)

                             aTCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAgg                                                            >  2:278633/1‑93 (MQ=255)

                              tCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGa                                                           >  2:1232372/1‑93 (MQ=255)

                                  cAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGt                                                       >  2:1817978/1‑93 (MQ=255)

                                  cAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGt                                                       >  1:592851/1‑93 (MQ=255)

                                  cAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGt                                                       >  1:8113568/1‑93 (MQ=255)

                                   aGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTa                                                      >  1:3281876/1‑93 (MQ=255)

                                     gACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTaaa                                                    <  1:1931498/93‑1 (MQ=255)

                                      aCCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAg                                                   <  1:3528415/93‑1 (MQ=255)

                                        cAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCg                                                 >  2:627154/1‑93 (MQ=255)

                                         aaaTTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGt                                                <  2:6689480/93‑1 (MQ=255)

                                         aaaTTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGt                                                >  1:4329602/1‑93 (MQ=255)

                                            ttCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTc                                             <  1:7671941/93‑1 (MQ=255)

                                            ttCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTc                                             >  1:4835812/1‑93 (MQ=255)

                                            ttCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGCCAAGCCAGGAACCTAAAGCGTGTc                                             >  2:2449653/1‑93 (MQ=255)

                                              ccTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCAc                                           >  2:5294572/1‑93 (MQ=255)

                                              ccTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCAc                                           <  1:1562017/93‑1 (MQ=255)

                                               cTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACt                                          <  1:6077937/93‑1 (MQ=255)

                                                   cgcTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGAt                                      <  2:2043178/93‑1 (MQ=255)

                                                    gcTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATg                                     <  1:1147888/93‑1 (MQ=255)

                                                    gcTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATg                                     <  2:2579545/93‑1 (MQ=255)

                                                    gcTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATg                                     <  2:3121432/93‑1 (MQ=255)

                                                      ttACGATGTTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGc                                   <  1:5418620/93‑1 (MQ=255)

                                                       tACGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAATCCAGGAACGTAAAGCGTTTCACTGGATGGCg                                  <  1:6373462/93‑1 (MQ=255)

                                                        aCGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGa                                 <  2:1682754/93‑1 (MQ=255)

                                                        aCGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGa                                 <  1:4281789/93‑1 (MQ=255)

                                                         cGATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAg                                <  1:3623804/93‑1 (MQ=255)

                                                          gATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGc                               <  2:1351471/93‑1 (MQ=255)

                                                          gATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGc                               <  2:3687837/93‑1 (MQ=255)

                                                          gATGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGc                               <  1:7516643/93‑1 (MQ=255)

                                                           aTGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGcc                              >  2:4524222/1‑93 (MQ=255)

                                                            tGGTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCa                             <  2:1615686/93‑1 (MQ=255)

                                                             ggTACGGTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCAt                            <  1:7240360/93‑1 (MQ=255)

                                                                  ggTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGcc                       >  2:729894/1‑93 (MQ=255)

                                                                  ggTGTTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGcc                       >  1:1499335/1‑93 (MQ=255)

                                                                   gtgtTCCCACCTAACGTTGAAGGCGGTCGTTTTGCCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCt                      <  2:4332157/93‑1 (MQ=255)

                                                                    tgtTCCCCCCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTg                     <  1:8589851/93‑1 (MQ=255)

                                                                     gtTCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGa                    <  1:5836005/93‑1 (MQ=255)

                                                                      ttCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGAt                   <  1:2393273/93‑1 (MQ=255)

                                                                      ttCCCACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGACTGAt                   >  1:7872433/1‑93 (MQ=255)

                                                                       tCCCCCCTAATGTGGAAGGAGGTTGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATc                  <  1:2781416/93‑1 (MQ=255)

                                                                        cccACCTACCGTTGAAGGCGGTCGTTTTCTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCCCTGGATGGCGAGCCATAAGCCTGATCt                 <  2:3338482/93‑1 (MQ=255)

                                                                          cACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTgg               <  1:1005375/93‑1 (MQ=255)

                                                                          cACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTgg               <  2:3593901/93‑1 (MQ=255)

                                                                          cACCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTgg               >  2:763690/1‑93 (MQ=255)

                                                                           cccTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGt              <  1:1205068/92‑1 (MQ=255)

                                                                           aCCTAACGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCTTGATCTGGt              >  1:4634363/1‑93 (MQ=255)

                                                                             cTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtca            >  2:8095231/1‑93 (MQ=255)

                                                                               aaCGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGAGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtcatc          >  2:3896170/1‑93 (MQ=255)

                                                                               aaCGTTGAAGGCGGTCGTTTTGTCTCATGGCGGGAAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtcatc          <  2:7058095/93‑1 (MQ=255)

                                                                                aCGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtcatca         >  1:1350955/1‑93 (MQ=255)

                                                                                aCGTTGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGtcatca         <  1:2522009/93‑1 (MQ=255)

                                                                                   ttGAAGGCGGTCGTTTTTTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAAc      <  2:5456182/93‑1 (MQ=255)

                                                                                    tGAAGGCGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACt     >  1:3805694/1‑93 (MQ=255)

                                                                                      aaGGCGGTTTTTTTTTCTCATGGCTGGCAATCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCa   <  1:6658468/93‑1 (MQ=255)

                                                                                       aGGCGGTCGTTTTGTTTCAGGGCGGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCaa  <  1:5625194/93‑1 (MQ=255)

TCGCGGGTTAATTGTTCTAACCGACCGTTATCCTCAGGACCAAATTCCTGGCGCTTACGATGGTACGGTGTTCCCACCTAACGTTGAAGGTGGTCGTTTTGTCTCATGGCTGGCAAGCCAGGAACGTAAAGCGTTTCACTGGATGGCGAGCCATAAGCCTGATCTGGTCATCAAACTCAA  >  NC_000913/3132460‑3132639

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