BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,370,182	A→G	N527S (AAC→AGC)	mopA →	GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,370,182	0	A	G	99.5%	1205.8 / ‑5.2	371	N527S (AAC→AGC)	mopA	GroEL
Reads supporting (aligned to +/- strand): ref base A (1/1); new base G (164/205); total (165/206)
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.94e-01

Only 100 of 393 total aligned reads displayed.

CCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAAT  >  NC_000913/4370092‑4370271

 caaAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGa                                                                                        <  1:2316838/93‑1 (MQ=255)

  aaaGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGAt                                                                                       <  2:8782917/93‑1 (MQ=255)

   aaGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATg                                                                                      <  1:6206905/93‑1 (MQ=255)

   aaGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGAAGATCACCACCGAATGCCTGGTTTTCGGCCTGCCGAAAAGAGATg                                                                                      >  2:1708716/1‑93 (MQ=255)

    aGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGc                                                                                     >  1:118323/1‑93 (MQ=255)

    aGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGc                                                                                     >  1:2459580/1‑93 (MQ=255)

    aGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGc                                                                                     <  1:9523878/93‑1 (MQ=255)

      tAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGCCCTGCCGACAAGCGATGCAg                                                                                   >  1:7845062/1‑93 (MQ=255)

      tAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAg                                                                                   <  1:4693420/93‑1 (MQ=255)

       aaCTCGTTCTGCTCTGCAGTATGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGc                                                                                  >  1:6335403/1‑93 (MQ=255)

           cGTTCTGCTCTGCAGTACGCAGCTTCGGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGAc                                                                              <  1:5342036/93‑1 (MQ=255)

             ttCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACtt                                                                            >  2:3726763/1‑93 (MQ=255)

             ttCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACtt                                                                            <  2:9188438/93‑1 (MQ=255)

             ttCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACtt                                                                            >  1:3911587/1‑93 (MQ=255)

              tctgctctgcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTa                                                                           <  1:508728/93‑1 (MQ=255)

              tctgctctgcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTa                                                                           <  1:7249658/93‑1 (MQ=255)

               ctgctctgcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAg                                                                          <  2:3859992/93‑1 (MQ=255)

                tgctctgcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAgg                                                                         <  1:5635229/93‑1 (MQ=255)

                 gctctgcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGgc                                                                        <  2:1417674/93‑1 (MQ=255)

                  ctctgcAGTACGCAGCTTCTGTTGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGgcg                                                                       >  1:7115959/1‑93 (MQ=255)

                  ctctgcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGGAAAGCGATGCAGCTGAATTAggtg                                                                       >  1:9200649/1‑91 (MQ=255)

                  ctctgcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGgcg                                                                       >  2:1583429/1‑93 (MQ=255)

                   tctgcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGgcgc                                                                      >  2:9394903/1‑93 (MQ=255)

                      gcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCgctgc                                                                   >  2:1146411/1‑93 (MQ=255)

                      gcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCgctgc                                                                   <  2:3294537/93‑1 (MQ=255)

                      gcAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCgctgc                                                                   <  1:7754906/93‑1 (MQ=255)

                          tACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTggg                                                               >  2:1279215/1‑92 (MQ=255)

                          tACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGc                                                               >  2:7861934/1‑93 (MQ=255)

                          tACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGc                                                               <  1:4096591/93‑1 (MQ=255)

                          tACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGc                                                               <  1:6852747/93‑1 (MQ=255)

                            cGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCgg                                                             <  1:3739339/93‑1 (MQ=255)

                            cGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCgg                                                             <  2:3473348/93‑1 (MQ=255)

                             gCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGt                                                            >  1:1412350/1‑93 (MQ=255)

                              cAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTa                                                           >  2:727444/1‑93 (MQ=255)

                              cAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTa                                                           <  2:7627164/93‑1 (MQ=255)

                              cAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTa                                                           <  2:280415/93‑1 (MQ=255)

                              cAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTa                                                           <  2:5349087/93‑1 (MQ=255)

                                gCTTCTGTGGCTGGCCTGATGATCGCCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATg                                                         <  1:6479647/93‑1 (MQ=255)

                                gCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATg                                                         >  2:8554494/1‑93 (MQ=255)

                                gCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATg                                                         <  2:4588175/93‑1 (MQ=255)

                                gCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATg                                                         >  2:7952201/1‑93 (MQ=255)

                                gCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGAGCTGCTGGCGGGATg                                                         >  2:5945047/1‑93 (MQ=255)

                                 cTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATgg                                                        >  1:8139794/1‑93 (MQ=255)

                                 cTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATgg                                                        <  1:9069302/93‑1 (MQ=255)

                                  ttCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATggg                                                       >  2:6396935/1‑93 (MQ=255)

                                  ttCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATggg                                                       <  1:9387215/93‑1 (MQ=255)

                                      gtgGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGggcggc                                                   <  1:742125/93‑1 (MQ=255)

                                      gtgGCTGGCCTGATGATCACCACCGAATGCATGGGTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGggcggc                                                   >  2:1088496/1‑93 (MQ=255)

                                       tggctgtcCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCa                                                  <  2:5075854/93‑1 (MQ=255)

                                       tggctggcCTGATGATCCCCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCa                                                  <  2:5860727/93‑1 (MQ=255)

                                       tggctggcCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCa                                                  <  1:3011507/93‑1 (MQ=255)

                                       tggctgacCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCa                                                  <  1:10023220/93‑1 (MQ=255)

                                        ggctggcCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCAt                                                 >  2:2702394/1‑93 (MQ=255)

                                        ggctggcCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCAt                                                 <  2:2052180/93‑1 (MQ=255)

                                        ggctggcCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCAt                                                 <  1:5918915/93‑1 (MQ=255)

                                          ctggcCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATgg                                               <  1:263018/93‑1 (MQ=255)

                                           tggcCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATggg                                              <  2:3068185/93‑1 (MQ=255)

                                           tggcCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATggg                                              <  2:6910625/93‑1 (MQ=255)

                                            ggcCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATggtt                                             >  2:6596707/1‑91 (MQ=255)

                                            ggcCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGt                                             <  1:6053106/93‑1 (MQ=255)

                                             gcCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTg                                            >  2:6168923/1‑93 (MQ=255)

                                             gcCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTg                                            >  2:76193/1‑93 (MQ=255)

                                                tgatgaTCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCa                                         >  1:325119/1‑93 (MQ=255)

                                                tgatgaTCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCa                                         <  2:673402/93‑1 (MQ=255)

                                                tgatgaTCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCa                                         >  1:1585295/1‑93 (MQ=255)

                                                tgatgaTCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCa                                         <  2:9748409/93‑1 (MQ=255)

                                                 gatgatCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCAt                                        >  2:6032509/1‑93 (MQ=255)

                                                  atgatCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATg                                       <  1:4234718/93‑1 (MQ=255)

                                                   tgatCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTATGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATgg                                      >  2:2251263/1‑93 (MQ=255)

                                                   tgatCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATgg                                      >  2:4168343/1‑93 (MQ=255)

                                                    gatCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATggg                                     >  2:6241424/1‑93 (MQ=255)

                                                     atCACCACCGAATGCATGGTTACCGACCTGCCGAGAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGggc                                    <  2:2531033/93‑1 (MQ=255)

                                                     atCACCACCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGggc                                    <  1:7887977/93‑1 (MQ=255)

                                                       caccacCGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGggcgg                                  <  1:8793111/93‑1 (MQ=255)

                                                        accaccGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGggcggc                                 >  1:4163087/1‑93 (MQ=255)

                                                          caccGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCat                               >  1:2782103/1‑93 (MQ=255)

                                                          caccGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCat                               >  1:8829234/1‑93 (MQ=255)

                                                           accGAATGCATGGTTACCGACCTGCCGTAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGTCTGGGAGGCatg                              >  2:4851603/1‑93 (MQ=255)

                                                            ccGAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCatga                             <  1:4294659/93‑1 (MQ=255)

                                                              gAATGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCatgatg                           <  1:5039010/93‑1 (MQ=255)

                                                                aTGCATGGTTACCGGCCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTa                         <  1:2940684/93‑1 (MQ=255)

                                                                aTGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTa                         >  2:2343478/1‑93 (MQ=255)

                                                                 tGCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGGGGCATGATGTaa                        >  1:3531122/1‑93 (MQ=255)

                                                                  gCATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAAt                       <  2:4081765/93‑1 (MQ=255)

                                                                   cATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAAtt                      >  2:9864616/1‑93 (MQ=255)

                                                                   cATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAAtt                      <  1:3338698/93‑1 (MQ=255)

                                                                   cATGGTTACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAAtt                      <  1:1848745/93‑1 (MQ=255)

                                                                        ttACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGTGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCt                 <  1:9966187/93‑1 (MQ=255)

                                                                        gtACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCt                 <  1:8271492/92‑1 (MQ=255)

                                                                         tACCGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTg                <  1:6593554/93‑1 (MQ=255)

                                                                           ccGACCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCa              <  2:6680706/93‑1 (MQ=255)

                                                                             gACCTGCCGAAAAGCGATGCAGCTGACGTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCTTGGGCGGCATGATGTAATTGCCCTGCAcc            <  1:8589471/93‑1 (MQ=255)

                                                                              aCCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCt           >  2:7485656/1‑93 (MQ=255)

                                                                              aCCTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCt           >  1:4381654/1‑93 (MQ=255)

                                                                               ccTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTc          >  2:1840131/1‑93 (MQ=255)

                                                                               ccTGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTc          <  2:8672444/93‑1 (MQ=255)

                                                                                 tGCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGc        >  1:9742674/1‑93 (MQ=255)

                                                                                  gCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCa       >  1:963446/1‑93 (MQ=255)

                                                                                  gCCGAAAAGCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCa       >  1:5800412/1‑93 (MQ=255)

                                                                                      aaaagCGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGaaa   <  1:986856/93‑1 (MQ=255)

CCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACCTCGCAGAAAT  >  NC_000913/4370092‑4370271

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