BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,817,156	T→C	S550G (AGT→GGT)	yicF ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,817,156	0	T	C	99.5%	1233.1 / ‑3.7	394	S550G (AGT→GGT)	yicF	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (209/182); total (210/184)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

Only 100 of 413 total aligned reads displayed.

TACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACT  >  NC_000913/3817066‑3817245

 aCGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGcc                                                                                        <  1:7218628/93‑1 (MQ=255)

 aCGCCTTAAGCTGGGATTACCGATCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGcc                                                                                        >  2:4981270/1‑93 (MQ=255)

  cGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGccc                                                                                       <  1:2739052/93‑1 (MQ=255)

   gCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCa                                                                                      >  1:154293/1‑93 (MQ=255)

   gCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCa                                                                                      >  1:5845896/1‑93 (MQ=255)

   gCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCa                                                                                      >  2:2801944/1‑93 (MQ=255)

     cTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAGCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGc                                                                                    <  1:7520271/93‑1 (MQ=255)

     cTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAgt                                                                                    >  2:2693846/1‑92 (MQ=255)

     cTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGc                                                                                    <  2:899972/93‑1 (MQ=255)

     cTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGc                                                                                    >  2:1767175/1‑93 (MQ=255)

      ttAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGct                                                                                   >  2:773011/1‑93 (MQ=255)

      ttAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGct                                                                                   >  2:2495744/1‑93 (MQ=255)

         aGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttct                                                                                <  2:5316300/93‑1 (MQ=255)

          gCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttctt                                                                               <  1:3535911/93‑1 (MQ=255)

          gCTGGGATTACCGGTCTTGTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttctt                                                                               <  2:5340738/93‑1 (MQ=255)

          agtGGGATTACCGGTCTTTTTCTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGcttctt                                                                               <  1:2588953/91‑1 (MQ=255)

            tGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGa                                                                             <  1:7115083/93‑1 (MQ=255)

             gggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAt                                                                            <  1:6645484/93‑1 (MQ=255)

             gggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAt                                                                            <  2:5925034/93‑1 (MQ=255)

             gggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGAt                                                                            >  1:2400273/1‑93 (MQ=255)

              ggATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGTCAGCCAGCCAACCGCCCAGCTTCTTGAtt                                                                           >  1:4477385/1‑93 (MQ=255)

                 ttACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgc                                                                        <  2:3701295/93‑1 (MQ=255)

                 ttACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgc                                                                        >  2:5039339/1‑93 (MQ=255)

                   aCCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgcgc                                                                      >  2:4195188/1‑93 (MQ=255)

                   aCCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTgcgc                                                                      >  2:3518726/1‑93 (MQ=255)

                    ccGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCa                                                                     <  1:5053374/93‑1 (MQ=255)

                     cGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAt                                                                    <  1:1955041/93‑1 (MQ=255)

                     cGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAt                                                                    >  2:7228131/1‑93 (MQ=255)

                      ggTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAtt                                                                   >  1:6944924/1‑93 (MQ=255)

                      ggTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAtt                                                                   <  1:4523910/93‑1 (MQ=255)

                       gTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAttt                                                                  >  1:6057640/1‑93 (MQ=255)

                        tCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGGCCCGCTTCTTGATTTGCGCCtttt                                                                 >  2:1871807/1‑93 (MQ=255)

                        tCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCAtttt                                                                 >  2:3165680/1‑93 (MQ=255)

                         cTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttc                                                                >  1:3896332/1‑93 (MQ=255)

                          ttttACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttcc                                                               >  2:6606801/1‑93 (MQ=255)

                          ttttACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttcc                                                               <  2:6725040/93‑1 (MQ=255)

                           tttACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttcct                                                              >  2:1503665/1‑93 (MQ=255)

                             tactacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttccttc                                                            >  2:6515470/1‑93 (MQ=255)

                             tactacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTttccttc                                                            <  1:4137665/93‑1 (MQ=255)

                               ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                          >  2:3343951/1‑93 (MQ=255)

                               ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                          <  1:6755776/93‑1 (MQ=255)

                               ctacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                          >  1:4694834/1‑93 (MQ=255)

                               ctacAGCCACGAACGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCa                                                          <  2:5663034/93‑1 (MQ=255)

                                tacAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCCACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCAt                                                         >  1:1239815/1‑93 (MQ=255)

                                  cAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTc                                                       <  2:4496300/93‑1 (MQ=255)

                                   aGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCa                                                      <  1:1022816/93‑1 (MQ=255)

                                      cACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTTAATa                                                   >  2:2245144/1‑93 (MQ=255)

                                       aCGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                                  >  2:3045582/1‑93 (MQ=255)

                                       aCGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                                  >  1:309568/1‑93 (MQ=255)

                                       aCGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATaa                                                  <  1:3758399/93‑1 (MQ=255)

                                       aCGAAGGACACTAAGGTTCAAAACCCGTGATCTGCTGGGCAGCCAGCCAACCGTCCAGCCGCCGGACGCGCGCATTTTCCTGCCAGTCAATaa                                                  >  2:7045026/1‑93 (MQ=255)

                                         gAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAAcc                                                <  1:822611/93‑1 (MQ=255)

                                         gAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAAcc                                                <  2:4858757/93‑1 (MQ=255)

                                           aGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTg                                              >  1:6301638/1‑93 (MQ=255)

                                            ggACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGt                                             <  1:4413681/93‑1 (MQ=255)

                                             gACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTc                                            <  2:4014576/93‑1 (MQ=255)

                                              acacTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCt                                           <  2:277103/93‑1 (MQ=255)

                                               cacTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTg                                          <  1:2543476/93‑1 (MQ=255)

                                               cacTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTg                                          <  2:3566787/93‑1 (MQ=255)

                                                acTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTgg                                         <  1:7322145/93‑1 (MQ=255)

                                                acTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTGCCTGCCATTCAATAACCTGTCTgg                                         <  2:3100172/93‑1 (MQ=255)

                                                 cTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgc                                        >  2:2059605/1‑93 (MQ=255)

                                                 cTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgc                                        >  2:5747276/1‑93 (MQ=255)

                                                 cTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgc                                        >  2:7377186/1‑93 (MQ=255)

                                                  tAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGgcg                                       >  1:4004791/1‑93 (MQ=255)

                                                     ggTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGt                                    <  2:5119740/93‑1 (MQ=255)

                                                          aaaaCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGa                               >  2:7900674/1‑93 (MQ=255)

                                                          aaaaCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGa                               >  2:1312831/1‑93 (MQ=255)

                                                            aaCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATc                             <  2:5242973/93‑1 (MQ=255)

                                                            aaCCTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATc                             <  1:337230/93‑1 (MQ=255)

                                                             aCCTGTGATCTGCTGTGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGTGTCCCTATct                            >  2:5525749/1‑92 (MQ=255)

                                                              ccTGTGATCTGCTGGGCAGCCAGCCAACCGCCCTGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCa                           <  1:7619512/93‑1 (MQ=255)

                                                              ccTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCa                           <  1:2246514/93‑1 (MQ=255)

                                                              ccTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCa                           >  1:2955794/1‑93 (MQ=255)

                                                               cTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAg                          >  1:5027474/1‑93 (MQ=255)

                                                               cTGTGATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAg                          >  2:6851741/1‑93 (MQ=255)

                                                                  tgATCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTcc                       >  1:7840085/1‑93 (MQ=255)

                                                                    aTCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTcccc                     >  1:6439109/1‑93 (MQ=255)

                                                                     tCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCg                    >  2:4318222/1‑93 (MQ=255)

                                                                     tCTGCTGGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCg                    <  1:205952/93‑1 (MQ=255)

                                                                      ctgctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCgg                   <  1:1752600/93‑1 (MQ=255)

                                                                       tgctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGc                  <  1:4748577/93‑1 (MQ=255)

                                                                       tgctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGc                  >  2:461183/1‑93 (MQ=255)

                                                                       tgctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGc                  >  1:4247946/1‑93 (MQ=255)

                                                                        gctgGGCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCa                 >  2:4156901/1‑93 (MQ=255)

                                                                           gcgCACCCAGCCAACCGCCCAGCTTCTTGATTTGCCCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgct              <  1:2858210/91‑1 (MQ=255)

                                                                             gCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgctgc            <  2:1692592/93‑1 (MQ=255)

                                                                             gCAGCCAGCCAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAgctgc            <  1:6515756/93‑1 (MQ=255)

                                                                               agccagccAACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGctgctg          >  1:5555787/1‑93 (MQ=255)

                                                                                gccagccaCCCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc         <  2:2940334/93‑1 (MQ=255)

                                                                                gccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc         >  2:2185321/1‑93 (MQ=255)

                                                                                gccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCtgctgc         >  1:1318561/1‑93 (MQ=255)

                                                                                 ccagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGcc        >  1:638391/1‑93 (MQ=255)

                                                                                  cagccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCa       <  2:6026122/93‑1 (MQ=255)

                                                                                     ccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGaa    >  1:7404318/1‑93 (MQ=255)

                                                                                     ccaACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGaa    <  1:5282553/93‑1 (MQ=255)

                                                                                      caACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAAc   <  2:7509035/93‑1 (MQ=255)

                                                                                      caACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAAc   <  1:1877715/93‑1 (MQ=255)

                                                                                      caACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAAc   <  1:5655996/93‑1 (MQ=255)

                                                                                      caACCGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAAc   <  2:637378/93‑1 (MQ=255)

TACGCCTTAAGCTGGGATTACCGGTCTTTTACTACAGCCACGAAGGACACTAAGGTTCAAAACCTGTGATCTGCTGGGCAGCCAGCCAACTGCCCAGCTTCTTGATTTGCGCATTTTCCTTCCATTCAATAACCTGTCTGGCGCGTCCCGATCCAGTCCCCGGCAGCTGCTGCCAGAACT  >  NC_000913/3817066‑3817245

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