BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,382,443	A→G	T357A (ACG→GCG)	yjeM →	putative transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,382,443	0	A	G	98.4%	573.5 / ‑0.8	183	T357A (ACG→GCG)	yjeM	putative transport
Reads supporting (aligned to +/- strand): ref base A (0/3); new base G (105/74); total (105/78)
Fisher's exact test for biased strand distribution p-value = 7.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 191 total aligned reads displayed.

CGCGTATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGACGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGGGTTGGTT  >  NC_000913/4382354‑4382532

cgcgTATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGgcgc                                                                                        <  1:745987/93‑1 (MQ=255)

   gTATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGa                                                                                     >  2:3135886/1‑93 (MQ=255)

   gTATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGa                                                                                     >  2:3001128/1‑93 (MQ=255)

   gTATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGa                                                                                     <  1:571213/93‑1 (MQ=255)

    tATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGaa                                                                                    >  2:3665852/1‑93 (MQ=255)

    tATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGaa                                                                                    >  1:23261/1‑93 (MQ=255)

    tATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGaa                                                                                    <  1:3632310/93‑1 (MQ=255)

       tACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGc                                                                                 >  1:4621033/1‑93 (MQ=255)

        aCTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCa                                                                                >  1:4137489/1‑93 (MQ=255)

         cTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCAt                                                                               >  2:161956/1‑93 (MQ=255)

          tGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCAtt                                                                              >  2:4513318/1‑93 (MQ=255)

          tGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCAtt                                                                              >  2:2318353/1‑93 (MQ=255)

            gACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATtgt                                                                            <  1:4288112/93‑1 (MQ=255)

            gACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATtgt                                                                            <  2:1715513/93‑1 (MQ=255)

             gcTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATtgtg                                                                           <  1:1433812/92‑1 (MQ=255)

             gcTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATtgtg                                                                           <  1:4011928/92‑1 (MQ=255)

             aCTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATtgtg                                                                           >  2:567032/1‑93 (MQ=255)

               tttCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGc                                                                         >  2:2700181/1‑93 (MQ=255)

                      gTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAAcc                                                                  >  1:4198090/1‑93 (MQ=255)

                       ttCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCg                                                                 <  1:2006012/93‑1 (MQ=255)

                       ttCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCg                                                                 <  1:2006011/93‑1 (MQ=255)

                        tCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAAccca                                                                >  2:1170989/1‑91 (MQ=255)

                        tCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGa                                                                >  2:1136985/1‑93 (MQ=255)

                        tCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGa                                                                >  2:3319662/1‑93 (MQ=255)

                        tCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGa                                                                <  2:1575975/93‑1 (MQ=255)

                        tCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGACGCCGAAAGCATTGTGGCCGGAACCGa                                                                <  1:1570888/93‑1 (MQ=255)

                          cTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATg                                                              >  1:3851631/1‑93 (MQ=255)

                           tCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGa                                                             <  1:1318785/93‑1 (MQ=255)

                           tCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGa                                                             >  1:2551435/1‑93 (MQ=255)

                           tCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGa                                                             >  1:3962573/1‑93 (MQ=255)

                           ccGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGa                                                             <  1:2349837/92‑1 (MQ=255)

                            cGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGAc                                                            <  1:2719078/93‑1 (MQ=255)

                             gCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGAcg                                                           <  2:674919/93‑1 (MQ=255)

                             gCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGAcg                                                           >  2:3309529/1‑93 (MQ=255)

                              ccTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGAcgc                                                          >  2:326552/1‑93 (MQ=255)

                              ccTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGAcgc                                                          <  1:882325/93‑1 (MQ=255)

                              ccTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGAcgc                                                          >  2:1386703/1‑93 (MQ=255)

                               gtataCCGGCGGGAGCTTTACGCGTTGGGATTGACCGTTGATATCCGTCATCCACGGGGCGCTGAAAGCATTGTGGCCGGAACCGATGAcgcg                                                         >  2:1214907/2‑93 (MQ=255)

                               cTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGAcgcg                                                         >  1:3718454/1‑93 (MQ=255)

                                tataCCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACgcgc                                                        >  1:2251313/1‑93 (MQ=255)

                                tataCCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACgcgc                                                        >  1:1696984/1‑93 (MQ=255)

                                 ataCCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGcc                                                       <  2:1897107/93‑1 (MQ=255)

                                 ataCCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGcc                                                       <  2:1527838/93‑1 (MQ=255)

                                 ataCCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGcc                                                       <  2:2741112/93‑1 (MQ=255)

                                 ataCCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGcc                                                       >  2:3507943/1‑93 (MQ=255)

                                      ggTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAAt                                                  <  2:1273501/93‑1 (MQ=255)

                                       gTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATg                                                 >  2:1087376/1‑93 (MQ=255)

                                        tGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGc                                                >  1:3896714/1‑93 (MQ=255)

                                         gCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCg                                               >  2:4379693/1‑93 (MQ=255)

                                         gCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCg                                               >  1:523444/1‑93 (MQ=255)

                                           gTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGAt                                             >  1:1869300/1‑93 (MQ=255)

                                           gTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGAt                                             <  2:4204937/93‑1 (MQ=255)

                                            ttCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATg                                            <  1:4504463/93‑1 (MQ=255)

                                            ttCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATg                                            <  2:3763922/93‑1 (MQ=255)

                                             tCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATgg                                           <  1:3508243/93‑1 (MQ=255)

                                               tttACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATgggg                                         <  2:3870437/93‑1 (MQ=255)

                                               tttACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATgggg                                         >  2:322160/1‑93 (MQ=255)

                                               tttACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATgggg                                         >  2:977533/1‑93 (MQ=255)

                                                ttACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGa                                        <  1:1120995/93‑1 (MQ=255)

                                                 tACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGAt                                       <  1:3619346/93‑1 (MQ=255)

                                                  aCGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATg                                      >  1:764867/1‑93 (MQ=255)

                                                    gCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGcc                                    >  2:2536091/1‑93 (MQ=255)

                                                       ttGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTc                                 >  2:4246155/1‑93 (MQ=255)

                                                        tGCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCt                                >  1:3705283/1‑93 (MQ=255)

                                                         gCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTa                               <  2:716436/93‑1 (MQ=255)

                                                         gCTATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTa                               <  1:3079519/93‑1 (MQ=255)

                                                           tATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATc                             >  1:778337/1‑93 (MQ=255)

                                                           tATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATc                             <  2:393944/93‑1 (MQ=255)

                                                           tATTCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATc                             >  1:1497496/1‑93 (MQ=255)

                                                             ttCACCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATCGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGc                           >  2:2733254/1‑93 (MQ=255)

                                                                aCCGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCAt                        <  1:2843771/93‑1 (MQ=255)

                                                                 ccGTTGAAAGCCATCATCGAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCAtg                       >  2:4643967/1‑93 (MQ=255)

                                                                 ccGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCAtg                       >  2:1134262/1‑93 (MQ=255)

                                                                  cGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCAtgt                      <  1:1650369/93‑1 (MQ=255)

                                                                  cGTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCAtgg                      >  1:2509582/1‑92 (MQ=255)

                                                                   gTTGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCAtgtg                     <  2:4586432/93‑1 (MQ=255)

                                                                    ttGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTgg                    >  1:3443201/1‑93 (MQ=255)

                                                                    ttGAAAGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTgg                    >  2:3619271/1‑93 (MQ=255)

                                                                         aGCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCa               >  1:1080108/1‑93 (MQ=255)

                                                                          gCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAg              >  1:3795328/1‑93 (MQ=255)

                                                                          gCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAg              >  1:2113587/1‑93 (MQ=255)

                                                                          gCCATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAg              >  1:4667568/1‑93 (MQ=255)

                                                                           ccATCATCCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGt             <  2:752937/93‑1 (MQ=255)

                                                                             atcatcCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGc           <  2:1813740/93‑1 (MQ=255)

                                                                              tcatcCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCg          >  1:2776519/1‑93 (MQ=255)

                                                                              tcatcCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCg          >  1:4348771/1‑93 (MQ=255)

                                                                               catcCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCgg         <  2:2455973/93‑1 (MQ=255)

                                                                               catcCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCgg         <  1:4477579/93‑1 (MQ=255)

                                                                               catcCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCgg         <  1:1816266/93‑1 (MQ=255)

                                                                                atcCAGGGGGCGCCGAAAGCATTGTGGCCGGACCCGATGACGCGCCTGAATGCGATGGGGATCCCGTGTATCTCCATGCGGATCCAGTGCggg        >  2:4058686/1‑93 (MQ=255)

                                                                                atcCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGTGGATGCCGTCTATCGCCATGTGGATGCAGTGCggg        <  2:4676536/93‑1 (MQ=255)

                                                                                 tcCAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggt       <  2:2984687/93‑1 (MQ=255)

                                                                                  ccAGCGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggtt      <  1:2883476/93‑1 (MQ=255)

                                                                                   cAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttg     >  1:1492103/1‑93 (MQ=255)

                                                                                   cAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttg     >  1:2696581/1‑93 (MQ=255)

                                                                                   cAGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttg     >  1:2477753/1‑93 (MQ=255)

                                                                                    aGGGGGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttgg    >  1:3299673/1‑93 (MQ=255)

                                                                                     ggggGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttggt   >  1:2409788/1‑93 (MQ=255)

                                                                                     ggggGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttggt   <  1:2434248/93‑1 (MQ=255)

                                                                                     ggggGCGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGggttggt   >  2:1225178/1‑93 (MQ=255)

CGCGTATTACTGGACTTTCGATGTTCCTCGCCTATACCGGTGCGTTCTTTACGCTTTGCTATTCACCGTTGAAAGCCATCATCCAGGGGACGCCGAAAGCATTGTGGCCGGAACCGATGACGCGCCTGAATGCGATGGGGATGCCGTCTATCGCCATGTGGATGCAGTGCGGGTTGGTT  >  NC_000913/4382354‑4382532

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