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	93.3%	1195.9 / 53.4	430	A35A (GCA→GCG)	yjcH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (16/13); new base C (212/189); total (228/202)
Fisher's exact test for biased strand distribution p-value = 8.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 443 total aligned reads displayed.

AATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAAT  >  NC_000913/4282598‑4282780

aaTACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGCAGGCGATCAGTAAAATCACACCGATATAAACCGc                                                                                            >  1:7771973/1‑93 (MQ=255)

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

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

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

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

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

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

        gtgtGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCataa                                                                                    >  2:7621212/1‑93 (MQ=255)

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

           tgACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCAtaataa                                                                                 <  1:3766149/93‑1 (MQ=255)

            gACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATaataat                                                                                <  2:708293/93‑1 (MQ=255)

             aCGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATc                                                                               <  1:4570619/93‑1 (MQ=255)

             aCGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATc                                                                               <  1:4738693/93‑1 (MQ=255)

             aCGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCACCATAATAATc                                                                               >  1:1913727/1‑93 (MQ=255)

              cGCTTGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCg                                                                              <  1:2429363/93‑1 (MQ=255)

               gCTGGCGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGa                                                                             <  2:2248026/93‑1 (MQ=255)

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

                cTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGAc                                                                            <  2:598880/93‑1 (MQ=255)

                 tGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACa                                                                           <  1:352992/93‑1 (MQ=255)

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

                   gtgtTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAgg                                                                         >  1:10008839/1‑93 (MQ=255)

                    tgtTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGCCAGGa                                                                        >  2:2023646/1‑93 (MQ=255)

                    tgtTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAggg                                                                        >  1:7444657/1‑92 (MQ=255)

                      ttCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGAtg                                                                      >  1:6435910/1‑93 (MQ=255)

                       tCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGAtgg                                                                     <  2:6921560/93‑1 (MQ=255)

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

                         ggATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGAtggtg                                                                   >  1:4150499/1‑93 (MQ=255)

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

                          gATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGAtggtgg                                                                  >  1:7945874/1‑93 (MQ=255)

                          gATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGAtggtgg                                                                  >  1:6298008/1‑93 (MQ=255)

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

                            ttCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCa                                                                <  1:2332915/93‑1 (MQ=255)

                            ttCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCa                                                                <  1:8779151/93‑1 (MQ=255)

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

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

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

                                gcggcgTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCCGCATAATAATCGACAGGATGGTGGCAAAcc                                                            >  1:4993781/1‑93 (MQ=255)

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

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

                                       gCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGCGGCAAACCGTTGCCg                                                     >  1:2765887/1‑93 (MQ=255)

                                        cccAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGt                                                    <  1:4309041/93‑1 (MQ=255)

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

                                           agccagccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGtttt                                                 >  1:5539328/1‑93 (MQ=255)

                                            gccagccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGttttt                                                <  2:4622764/93‑1 (MQ=255)

                                              cagccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCg                                              >  1:5154713/1‑93 (MQ=255)

                                              cagccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCg                                              >  1:283940/1‑93 (MQ=255)

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

                                                gccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGAc                                            <  2:2884959/93‑1 (MQ=255)

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

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

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

                                                    ggCGCGAAGGCGAGCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAAc                                        <  1:2636133/93‑1 (MQ=255)

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

                                                       gcgAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTcc                                     >  2:6548934/1‑93 (MQ=255)

                                                       gcgAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTcc                                     >  1:7227793/1‑93 (MQ=255)

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

                                                         gAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCt                                   >  1:2127840/1‑93 (MQ=255)

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

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

                                                         gAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGGTTTTCGACTAACTCCCt                                   >  1:6549551/1‑93 (MQ=255)

                                                         gAAGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATCGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCt                                   >  2:1844441/1‑93 (MQ=255)

                                                          aaGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGTCTAACTCCCTg                                  >  2:1666515/1‑93 (MQ=255)

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

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

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

                                                           aGGCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGa                                 >  2:1084469/1‑93 (MQ=255)

                                                             gCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGaaa                               <  2:1941468/93‑1 (MQ=255)

                                                             gCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGaaa                               <  1:7022402/93‑1 (MQ=255)

                                                             gCGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGaaa                               <  1:322274/93‑1 (MQ=255)

                                                              cGATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTCACTCCCTGCAAt                              >  2:3936304/1‑93 (MQ=255)

                                                               gATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATg                             <  2:930231/93‑1 (MQ=255)

                                                               gATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATg                             >  2:5696563/1‑93 (MQ=255)

                                                               gATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATg                             >  1:68847/1‑93 (MQ=255)

                                                               gATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATg                             <  1:5858141/93‑1 (MQ=255)

                                                               gATCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATg                             >  2:9157794/1‑93 (MQ=255)

                                                                aTCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATgc                            >  2:9645210/1‑93 (MQ=255)

                                                                aTCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATgc                            <  1:5697422/93‑1 (MQ=255)

                                                                 tCAGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATgcg                           >  2:3051442/1‑93 (MQ=255)

                                                                   aGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCa                         <  2:3859718/93‑1 (MQ=255)

                                                                   aGTAAAATAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCa                         <  1:8590398/93‑1 (MQ=255)

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

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

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

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

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

                                                                        aaTAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGtc                    <  2:8845877/93‑1 (MQ=255)

                                                                        aaTAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGtc                    <  2:64093/93‑1 (MQ=255)

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

                                                                        aaTAAAGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGtc                    <  1:9461019/93‑1 (MQ=255)

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

                                                                           aaaGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGtcttc                 <  2:6694623/93‑1 (MQ=255)

                                                                             aGCCGATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGGCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTa               <  2:4273320/93‑1 (MQ=255)

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

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

                                                                                 gATATAAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCg           >  2:4183377/1‑93 (MQ=255)

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

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

                                                                                    atgaaCCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTg        <  1:1631053/90‑1 (MQ=255)

                                                                                    ataAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTAGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTg        >  1:4056226/1‑93 (MQ=255)

                                                                                     taAACCGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGCAATGCGCATTGTCTTCTATCCGCTGa       >  1:5368843/1‑93 (MQ=255)

AATACCGCGTGTGACGCTGGTGTTCGGATTCAGCGGCGTGCCCAGCCAGCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAAT  >  NC_000913/4282598‑4282780

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