BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,597,072	T→C	D73G (GAC→GGC)	b3460 ←	LivJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,597,072	0	T	C	99.6%	870.8 / ‑5.0	274	D73G (GAC→GGC)	b3460	LivJ
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (154/119); total (155/119)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 276 total aligned reads displayed.

GAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGTCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGG  >  NC_000913/3596982‑3597161

   gaaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCatat                                                                                      <  2:712330/93‑1 (MQ=255)

    aaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAtt                                                                                     <  1:1215933/93‑1 (MQ=255)

    aaCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAtt                                                                                     <  1:319205/93‑1 (MQ=255)

     aCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAttt                                                                                    >  1:2721492/1‑93 (MQ=255)

     aCAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACAGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATAttt                                                                                    >  2:1566796/1‑93 (MQ=255)

       agagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTa                                                                                  >  2:750153/1‑93 (MQ=255)

       agagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTa                                                                                  >  2:2407740/1‑93 (MQ=255)

       agagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTa                                                                                  <  1:3440741/93‑1 (MQ=255)

        gagGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTAc                                                                                 >  2:1655134/1‑93 (MQ=255)

          ggTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACgt                                                                               >  2:2923817/1‑92 (MQ=255)

          ggTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGa                                                                               <  1:1970698/93‑1 (MQ=255)

             gACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGAttt                                                                            <  2:1087309/93‑1 (MQ=255)

             gACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGAttt                                                                            <  1:558552/93‑1 (MQ=255)

             gACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGACACAAGCATCGCCATATTTTACGAttt                                                                            >  2:1540370/1‑93 (MQ=255)

               ccAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGc                                                                          >  1:2794111/1‑93 (MQ=255)

                cAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCa                                                                         <  1:1007301/93‑1 (MQ=255)

                cAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCa                                                                         >  2:581957/1‑93 (MQ=255)

                cAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCa                                                                         >  2:4163927/1‑93 (MQ=255)

                 aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCCACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAg                                                                        >  2:3783313/1‑93 (MQ=255)

                 aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAg                                                                        >  1:4593903/1‑93 (MQ=255)

                 aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAg                                                                        >  2:1753443/1‑93 (MQ=255)

                 aaTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAg                                                                        <  2:4042342/93‑1 (MQ=255)

                  aTCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGt                                                                       >  2:3726348/1‑93 (MQ=255)

                   tCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGtt                                                                      <  1:2521289/93‑1 (MQ=255)

                   tCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGtt                                                                      >  2:4096270/1‑93 (MQ=255)

                    cacaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGttt                                                                     <  2:11451/93‑1 (MQ=255)

                     acaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttg                                                                    >  2:882618/1‑93 (MQ=255)

                     acaTATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttg                                                                    >  1:2020022/1‑93 (MQ=255)

                       atatTTCATTCCGTCGTTAACGACATTGAACACATCCGCAAACGCCTATTTAGAGACACAGGCATCGCCATAAAATACGATTTGCAGTatgtt                                                                  <  2:4686637/93‑1 (MQ=255)

                       atatTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgtt                                                                  <  1:2231607/93‑1 (MQ=255)

                        tatTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgttg                                                                 >  2:1296103/1‑93 (MQ=255)

                        tatTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTttgttg                                                                 >  1:3823973/1‑93 (MQ=255)

                          tttAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGcc                                                               >  2:3140971/1‑93 (MQ=255)

                            tAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCtt                                                             >  1:2342785/1‑93 (MQ=255)

                             aaTGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCttt                                                            >  1:4431846/1‑93 (MQ=255)

                             aaTGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCttt                                                            >  1:2161068/1‑93 (MQ=255)

                              aTGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTa                                                           >  1:3512308/1‑93 (MQ=255)

                               tGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTaa                                                          <  2:2173112/93‑1 (MQ=255)

                               tGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTaa                                                          <  1:2408132/93‑1 (MQ=255)

                               tGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTaa                                                          <  1:2906799/93‑1 (MQ=255)

                               tGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTaa                                                          <  2:1583879/93‑1 (MQ=255)

                                gCCGTCGTTAACGACTTTGTTCGCCACCGCAACCTCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAAt                                                         <  1:627102/93‑1 (MQ=255)

                                gCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAAt                                                         >  2:1189244/1‑93 (MQ=255)

                                gCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAAt                                                         <  1:4543835/93‑1 (MQ=255)

                                 ccGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATg                                                        <  2:3386975/93‑1 (MQ=255)

                                  cgtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgc                                                       >  1:932973/1‑93 (MQ=255)

                                   gtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgcc                                                      <  1:4272746/93‑1 (MQ=255)

                                   gtcgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgcc                                                      >  2:1207194/1‑93 (MQ=255)

                                     cgtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgccgc                                                    >  2:2005657/1‑93 (MQ=255)

                                      gtTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATgccgcc                                                   <  2:4272601/93‑1 (MQ=255)

                                            taCTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGc                                             <  1:3752259/92‑1 (MQ=255)

                                            gACTTTGTTCGCCACCGCACCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGc                                             >  2:4568043/1‑93 (MQ=255)

                                             aCTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCg                                            >  2:2504082/1‑93 (MQ=255)

                                             aCTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCg                                            >  2:3035588/1‑93 (MQ=255)

                                              cTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGt                                           >  2:911444/1‑93 (MQ=255)

                                              cTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGt                                           >  2:2205368/1‑93 (MQ=255)

                                              cTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGt                                           <  1:720843/93‑1 (MQ=255)

                                               tttGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGtt                                          >  1:1602486/1‑93 (MQ=255)

                                               tttGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGtt                                          <  2:4080301/93‑1 (MQ=255)

                                                 tGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGa                                        <  1:4106942/93‑1 (MQ=255)

                                                  gTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGat                                       <  1:1473484/93‑1 (MQ=255)

                                                   ttCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGata                                      <  1:4390908/93‑1 (MQ=255)

                                                   ttCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGata                                      <  2:4455757/93‑1 (MQ=255)

                                                    tCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGatat                                     <  2:3828950/93‑1 (MQ=255)

                                                    tCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGatat                                     <  2:3250331/93‑1 (MQ=255)

                                                     cGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATc                                    <  2:2662110/93‑1 (MQ=255)

                                                      gCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATcc                                   <  1:3770361/93‑1 (MQ=255)

                                                         aCCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCa                                >  1:2945024/1‑93 (MQ=255)

                                                         aCCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCa                                >  1:3103438/1‑93 (MQ=255)

                                                         aCCGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCa                                >  2:3082065/1‑93 (MQ=255)

                                                          ccGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCaa                               <  2:3670995/93‑1 (MQ=255)

                                                          ccGCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCaa                               >  1:4130790/1‑93 (MQ=255)

                                                           cGCAACCGCCTGTTTCGTGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAAc                              <  2:2795660/93‑1 (MQ=255)

                                                            gCAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAAcc                             >  1:1381248/1‑93 (MQ=255)

                                                             cAACCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCg                            >  1:1303747/1‑93 (MQ=255)

                                                               aCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGcc                          <  1:3393881/93‑1 (MQ=255)

                                                               aCCGCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGcc                          >  2:4500668/1‑93 (MQ=255)

                                                                  gCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGc                       <  1:1022972/93‑1 (MQ=255)

                                                                  gCCTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGc                       >  1:3254411/1‑93 (MQ=255)

                                                                    cTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctc                     >  1:1725681/1‑93 (MQ=255)

                                                                    cTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctc                     >  1:1678250/1‑93 (MQ=255)

                                                                    cTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctc                     >  2:816551/1‑93 (MQ=255)

                                                                    cTGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctc                     <  1:1555163/93‑1 (MQ=255)

                                                                    cTGTTTCGGGACACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGACGCCTTTAGCGTTGATATCAGAAACCGCCTGctc                     >  2:1420652/1‑93 (MQ=255)

                                                                     tGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctct                    >  1:2837485/1‑93 (MQ=255)

                                                                     tGTTTCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGctct                    >  2:3285560/1‑93 (MQ=255)

                                                                       tttCGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTgc                  >  2:1493240/1‑93 (MQ=255)

                                                                          cGGGTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGcc               >  1:1141021/1‑93 (MQ=255)

                                                                             gTCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGt            <  2:524351/93‑1 (MQ=255)

                                                                              tCACAGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTa           >  1:369805/1‑93 (MQ=255)

                                                                                acaGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTaaa         >  2:1889217/1‑93 (MQ=255)

                                                                                 caGGCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAAc        <  2:2830142/93‑1 (MQ=255)

                                                                                   ggCATCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTc      >  1:3937274/1‑93 (MQ=255)

                                                                                      aTCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTg   <  1:431031/93‑1 (MQ=255)

                                                                                      aTCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTg   >  1:1299446/1‑93 (MQ=255)

                                                                                      aTCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTg   <  2:3516995/93‑1 (MQ=255)

                                                                                       tCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACTGCCTGCTCTGCGACGGTAAACTCCTgg  >  1:1452135/1‑93 (MQ=255)

                                                                                       tCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTgg  >  2:4365724/1‑93 (MQ=255)

                                                                                       tCGCCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTgg  >  1:1660319/1‑93 (MQ=255)

                                                                                       tCGCCATATTTTACGATTTGCAGTTTGTTGCCGTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTgg  >  2:698004/1‑93 (MQ=255)

GAAGAACAGAGGTGACCAATCACATATTTAATGCCGTCGTTAACGACTTTGTTCGCCACCGCAACCGCCTGTTTCGGGTCACAGGCATCGTCATATTTTACGATTTGCAGTTTGTTGCCTTTAATGCCGCCTTTAGCGTTGATATCCGCAACCGCCTGCTCTGCGCCGGTAAACTCCTGG  >  NC_000913/3596982‑3597161

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