BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,282,688	T→C	A35A (GCA→GCG)	yjcH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,282,688	0	T	C	99.5%	702.8 / NA	218	A35A (GCA→GCG)	yjcH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (92/125); minor base A (1/0); total (93/125)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 222 total aligned reads displayed.

AATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATA  >  NC_000913/4282598‑4282777

aaTACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGc                                                                                         <  2:4740910/93‑1 (MQ=255)

 aTACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGcc                                                                                        >  2:5969553/1‑93 (MQ=255)

 aTACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGcc                                                                                        <  1:1884708/93‑1 (MQ=255)

 aTACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGcc                                                                                        <  2:616084/93‑1 (MQ=255)

  tACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCa                                                                                       <  1:4987812/93‑1 (MQ=255)

  tACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCa                                                                                       >  2:2500969/1‑93 (MQ=255)

   aCCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAg                                                                                      >  2:631304/1‑93 (MQ=255)

    ccGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGc                                                                                     >  1:224003/1‑93 (MQ=255)

     cgcgTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCa                                                                                    <  1:4500583/93‑1 (MQ=255)

      gcgTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCat                                                                                   >  2:4405360/1‑93 (MQ=255)

       cgTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCata                                                                                  >  1:4467181/1‑93 (MQ=255)

       cgTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCata                                                                                  >  1:2058799/1‑93 (MQ=255)

         tgtgACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCataat                                                                                <  1:5511175/93‑1 (MQ=255)

         tgtgACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCataat                                                                                <  1:6241132/93‑1 (MQ=255)

         ggtgCCGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCataat                                                                                <  1:3347873/92‑1 (MQ=255)

          gtgACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCataata                                                                               <  2:1478589/93‑1 (MQ=255)

          gtgACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCataata                                                                               >  1:272411/1‑93 (MQ=255)

            gACGCTGGTGTTCGGATTCAGCGGCGTGGCCAGCCCGCCGGGCGCGAAGGCGGTCAGTAAAATAAAAACGATATAAACCGCCAGCATaataat                                                                             >  2:3234842/1‑93 (MQ=255)

             aCGCTGGTGTTCGGATTCAGCGGCGTGCCGAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATc                                                                            <  1:4410123/93‑1 (MQ=255)

               gCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGa                                                                          >  2:903147/1‑93 (MQ=255)

                cTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGAc                                                                         >  2:860107/1‑93 (MQ=255)

                  ggTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAg                                                                       <  1:1907308/93‑1 (MQ=255)

                      ttCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGAtg                                                                   <  2:1304571/93‑1 (MQ=255)

                        cGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGAtggt                                                                 <  2:529796/93‑1 (MQ=255)

                        cGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGAtggt                                                                 <  2:5620423/93‑1 (MQ=255)

                         ggATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGAtggtg                                                                <  1:1783944/93‑1 (MQ=255)

                         ggATCCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGAtggtg                                                                <  2:380738/93‑1 (MQ=255)

                          gATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGAtggtgg                                                               <  1:2490113/93‑1 (MQ=255)

                           aTTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGc                                                              <  2:5377387/93‑1 (MQ=255)

                             tCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCaa                                                            >  1:227930/1‑93 (MQ=255)

                             tCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCaa                                                            >  2:2873649/1‑93 (MQ=255)

                              cAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCaaa                                                           >  1:706445/1‑93 (MQ=255)

                               aGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAAc                                                          <  2:5825474/93‑1 (MQ=255)

                                 cggcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCg                                                        <  1:3548882/93‑1 (MQ=255)

                                  ggcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGt                                                       <  2:3245201/93‑1 (MQ=255)

                                  ggcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGt                                                       <  1:2740631/93‑1 (MQ=255)

                                  ggcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGt                                                       <  1:5170649/93‑1 (MQ=255)

                                  ggcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGt                                                       <  1:5302125/93‑1 (MQ=255)

                                  ggcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGt                                                       <  1:2781021/93‑1 (MQ=255)

                                   gcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGtt                                                      <  1:3010040/93‑1 (MQ=255)

                                   gcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGtt                                                      <  1:1077903/93‑1 (MQ=255)

                                    cgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTg                                                     <  1:381121/93‑1 (MQ=255)

                                     gTGCCCAGCCAGCTGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAGCCGTTGc                                                    <  1:3605899/93‑1 (MQ=255)

                                     gTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGc                                                    <  1:2397809/93‑1 (MQ=255)

                                     gTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGc                                                    >  2:4558003/1‑93 (MQ=255)

                                     gTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGc                                                    >  1:2915612/1‑93 (MQ=255)

                                      tGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGcc                                                   >  1:3168396/1‑93 (MQ=255)

                                         ccagccagCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGtt                                                >  2:3680199/1‑93 (MQ=255)

                                          cagccagcCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGttt                                               >  1:4366934/1‑93 (MQ=255)

                                             ccagccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTc                                            <  2:2989738/93‑1 (MQ=255)

                                               agccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGa                                          <  1:5845803/93‑1 (MQ=255)

                                               agccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGa                                          <  2:1634701/93‑1 (MQ=255)

                                                gccGGGCGCGAAGGCGAGCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGAc                                         <  2:1971838/93‑1 (MQ=255)

                                                   gggCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTaa                                      >  1:4009045/1‑93 (MQ=255)

                                                   gggCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTaa                                      <  2:5940609/93‑1 (MQ=255)

                                                   gggCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTaa                                      >  2:3532578/1‑93 (MQ=255)

                                                    ggCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAAc                                     <  1:128520/93‑1 (MQ=255)

                                                      cgcgAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTc                                   <  1:3215600/93‑1 (MQ=255)

                                                        cgAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTccc                                 <  1:3819505/93‑1 (MQ=255)

                                                         gAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCt                                >  2:1835277/1‑93 (MQ=255)

                                                         gAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCt                                >  1:4282751/1‑93 (MQ=255)

                                                          aaGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTg                               >  2:4169020/1‑93 (MQ=255)

                                                            ggCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGaa                             >  2:6069982/1‑93 (MQ=255)

                                                               gATCAGTAAAATAAAACCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATg                          <  1:866503/93‑1 (MQ=255)

                                                                  cAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATgcgc                       <  1:3772748/93‑1 (MQ=255)

                                                                   aGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCa                      >  1:4906187/1‑93 (MQ=255)

                                                                    gTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCAt                     >  2:5772777/1‑93 (MQ=255)

                                                                    gTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCAt                     >  2:135245/1‑93 (MQ=255)

                                                                     tAAAATAAAGCCGATATAACCCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCAtt                    <  2:3412416/93‑1 (MQ=255)

                                                                     tAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCAtt                    <  2:4230087/93‑1 (MQ=255)

                                                                      aaaaTAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTg                   <  2:4989958/93‑1 (MQ=255)

                                                                      aaaaTAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTg                   >  2:2319906/1‑93 (MQ=255)

                                                                      aaaaTAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTg                   >  1:681840/1‑93 (MQ=255)

                                                                      aaaaTAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTg                   <  2:5426179/93‑1 (MQ=255)

                                                                       aaaTAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGt                  <  2:2071887/93‑1 (MQ=255)

                                                                       aaaTAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGt                  >  1:3473617/1‑93 (MQ=255)

                                                                        aaTAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGtc                 >  2:2807/1‑93 (MQ=255)

                                                                        aaTAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGtc                 >  2:1686073/1‑93 (MQ=255)

                                                                         aTAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGtct                >  1:6011173/1‑93 (MQ=255)

                                                                         aTAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGtct                <  1:4794659/93‑1 (MQ=255)

                                                                          tAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGtctt               >  1:1471379/1‑93 (MQ=255)

                                                                          tAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGtctt               >  1:3318932/1‑93 (MQ=255)

                                                                          tAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGtctt               >  2:4500079/1‑93 (MQ=255)

                                                                            aaGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGtcttct             <  2:2169650/93‑1 (MQ=255)

                                                                            aaGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGtcttct             >  1:5734077/1‑93 (MQ=255)

                                                                              gCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTAt           >  1:4587805/1‑93 (MQ=255)

                                                                               ccGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATc          <  1:2112773/93‑1 (MQ=255)

                                                                                cGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATcc         >  1:912045/1‑93 (MQ=255)

                                                                                cGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATcc         <  1:3869291/93‑1 (MQ=255)

                                                                                  atatAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATcctc       >  1:614572/1‑91 (MQ=255)

                                                                                  atatAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGc       >  1:1393109/1‑93 (MQ=255)

                                                                                  atatAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGc       <  2:3489749/93‑1 (MQ=255)

                                                                                  atatAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGc       <  1:996555/93‑1 (MQ=255)

                                                                                   tataAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCt      >  1:5597459/1‑93 (MQ=255)

                                                                                   tataAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCt      >  1:1209438/1‑93 (MQ=255)

                                                                                   tataAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTACCTGAAATGCGCATTGTCTTCTATCCGCt      >  1:5097000/1‑93 (MQ=255)

                                                                                     taAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGa    <  1:2678974/93‑1 (MQ=255)

                                                                                     aaaaCCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGa    <  1:2690855/92‑1 (MQ=255)

                                                                                      aaaCCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGAt   <  1:5994768/93‑1 (MQ=255)

                                                                                      aaaCCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGAt   >  2:2686755/1‑93 (MQ=255)

AATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATA  >  NC_000913/4282598‑4282777

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