BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,334,367	G→A	S31S (AGC→AGT)	yrbA ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,334,367	0	G	A	83.1%	928.4 / 133.4	360	S31S (AGC→AGT)	yrbA	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (33/28); new base A (149/150); total (182/178)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

Only 100 of 362 total aligned reads displayed.

ATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATGGGATCTTCAATCAT  >  NC_000913/3334275‑3334459

atatatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTgg                                                                                              <  1:1210696/93‑1 (MQ=255)

  atatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACt                                                                                            >  2:543860/1‑93 (MQ=255)

     ttCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGcc                                                                                         <  2:5299452/93‑1 (MQ=255)

      tCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCa                                                                                        >  2:1162972/1‑93 (MQ=255)

       ccATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCAt                                                                                       >  1:2979535/1‑93 (MQ=255)

        cATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATc                                                                                      >  2:6675815/1‑93 (MQ=255)

          tCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGc                                                                                    >  2:1754878/1‑93 (MQ=255)

          tCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGc                                                                                    <  1:2172111/93‑1 (MQ=255)

            aGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCg                                                                                  >  2:6697227/1‑93 (MQ=255)

            aGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCAACGCCg                                                                                  >  2:4836508/1‑93 (MQ=255)

             gCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCgg                                                                                 >  1:4089737/1‑93 (MQ=255)

                gACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCGCGGCAATAACCTGAAAGTGACTGCCATCGCCGGaaa                                                                              <  1:5825918/93‑1 (MQ=255)

                gACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGaaa                                                                              <  1:1080264/93‑1 (MQ=255)

                  ccATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACg                                                                            >  2:1482032/1‑93 (MQ=255)

                   cATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGt                                                                           >  1:1015485/1‑93 (MQ=255)

                   cATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGt                                                                           >  1:6262080/1‑93 (MQ=255)

                      aGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGa                                                                        <  2:3946213/93‑1 (MQ=255)

                       gACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGAc                                                                       >  1:3926005/1‑93 (MQ=255)

                       gACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGAc                                                                       <  1:3217261/93‑1 (MQ=255)

                        aCCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACt                                                                      <  2:2302409/93‑1 (MQ=255)

                         ccGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACtt                                                                     >  2:2693541/1‑93 (MQ=255)

                         ccGTCTGCTGTTTTTTAACCCGACCCATGCCGTCAACCAACCCACCCCCGGCAAAACCCTGAAAGGGACTGCCCTCGCCGGAAACGTGGACtt                                                                     <  1:3628970/93‑1 (MQ=255)

                          cGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCTGAAACGTGGACTTc                                                                    <  2:5182297/93‑1 (MQ=255)

                          cGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTc                                                                    <  1:4231260/93‑1 (MQ=255)

                           gTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTcc                                                                   >  1:1810636/1‑93 (MQ=255)

                           gTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTcc                                                                   <  1:5240381/93‑1 (MQ=255)

                            tCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCt                                                                  >  1:3078593/1‑93 (MQ=255)

                             ctgctgTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTg                                                                 <  1:6337347/93‑1 (MQ=255)

                              tgctgTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTgg                                                                >  2:7047307/1‑93 (MQ=255)

                                ctgTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTggag                                                              >  1:3315614/1‑93 (MQ=255)

                                ctgTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTggag                                                              <  2:444925/93‑1 (MQ=255)

                                ctgTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTggag                                                              >  1:7165067/1‑93 (MQ=255)

                                 tgTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTggagg                                                             <  2:4037255/93‑1 (MQ=255)

                                 tgTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTggagg                                                             <  2:6300263/93‑1 (MQ=255)

                                 tgTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTggagg                                                             >  2:7379652/1‑93 (MQ=255)

                                 tgTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTggagg                                                             <  2:5046891/93‑1 (MQ=255)

                                  gTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTggaggg                                                            >  2:5663737/1‑93 (MQ=255)

                                   ttttttAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTggaggga                                                           <  2:4956636/93‑1 (MQ=255)

                                    tttttAACCCGACTCATGCCGTCAAACAACTCACCCCCGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTggagggag                                                          >  1:2518107/1‑93 (MQ=255)

                                     ttttAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGgaga                                                         <  2:652713/93‑1 (MQ=255)

                                       ttAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGgagaga                                                       <  1:2767904/93‑1 (MQ=255)

                                        tAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGagagag                                                      <  1:4367136/93‑1 (MQ=255)

                                         aaCCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGc                                                     >  2:155189/1‑93 (MQ=255)

                                          aCCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCg                                                    <  2:3647196/93‑1 (MQ=255)

                                          aCCCGACTCATGCCGTCAAACAACTCACCCACAGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCg                                                    >  2:1048731/1‑93 (MQ=255)

                                           cccGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGt                                                   <  2:4712052/93‑1 (MQ=255)

                                           cccGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGt                                                   <  2:6051613/93‑1 (MQ=255)

                                           cccGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGt                                                   >  1:2841612/1‑93 (MQ=255)

                                           cccGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGt                                                   >  2:7073831/1‑93 (MQ=255)

                                           cccGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGt                                                   <  2:2590127/93‑1 (MQ=255)

                                            ccGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGtt                                                  >  1:7059648/1‑93 (MQ=255)

                                             cGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtc                                                 <  1:2650560/93‑1 (MQ=255)

                                             cGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtc                                                 <  2:3412484/93‑1 (MQ=255)

                                              gACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtca                                                >  2:1842054/1‑93 (MQ=255)

                                                cTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtcatc                                              <  2:794208/93‑1 (MQ=255)

                                                cTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtcatc                                              >  2:7155711/1‑93 (MQ=255)

                                                cTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtcatc                                              <  1:2319452/93‑1 (MQ=255)

                                                 tCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtcatca                                             >  1:5453718/1‑93 (MQ=255)

                                                 tCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtcatca                                             >  2:2891224/1‑93 (MQ=255)

                                                  cATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCaa                                            <  2:5361607/93‑1 (MQ=255)

                                                  cATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCaa                                            <  2:5986527/93‑1 (MQ=255)

                                                  cATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCaa                                            <  2:1332532/93‑1 (MQ=255)

                                                  cATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCaa                                            <  2:2095415/93‑1 (MQ=255)

                                                      ccGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACg                                        <  2:6603365/93‑1 (MQ=255)

                                                       cGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGc                                       >  1:138962/1‑93 (MQ=255)

                                                       cGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGc                                       <  1:1274704/93‑1 (MQ=255)

                                                       cGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGc                                       <  2:3457093/93‑1 (MQ=255)

                                                       cGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGc                                       <  1:987282/93‑1 (MQ=255)

                                                       cGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGc                                       <  1:3701022/93‑1 (MQ=255)

                                                        gTCAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGct                                      >  2:4582564/1‑93 (MQ=255)

                                                          cAAACAACTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGctct                                    <  1:198411/93‑1 (MQ=255)

                                                             acaacTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGaa                                 <  2:5812082/93‑1 (MQ=255)

                                                             acaacTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGaa                                 >  1:6903906/1‑93 (MQ=255)

                                                              caacTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAAt                                <  1:4585932/93‑1 (MQ=255)

                                                                 cTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTc                             <  1:6047505/93‑1 (MQ=255)

                                                                 cTCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTc                             >  2:2509508/1‑93 (MQ=255)

                                                                  tCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCa                            <  1:4878929/93‑1 (MQ=255)

                                                                  tCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTTGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCa                            <  2:296573/93‑1 (MQ=255)

                                                                  tCACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCa                            >  2:1780509/1‑93 (MQ=255)

                                                                   cACCCACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCat                           <  2:1136829/93‑1 (MQ=255)

                                                                     cccACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCatta                         >  1:475800/1‑93 (MQ=255)

                                                                     cccACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCatta                         <  1:563992/93‑1 (MQ=255)

                                                                      ccACGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCattat                        <  2:2249923/93‑1 (MQ=255)

                                                                        aCGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTAttt                      >  1:1641127/1‑93 (MQ=255)

                                                                        aCGGCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTAttt                      >  1:1547262/1‑93 (MQ=255)

                                                                           gCAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTcc                   <  2:2006697/93‑1 (MQ=255)

                                                                            cAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCa                  <  2:4398191/93‑1 (MQ=255)

                                                                            cAATAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCCACACGCTCTGAATTTCATTATTTTCCa                  <  2:6417576/93‑1 (MQ=255)

                                                                              aTAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATg                <  1:6696669/93‑1 (MQ=255)

                                                                              aTAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATg                >  1:4918813/1‑93 (MQ=255)

                                                                              aTAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATg                <  2:2390031/93‑1 (MQ=255)

                                                                               tAACCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATgg               <  2:3282060/93‑1 (MQ=255)

                                                                                aaCCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATggg              >  1:7529195/1‑93 (MQ=255)

                                                                                aaCCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATggg              <  1:4471930/93‑1 (MQ=255)

                                                                                aaCCTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATggg              <  2:7472272/93‑1 (MQ=255)

                                                                                   cTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATGGGATc           <  1:4839157/93‑1 (MQ=255)

                                                                                   cTGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATGGGATc           <  1:644653/93‑1 (MQ=255)

                                                                                    tGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATGGGATCt          >  2:3841220/1‑93 (MQ=255)

                                                                                    tGAAAGTGACTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATGGGATCt          <  2:3144854/93‑1 (MQ=255)

                                                                                          tGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATGGGATCTTCAatc    >  1:6518897/1‑93 (MQ=255)

ATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGAATTTCATTATTTTCCATGGGATCTTCAATCAT  >  NC_000913/3334275‑3334459

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