BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	675,063	A→G	A236A (GCT→GCC)	ybeQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	675,063	0	A	G	97.7%	406.0 / ‑3.6	131	A236A (GCT→GCC)	ybeQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/0); new base G (76/52); total (79/52)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

CCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTG  >  NC_000913/674973‑675152

ccATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGc                                                                                         >  2:646365/1‑93 (MQ=255)

  aTTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCgg                                                                                       <  1:3763783/93‑1 (MQ=255)

    tCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCa                                                                                     >  2:1825509/1‑93 (MQ=255)

    tCAAAGGCAGCCTGATAATCTTTATCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCa                                                                                     >  1:3559018/1‑93 (MQ=255)

      aaaGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGc                                                                                   >  1:4161520/1‑93 (MQ=255)

       aaGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGcc                                                                                  >  1:1521067/1‑93 (MQ=255)

        aGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCt                                                                                 >  2:1129848/1‑93 (MQ=255)

        aGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCt                                                                                 >  1:3977271/1‑93 (MQ=255)

        aGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCt                                                                                 >  2:3853890/1‑93 (MQ=255)

        aGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGGGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCt                                                                                 >  2:3562029/1‑93 (MQ=255)

         ggCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCtt                                                                                >  1:4658644/1‑93 (MQ=255)

         ggCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCtt                                                                                >  2:3536612/1‑93 (MQ=255)

          gCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTc                                                                               <  2:1179829/93‑1 (MQ=255)

          gCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTc                                                                               <  1:3488843/93‑1 (MQ=255)

          gCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTc                                                                               >  1:3142135/1‑93 (MQ=255)

           cAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCa                                                                              <  1:2083441/93‑1 (MQ=255)

             gCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGa                                                                            >  1:1096582/1‑93 (MQ=255)

                 gATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATAcc                                                                        >  1:2798688/1‑93 (MQ=255)

                  aTAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCa                                                                       <  2:909897/93‑1 (MQ=255)

                  aTAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCa                                                                       <  1:2431937/93‑1 (MQ=255)

                   tAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAg                                                                      >  2:1338652/1‑93 (MQ=255)

                   tAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAg                                                                      >  2:4721300/1‑93 (MQ=255)

                     aTCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTa                                                                    >  2:287954/1‑93 (MQ=255)

                     aTCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTa                                                                    >  2:3743776/1‑93 (MQ=255)

                      tCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTAc                                                                   >  1:4372000/1‑93 (MQ=255)

                       cTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACa                                                                  <  1:4740283/93‑1 (MQ=255)

                        tttCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACAt                                                                 >  2:1170343/1‑93 (MQ=255)

                        tttCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACAt                                                                 <  2:2785468/93‑1 (MQ=255)

                        tttCTCAACGCCTTGTCCGTGTTAATACATATAACCGATGATGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACAt                                                                 <  1:4252597/93‑1 (MQ=255)

                         ttCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATc                                                                >  2:2734211/1‑93 (MQ=255)

                         ttCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATc                                                                >  1:819883/1‑93 (MQ=255)

                         ttCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATc                                                                <  2:165767/93‑1 (MQ=255)

                          tctcAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCg                                                               >  1:3330524/1‑93 (MQ=255)

                            tcAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgct                                                             >  2:3797304/1‑93 (MQ=255)

                            tcAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgct                                                             >  1:1096094/1‑93 (MQ=255)

                            tcAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgct                                                             >  2:679958/1‑93 (MQ=255)

                             cAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgctt                                                            >  1:3101631/1‑93 (MQ=255)

                              aaCGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgcttg                                                           >  2:4618485/1‑93 (MQ=255)

                              aaCGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCgcttg                                                           >  2:4531684/1‑93 (MQ=255)

                               aCGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgcttgc                                                          <  1:3225236/93‑1 (MQ=255)

                                cGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgcttgct                                                         <  2:102392/93‑1 (MQ=255)

                                cGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgcttgct                                                         <  1:1189229/93‑1 (MQ=255)

                                cGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCgcttgct                                                         <  1:4190403/93‑1 (MQ=255)

                                  ccTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTta                                                       <  1:3179096/93‑1 (MQ=255)

                                   cTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtat                                                      <  1:3762642/93‑1 (MQ=255)

                                    ttGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtata                                                     <  1:1152015/93‑1 (MQ=255)

                                    ttGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTtata                                                     >  2:750555/1‑93 (MQ=255)

                                     tGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAg                                                    >  1:2699607/1‑93 (MQ=255)

                                      gTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGt                                                   >  2:2195939/1‑93 (MQ=255)

                                       tCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTc                                                  <  1:3435934/93‑1 (MQ=255)

                                       tCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTc                                                  >  2:3845536/1‑93 (MQ=255)

                                        ccGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTct                                                 >  2:2743370/1‑92 (MQ=255)

                                         cGTGTTTATACATATACTCGATGTTGACTTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCt                                                <  1:290979/93‑1 (MQ=255)

                                         cGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCt                                                <  2:4693561/93‑1 (MQ=255)

                                         cGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCt                                                <  1:1688040/93‑1 (MQ=255)

                                          gtgtTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTg                                               >  2:1295406/1‑93 (MQ=255)

                                          gtgtTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTg                                               >  1:837929/1‑93 (MQ=255)

                                          gtgtTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTg                                               <  2:657292/93‑1 (MQ=255)

                                           tgtTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTgg                                              <  2:304187/93‑1 (MQ=255)

                                           tgtTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTgg                                              >  1:866638/1‑93 (MQ=255)

                                           tgtTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTgg                                              >  2:1853073/1‑93 (MQ=255)

                                              ttATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTcc                                           <  1:2188216/93‑1 (MQ=255)

                                                 taCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAcc                                        <  2:1988011/93‑1 (MQ=255)

                                                 taCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAcc                                        >  1:982025/1‑93 (MQ=255)

                                                 taCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAcc                                        >  1:1113209/1‑93 (MQ=255)

                                                 taCATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAcc                                        >  1:4110315/1‑93 (MQ=255)

                                                   cATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCAcccc                                      <  1:1609922/93‑1 (MQ=255)

                                                    atatAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCt                                     <  1:2041044/93‑1 (MQ=255)

                                                    atatAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTGGCTTATAGTCCTGGTCCACCCCt                                     >  2:2591149/1‑93 (MQ=255)

                                                      ataACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTc                                   >  2:3641494/1‑93 (MQ=255)

                                                      ataACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTc                                   >  2:3641493/1‑93 (MQ=255)

                                                       taACCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCg                                  >  2:4316907/1‑93 (MQ=255)

                                                        aaCCGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGc                                 <  2:3960932/93‑1 (MQ=255)

                                                          ccGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCa                               <  1:2317010/93‑1 (MQ=255)

                                                           cGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAg                              >  1:1273270/1‑93 (MQ=255)

                                                           cGATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAg                              <  2:189674/93‑1 (MQ=255)

                                                            gATGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGc                             >  2:153690/1‑93 (MQ=255)

                                                             aTGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCg                            >  1:3302509/1‑93 (MQ=255)

                                                              tGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATTCCAGTACATCGCTTGCTTATAGTCCTGGTCCTCCCCTTCGCCAGCGt                           <  1:1871354/93‑1 (MQ=255)

                                                              tGTTGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGt                           >  2:1254936/1‑93 (MQ=255)

                                                                ttGACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTta                         >  1:3005564/1‑93 (MQ=255)

                                                                ttGACGTAAGCGCCGACGCTCCCCTGGGCGGAAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTta                         >  1:1693154/1‑93 (MQ=255)

                                                                  gACGTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTtata                       >  2:2839520/1‑93 (MQ=255)

                                                                    agTAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAAt                     <  1:2790105/92‑1 (MQ=255)

                                                                      tAAGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATcc                   >  2:428193/1‑93 (MQ=255)

                                                                        aGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCa                 <  2:415542/93‑1 (MQ=255)

                                                                        aGCGCCGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCa                 <  2:3591531/93‑1 (MQ=255)

                                                                             cGACGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCaa            <  2:4056380/93‑1 (MQ=255)

                                                                              gACGCTCCCCTGGGCGGCAGCCTTCATATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAAt           >  1:3507862/1‑93 (MQ=255)

                                                                              caCGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAAt           <  1:2808051/92‑1 (MQ=255)

                                                                               aCGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAAtt          >  2:322891/1‑93 (MQ=255)

                                                                               aCGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAAtt          >  2:3776560/1‑93 (MQ=255)

                                                                                cGCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAAttt         >  1:1781732/1‑93 (MQ=255)

                                                                                 gCTCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTg        >  2:2871130/1‑93 (MQ=255)

                                                                                   tCCCCTGGGCGGCAGCCTTCAGATACCAGTACATCGGTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGaa      >  2:4460318/1‑93 (MQ=255)

                                                                                      ccTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACt   >  2:2063238/1‑93 (MQ=255)

                                                                                      ccTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACt   >  1:4411867/1‑93 (MQ=255)

                                                                                       cTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTg  >  1:542501/1‑93 (MQ=255)

                                                                                       cTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTg  >  2:4655206/1‑93 (MQ=255)

                                                                                       cTGGGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTg  <  2:2545186/93‑1 (MQ=255)

CCATTCAAAGGCAGCCTGATAATCTTTCTCAACGCCTTGTCCGTGTTTATACATATAACCGATGTTGACGTAAGCGCCGACGCTCCCCTGAGCGGCAGCCTTCAGATACCAGTACATCGCTTGCTTATAGTCCTGGTCCACCCCTTCGCCAGCGTTATAATCCCATGCAATTTGAAACTG  >  NC_000913/674973‑675152

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: `‑`