BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,764,283	T→C	V107A (GTT→GCT)	yibA →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,764,283	0	T	C	84.1%	726.8 / 140.3	339	V107A (GTT→GCT)	yibA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (26/28); new base C (138/147); total (164/175)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

Only 100 of 341 total aligned reads displayed.

TACTCGGGCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAATCTCAAATTACTGCTT  >  NC_000913/3764193‑3764375

taCTCGGGCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTc                                                                                            >  2:3108954/1‑93 (MQ=255)

taCTCGGGCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTc                                                                                            >  1:3186409/1‑93 (MQ=255)

 aCTCGGGCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCg                                                                                           <  2:7331567/93‑1 (MQ=255)

  cTCGGGCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGa                                                                                          >  1:3821143/1‑93 (MQ=255)

   tCGGGCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAg                                                                                         >  2:554182/1‑93 (MQ=37)

    cGGGCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGc                                                                                        <  1:6492050/93‑1 (MQ=255)

     gggCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCt                                                                                       <  2:6485714/93‑1 (MQ=255)

     gggCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCt                                                                                       <  1:436640/93‑1 (MQ=255)

     gggCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCt                                                                                       >  2:5999654/1‑93 (MQ=255)

      ggCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTa                                                                                      <  2:1606830/93‑1 (MQ=255)

      ggCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTa                                                                                      <  2:5392462/93‑1 (MQ=255)

        cAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACg                                                                                    <  2:4339453/93‑1 (MQ=255)

        cAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACg                                                                                    >  2:4668274/1‑93 (MQ=37)

         aaaTAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACgg                                                                                   <  1:2694558/93‑1 (MQ=255)

          aaTAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGc                                                                                  >  1:340964/1‑93 (MQ=37)

          aaTAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGc                                                                                  >  2:366498/1‑93 (MQ=37)

           aTAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCa                                                                                 >  2:1996821/1‑93 (MQ=37)

             aaaaaTTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAAt                                                                               <  1:2212191/93‑1 (MQ=255)

             aaaaaTTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAAt                                                                               <  2:4091234/93‑1 (MQ=255)

                aaTTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGa                                                                            <  1:7338433/93‑1 (MQ=21)

                aaTTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGa                                                                            >  2:5137108/1‑93 (MQ=18)

                 aTTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAg                                                                           <  2:2203147/93‑1 (MQ=21)

                 aTTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAg                                                                           <  1:4084040/93‑1 (MQ=21)

                 aTTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAg                                                                           <  2:5895985/93‑1 (MQ=18)

                  tttGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGt                                                                          >  2:5748578/1‑93 (MQ=18)

                    tGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCa                                                                        >  1:7263067/1‑93 (MQ=33)

                    tGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCa                                                                        <  1:3550802/93‑1 (MQ=34)

                     gCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCaa                                                                       <  1:7414774/93‑1 (MQ=21)

                     gCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCaa                                                                       <  2:4495199/93‑1 (MQ=18)

                      cAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAAc                                                                      <  1:3744070/93‑1 (MQ=34)

                       aaaaaaTGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACg                                                                     >  2:3751343/1‑93 (MQ=18)

                       aaaaaaTGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACg                                                                     >  2:6587397/1‑93 (MQ=18)

                         aaaaTGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGc                                                                   >  2:5843621/1‑93 (MQ=21)

                          aaaTGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCt                                                                  >  2:119269/1‑93 (MQ=25)

                           aaTGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTc                                                                 >  1:6185477/1‑93 (MQ=25)

                            aTGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCa                                                                <  2:5223739/93‑1 (MQ=255)

                              gCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGc                                                              <  2:2562985/93‑1 (MQ=255)

                              gCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGc                                                              <  1:4885753/93‑1 (MQ=255)

                                gAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGa                                                            >  2:6865739/1‑93 (MQ=255)

                                 aaGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGAt                                                           >  2:3341720/1‑93 (MQ=255)

                                 aaGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGAt                                                           >  2:883348/1‑93 (MQ=255)

                                   gATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGc                                                         >  2:5832378/1‑93 (MQ=255)

                                     tAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCaa                                                       <  2:6230956/93‑1 (MQ=255)

                                      aaTGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAg                                                      >  2:4158734/1‑93 (MQ=255)

                                        tGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGaa                                                    >  1:5489980/1‑93 (MQ=255)

                                        tGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGaa                                                    <  2:4933827/93‑1 (MQ=255)

                                        tGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGaa                                                    >  1:2722889/1‑93 (MQ=255)

                                        tGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGaa                                                    <  2:2425326/93‑1 (MQ=255)

                                          tttttAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGaaaa                                                  <  2:751431/93‑1 (MQ=255)

                                           ttttAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGaaaaa                                                 <  1:637291/93‑1 (MQ=255)

                                           ttttAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGaaaaa                                                 >  1:6035834/1‑93 (MQ=255)

                                            tttAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAAc                                                <  1:6937214/93‑1 (MQ=255)

                                            tttAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAAc                                                >  2:1082605/1‑93 (MQ=255)

                                              tAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAAccc                                              >  1:6285933/1‑93 (MQ=255)

                                                atatTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCaa                                            <  1:4012068/93‑1 (MQ=255)

                                                 tatTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAAt                                           >  2:7086381/1‑93 (MQ=255)

                                                  atTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAAtt                                          >  2:4764854/1‑93 (MQ=255)

                                                  atTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAAtt                                          <  1:3699180/93‑1 (MQ=255)

                                                   ttttGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAAttt                                         >  2:6164380/1‑93 (MQ=255)

                                                    tttGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTa                                        <  1:1760573/93‑1 (MQ=255)

                                                     ttGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTAt                                       <  1:4468635/93‑1 (MQ=255)

                                                      tGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTAtt                                      >  1:4933999/1‑93 (MQ=255)

                                                      tGAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTAtt                                      <  2:6314156/93‑1 (MQ=255)

                                                       gAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTGAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTc                                     >  1:2907252/1‑93 (MQ=255)

                                                       gAACAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTc                                     <  1:2089532/93‑1 (MQ=255)

                                                        aaCAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCa                                    <  1:3944934/93‑1 (MQ=255)

                                                        aaCAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCa                                    >  2:2803773/1‑93 (MQ=255)

                                                        aaCAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCa                                    <  1:5891335/93‑1 (MQ=255)

                                                          cAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCAcc                                  <  2:2076751/93‑1 (MQ=255)

                                                          cAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCAcc                                  <  2:1837178/93‑1 (MQ=255)

                                                          cAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCAcc                                  <  2:1362620/93‑1 (MQ=255)

                                                          cAATATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCAcc                                  <  2:4550220/93‑1 (MQ=255)

                                                           aaTATGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCt                                 >  2:2881818/1‑93 (MQ=255)

                                                              atGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTaaa                              <  2:683067/93‑1 (MQ=255)

                                                              atGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTaaa                              <  2:4878565/93‑1 (MQ=255)

                                                               tGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTGCGGCAATCGAGTCAACAGCTCAGCGATGCAAGTAAAACCCAATTTATTCACCTaaaa                             >  2:4433111/1‑93 (MQ=255)

                                                               tGGCATTGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTaaaa                             <  1:2477750/93‑1 (MQ=255)

                                                                 gCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAAta                           <  1:1156966/93‑1 (MQ=255)

                                                                 gCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAAta                           >  1:1303019/1‑93 (MQ=255)

                                                                    ttGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAAtagta                        <  1:911310/93‑1 (MQ=255)

                                                                    ttGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAAtagta                        <  1:1920775/93‑1 (MQ=255)

                                                                     tGAATGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAAtagtag                       >  1:1047933/1‑93 (MQ=255)

                                                                       aaTGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGaa                     <  1:7142315/93‑1 (MQ=255)

                                                                        aTGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAAc                    <  1:525512/93‑1 (MQ=255)

                                                                         tGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACa                   >  1:4127293/1‑93 (MQ=255)

                                                                         tGATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACa                   >  1:4608466/1‑93 (MQ=255)

                                                                          gATAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACaa                  >  2:7170862/1‑93 (MQ=255)

                                                                           aTAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAAt                 <  1:3322254/93‑1 (MQ=255)

                                                                            tAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAAtc                >  2:2728933/1‑93 (MQ=255)

                                                                            tAAGAGCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAAtc                >  2:3445938/1‑93 (MQ=255)

                                                                              agagCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAAtctc              >  1:6055903/1‑93 (MQ=255)

                                                                               gagCGCTTGCCCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTCTTCACCTAAAATAGTAGAACAATCTCa             <  2:6425837/93‑1 (MQ=255)

                                                                                agCGCTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAATCTCaa            <  2:3540366/93‑1 (MQ=255)

                                                                                 gcgcTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAATCTCaaa           <  2:3182206/93‑1 (MQ=255)

                                                                                  cgcTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAATCTCAAAt          <  1:5605737/93‑1 (MQ=255)

                                                                                   gcTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAATCTCAAAtt         >  1:6859359/1‑93 (MQ=255)

                                                                                   gcTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAATCTCAAAtt         >  2:2700945/1‑93 (MQ=255)

                                                                                   gcTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAATCTCAAAtt         >  2:4107352/1‑93 (MQ=255)

                                                                                    cTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAATCTCAAATTa        <  2:6689186/93‑1 (MQ=255)

                                                                                    cTTGCGCTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAATCTCAAGTTa        >  2:1544966/1‑93 (MQ=255)

TACTCGGGCAAATAAAAATTTGCAAAAAATGCGAAGATAATGTTTTTAATATTTTGAACAATATGGCATTGAATGATAAGAGCGCTTGCGTTCGAGCTACGGCAATCGAGTCAACGGCTCAGCGATGCAAGAAAAACCCAATTTATTCACCTAAAATAGTAGAACAATCTCAAATTACTGCTT  >  NC_000913/3764193‑3764375

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