BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,366,832	T→C	M387V (ATG→GTG)	yjeH ←	putative transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,366,832	0	T	C	98.3%	371.0 / ‑2.8	119	M387V (ATG→GTG)	yjeH	putative transport
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (65/52); total (66/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

Only 100 of 128 total aligned reads displayed.

TTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTGCCAGCATGCATAAC  >  NC_000913/4366742‑4366923

ttATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTg                                                                                           <  2:4575906/93‑1 (MQ=255)

ttATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTg                                                                                           >  1:99104/1‑93 (MQ=255)

   tGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCa                                                                                        >  1:2655239/1‑93 (MQ=255)

   tGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCa                                                                                        >  1:3338083/1‑93 (MQ=255)

   tGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCa                                                                                        >  1:2201640/1‑93 (MQ=255)

     ccATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGt                                                                                      <  2:3469589/93‑1 (MQ=255)

     ccATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAACGCATAGAGACTTTTCCAGCCGACCACTGCCTGt                                                                                      >  1:1495875/1‑93 (MQ=255)

       aTTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTaa                                                                                    >  2:2932656/1‑93 (MQ=255)

       aTTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTaa                                                                                    >  2:3518589/1‑93 (MQ=255)

        ttttCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAAc                                                                                   >  1:514348/1‑93 (MQ=255)

        ttttCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACGTTTCCAGCCGACCACTGCCAGTAAc                                                                                   >  2:1482846/1‑93 (MQ=255)

           tCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGa                                                                                <  2:3574245/93‑1 (MQ=255)

            ccGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCTACCACTGCCAGTAACAGaa                                                                               >  1:816247/1‑93 (MQ=255)

            ccGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGaa                                                                               <  2:1715546/93‑1 (MQ=255)

            ccGGCGTTTTTCGTTTGGGCGTCAACAGCCATAACCACGCCAGCATGATCAGCGCATAGAGCCTTTTCCAGCCGACCTCTGCCAGTAACAGaa                                                                               <  2:900489/93‑1 (MQ=255)

             cGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCCTAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAAc                                                                              >  2:2840730/1‑93 (MQ=255)

             cGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACGGAAc                                                                              <  2:289594/93‑1 (MQ=255)

             cGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAAc                                                                              >  2:4187271/1‑93 (MQ=255)

               gCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGc                                                                            >  1:2091848/1‑93 (MQ=255)

               gCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGc                                                                            >  2:2973150/1‑93 (MQ=255)

                cGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCa                                                                           <  2:3084904/93‑1 (MQ=255)

                 gTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCAt                                                                          >  2:1583981/1‑93 (MQ=255)

                 ctttttCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCAt                                                                          <  2:4042892/92‑1 (MQ=255)

                    tttCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATCGAGACTTTTCCAACCGACCACTGCCAGTAACAGAACGCATaca                                                                       >  1:2553413/1‑91 (MQ=255)

                    tttCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAAc                                                                       >  1:2553412/1‑93 (MQ=255)

                    tttCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAAc                                                                       <  2:449516/93‑1 (MQ=255)

                     ttCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACa                                                                      <  1:851897/93‑1 (MQ=255)

                       cGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGc                                                                    <  2:2403701/93‑1 (MQ=255)

                       cGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGc                                                                    <  1:4180822/93‑1 (MQ=255)

                       cGTTTGGGCAGCAACAGCCATACCCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGc                                                                    <  2:1267244/93‑1 (MQ=255)

                        gTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGcc                                                                   >  1:2678188/1‑93 (MQ=255)

                        gTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGcc                                                                   <  1:2327127/93‑1 (MQ=255)

                        gTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGcc                                                                   >  1:1572964/1‑93 (MQ=255)

                        gTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGcc                                                                   >  1:1179931/1‑93 (MQ=255)

                         ttttGGCAGAAACAGCCATAGCCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAAAGCATAACAGccc                                                                  <  2:2327117/93‑1 (MQ=255)

                           ttGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGc                                                                >  1:3043368/1‑93 (MQ=255)

                           ttGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGc                                                                <  1:3404023/93‑1 (MQ=255)

                            tGGCAGCAACAGCCATACCCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGcc                                                               <  2:1375024/93‑1 (MQ=255)

                            gggCAGCACCACCCATACCCCCGCCAGTATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGcc                                                               <  1:4249345/92‑1 (MQ=255)

                               cagcaACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAac                                                            <  2:2513439/93‑1 (MQ=255)

                                 gcaACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAacca                                                          >  1:328738/1‑93 (MQ=255)

                                  caACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAaccac                                                         >  1:251/1‑93 (MQ=255)

                                   aaCAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAaccacc                                                        <  1:1718999/93‑1 (MQ=255)

                                    aCAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCg                                                       >  2:3518893/1‑93 (MQ=255)

                                     cAGCCATACCCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGc                                                      <  1:4206652/93‑1 (MQ=255)

                                      aGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGcc                                                     <  2:93899/93‑1 (MQ=255)

                                        ccATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCag                                                   >  1:885003/1‑93 (MQ=255)

                                        ccATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCag                                                   >  2:4170352/1‑93 (MQ=255)

                                          aTAACCCCGCCCGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCCGCCACCCCCGACagta                                                 >  2:2360919/1‑93 (MQ=255)

                                          aTAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCagta                                                 >  1:3410161/1‑93 (MQ=255)

                                           tAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCagtag                                                >  2:2977882/1‑93 (MQ=255)

                                           tAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCagtag                                                <  2:3552954/93‑1 (MQ=255)

                                            acccccGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCagtagt                                               <  1:102898/91‑1 (MQ=255)

                                            acccccGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCagtagt                                               <  2:2618720/91‑1 (MQ=255)

                                             aCCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTc                                              >  2:4174734/1‑93 (MQ=255)

                                               cccGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGa                                            >  1:818979/1‑93 (MQ=255)

                                                ccGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGAt                                           >  2:1800341/1‑93 (MQ=255)

                                                ccGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGAt                                           >  2:1225485/1‑93 (MQ=255)

                                                   ccAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAAc                                        >  2:1907827/1‑93 (MQ=255)

                                                   ccAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAAc                                        >  1:2203339/1‑93 (MQ=255)

                                                       cATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTcc                                    <  2:3536383/93‑1 (MQ=255)

                                                       cATGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTcc                                    <  2:3536523/93‑1 (MQ=255)

                                                        aTGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCt                                   <  1:3124580/93‑1 (MQ=255)

                                                        aTGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCt                                   >  2:343794/1‑93 (MQ=255)

                                                        aTGATCAGCGCATAGAGACTTTCCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATATCGTCCt                                   <  2:3380348/93‑1 (MQ=255)

                                                         tGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAGCCACCGCCAGTAGTCGATAACGTCCtt                                  >  1:3279834/1‑93 (MQ=255)

                                                         tGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCtt                                  >  1:1442182/1‑93 (MQ=255)

                                                         tGATCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCtt                                  >  1:2508009/1‑93 (MQ=255)

                                                           aTCAGCGCATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGc                                >  1:2977999/1‑93 (MQ=255)

                                                                 gcATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCaataat                          >  1:3629154/1‑93 (MQ=255)

                                                                 gatgagagAATTTTCCAGCCCACCACTGCGAGTAACAGAACGCATAACAGCCCGCCAACCTCCGCCAGTAGTCGATAACGTCCTTGCaataat                          <  1:3642982/89‑1 (MQ=255)

                                                                  cATAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAAtt                         >  2:3931701/1‑93 (MQ=255)

                                                                   aTAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAAttt                        >  1:1735923/1‑93 (MQ=255)

                                                                    tAGAGACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGACAACCAACGCCAGTAGTCGATAACGTCTTTGCAATAATTTa                       >  2:3777575/1‑93 (MQ=255)

                                                                     agagACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTAc                      <  2:467752/93‑1 (MQ=255)

                                                                     agagACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTAc                      <  1:3649633/93‑1 (MQ=255)

                                                                     agagACTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTAc                      >  1:998978/1‑93 (MQ=255)

                                                                      gagaCTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACa                     >  1:1357009/1‑93 (MQ=255)

                                                                      gagaCTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACa                     <  1:2141200/93‑1 (MQ=255)

                                                                      gagaCTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACa                     <  2:3838765/93‑1 (MQ=255)

                                                                      aagaCTTTTCCAGCCGACCACTGCCAGTACCCGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACa                     <  1:1181046/92‑1 (MQ=255)

                                                                       agaCTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAg                    <  2:4658188/93‑1 (MQ=255)

                                                                         gcTTTTCCAGACGACCACTGCCAGGAACAGAACGCATAACAGCCCGCCAACCGCCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGcc                  <  2:4467671/92‑1 (MQ=255)

                                                                         aCTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGcc                  <  2:782179/93‑1 (MQ=255)

                                                                         aCTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGcc                  >  1:1388068/1‑93 (MQ=255)

                                                                          cTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCt                 <  2:3547527/93‑1 (MQ=255)

                                                                          cTTTTCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCt                 >  2:1828831/1‑93 (MQ=255)

                                                                           ttttCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTg                >  2:2167510/1‑93 (MQ=255)

                                                                           tgtgccAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTg                <  2:3318844/89‑1 (MQ=255)

                                                                            ttgccAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTGc               <  1:3086968/90‑1 (MQ=255)

                                                                             ttCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTGcc              <  2:778828/93‑1 (MQ=255)

                                                                              tCCAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTGCCa             <  1:3398021/93‑1 (MQ=255)

                                                                               ccAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTGCCac            >  2:2019352/1‑92 (MQ=255)

                                                                                cAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTGCCAgc           >  2:94518/1‑93 (MQ=255)

                                                                                cAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTGCCAgc           <  2:2116742/93‑1 (MQ=255)

                                                                                cAGCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTGCCAgc           >  2:1222330/1‑93 (MQ=255)

                                                                                  gCCGACCACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTGCCAgcat         <  2:41908/93‑1 (MQ=255)

                                                                                      aCCACTGCAAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTGCCAgcatgcat     <  1:3766860/93‑1 (MQ=255)

                                                                                       ccACTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTGCCAGCATGCATa    <  1:4337244/93‑1 (MQ=255)

                                                                                         aTTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTGCCAGCATGCATAAc  >  1:2218659/1‑93 (MQ=255)

TTATGCCATTTTCCGGCGTTTTTCGTTTTGGCAGCAACAGCCATAACCCCGCCAGCATGATCAGCGCATAGAGACTTTTCCAGCCGACCATTGCCAGTAACAGAACGCATAACAGCCCGCCAACCACCGCCAGTAGTCGATAACGTCCTTGCAATAATTTACAGCCTGCCAGCATGCATAAC  >  NC_000913/4366742‑4366923

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