BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	260,537	G→C	5.9%	S50S (TCG→TCC)	proB →	gamma‑glutamate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	260,537	0	G	C	5.9%	330.4 / 6.9	97	S50S (TCG→TCC)	proB	gamma‑glutamate kinase
Reads supporting (aligned to +/- strand): ref base G (44/46); new base C (2/4); total (46/51)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 101 total aligned reads displayed.

GAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCGGG  >  NC_000913/260390‑260684

gAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCg                                                                                                                                                      <  1:571695/148‑1 (MQ=255)

  gTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCggg                                                                                                                                                    <  1:178983/148‑1 (MQ=255)

      cAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGa                                                                                                                                              <  3:391171/149‑1 (MQ=255)

       aGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGgc‑gc                                                                                                                                                >  4:830383/1‑146 (MQ=255)

       aGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGgc‑gc                                                                                                                                                <  3:830383/146‑1 (MQ=255)

       aGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCAGGC‑GCGAt                                                                                                                                             <  4:259976/149‑1 (MQ=255)

        gCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGAt                                                                                                                                             >  1:22324/1‑148 (MQ=255)

        gCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGAt                                                                                                                                             <  1:749577/148‑1 (MQ=255)

        gCCAGACGCTGGTGGTAAAACTCGGCACCAGGGTGCTAACAGGCGGAGCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCCCGCCGCTACATGCCGCCGGGCATCGGATTGTTAATGTGACGTCCGgc‑ggcga                                                                                                                                             >  2:342772/1‑144 (MQ=255)

         ccAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATc                                                                                                                                            >  2:888255/1‑148 (MQ=255)

          cAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCggg                                                                                                                                                    >  1:181521/1‑140 (MQ=255)

          cAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCggg                                                                                                                                                    <  2:181521/140‑1 (MQ=255)

          cAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcgc                                                                                                                                          <  3:167235/149‑1 (MQ=255)

           aGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcgc                                                                                                                                          >  1:892042/1‑148 (MQ=255)

             aCGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcgccg                                                                                                                                        <  2:516439/148‑1 (MQ=255)

                  ggtggtAAGACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGa                                                                                                                                   <  2:786362/148‑1 (MQ=255)

                    tggtAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACg                                                                                                                                 <  2:205121/148‑1 (MQ=255)

                    tgggAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACg                                                                                                                                 >  1:309507/1‑148 (MQ=255)

                     ggtAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGt                                                                                                                                <  2:246034/148‑1 (MQ=255)

                     ggtAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCAGATTGTTATTGTGCCGTCCGGCGGCGATCG‑CGCCCGACGtc                                                                                                                               <  4:131035/149‑2 (MQ=255)

                      gtAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTg                                                                                                                               <  2:828097/148‑1 (MQ=255)

                      gtAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGa                                                                                                                              <  3:565008/149‑1 (MQ=255)

                       tAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATCG‑CGCCCGACGtca                                                                                                                              <  1:342772/148‑3 (MQ=255)

                         aaaCTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGc                                                                                                                            >  1:488741/1‑148 (MQ=255)

                         aaaCTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCa                                                                                                                           >  3:156774/1‑149 (MQ=255)

                             tCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTg                                                                                                                       <  3:531929/149‑1 (MQ=255)

                              cGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGTACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCtt                                                                                                                       >  2:237715/1‑147 (MQ=255)

                              cGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTg                                                                                                                       >  1:393601/1‑148 (MQ=255)

                                  aCCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTa                                                                                                                  >  3:367525/1‑149 (MQ=255)

                                   ccAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGGAcc                                                                                                                         <  1:776145/141‑1 (MQ=255)

                                   ccAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCAcc                                                                                                                         >  2:776145/1‑141 (MQ=255)

                                   ccAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATcgcgcccga                                                                                                                                    <  1:867679/131‑7 (MQ=255)

                                   ccAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCCGGCGGCGATcgcgcccga                                                                                                                                    >  2:867679/1‑125 (MQ=255)

                                    cAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTAcc                                                                                                                >  4:242417/1‑149 (MQ=255)

                                       tgtgCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCg                                                                                                              >  2:441225/1‑148 (MQ=255)

                                         tgCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGaa                                                                                                           >  3:193802/1‑149 (MQ=255)

                                         tgCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGa                                                                                                            <  2:409563/148‑1 (MQ=255)

                                          gCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCt                                                                                                                        <  3:18635/135‑1 (MQ=255)

                                          gCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCt                                                                                                                        >  4:18635/1‑135 (MQ=255)

                                            tAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACt                                                                                                         >  1:797826/1‑148 (MQ=255)

                                             aaCAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTg                                                                                                        >  2:896864/1‑148 (MQ=255)

                                                     gATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGAcc                                                                                               >  3:904999/1‑149 (MQ=255)

                                                        cgcgCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGAGACCAt                                                                                             <  2:488741/148‑1 (MQ=255)

                                                            ccGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGcc                                                                                         <  2:821002/148‑1 (MQ=255)

                                                            ccGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCt                                                                                        <  4:142608/149‑1 (MQ=255)

                                                             cGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCt                                                                                        <  2:505660/148‑1 (MQ=255)

                                                             cGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTc                                                                                       <  4:193802/149‑1 (MQ=255)

                                                                    aCCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAaaca                                                                                 >  1:399208/1‑148 (MQ=255)

                                                                        tGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAActg                                                                             <  2:330053/148‑1 (MQ=255)

                                                                           ccATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCTGCGACCATCGCCTCGAAACAACTGCTgg                                                                         >  4:578998/1‑149 (MQ=255)

                                                                            cATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTggc                                                                        >  4:712086/1‑149 (MQ=255)

                                                                              tatCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTggcg                                                                       <  2:261837/148‑1 (MQ=255)

                                                                                   ttGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcgc                                                                                                                                          <  3:163128/76‑1 (MQ=255)

                                                                                   ttGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATcgc                                                                                                                                          >  4:163128/1‑76 (MQ=255)

                                                                                   ttGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTa                                                                 >  3:381524/1‑149 (MQ=255)

                                                                                     gAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAg                                                                <  2:393601/148‑1 (MQ=255)

                                                                                      aaCTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCt                                                                                                                        >  2:355104/1‑91 (MQ=255)

                                                                                      aaCTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCt                                                                                                                        <  1:355104/91‑1 (MQ=255)

                                                                                        cTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCa                                                            >  4:107649/1‑149 (MQ=255)

                                                                                         ttGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCag                                                           >  4:692837/1‑149 (MQ=255)

                                                                                         ttGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCag                                                           >  4:203588/1‑149 (MQ=255)

                                                                                         ttGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCag                                                           >  4:120403/1‑149 (MQ=255)

                                                                                            ttCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCagag                                                         <  1:408369/148‑1 (MQ=255)

                                                                                            ttCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCagag                                                         >  2:812455/1‑148 (MQ=255)

                                                                                             tCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGt                                                        >  1:255104/1‑148 (MQ=255)

                                                                                             tCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTc                                                       <  4:838239/149‑1 (MQ=255)

                                                                                               gCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGa                                                     <  3:242417/149‑1 (MQ=255)

                                                                                                ccAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGa                                                     >  1:645075/1‑148 (MQ=255)

                                                                                                ccAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGAc                                                    >  3:241253/1‑149 (MQ=255)

                                                                                                ccAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGAc                                                    <  4:381524/149‑1 (MQ=255)

                                                                                                   gTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGa                                                 <  3:712086/149‑1 (MQ=255)

                                                                                                    tGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGAt                                                <  4:885968/149‑1 (MQ=255)

                                                                                                     gcgcgcAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCTGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGAt                                                <  2:703831/148‑1 (MQ=255)

                                                                                                        cgcAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCa                                             >  1:623525/1‑148 (MQ=255)

                                                                                                           aGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACtg                                         >  3:235066/1‑149 (MQ=255)

                                                                                                            gTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACtg                                         <  2:22324/148‑1 (MQ=255)

                                                                                                             ttACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCaa                                            >  3:115755/1‑144 (MQ=255)

                                                                                                             ttACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCaa                                            <  4:115755/144‑1 (MQ=255)

                                                                                                             ttACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACtgtg                                       <  3:300676/149‑1 (MQ=255)

                                                                                                                 aTGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCGGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGaa                                   >  4:890564/1‑149 (MQ=255)

                                                                                                                 aTGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGaa                                   >  4:275251/1‑149 (MQ=255)

                                                                                                                 aTGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGa                                    <  1:812455/148‑1 (MQ=255)

                                                                                                                  tGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAAc                                  <  3:120403/149‑1 (MQ=255)

                                                                                                                    ccgccgGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACa                                 <  1:842170/148‑1 (MQ=255)

                                                                                                                           gCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGttt                          <  2:507408/148‑1 (MQ=255)

                                                                                                                              tCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGa                      <  3:874670/149‑1 (MQ=255)

                                                                                                                              tCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGa                      >  3:608881/1‑149 (MQ=255)

                                                                                                                                  aTTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGAttt                   <  2:309507/148‑1 (MQ=255)

                                                                                                                                   ttGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTAt                 <  3:111750/149‑1 (MQ=255)

                                                                                                                                      ttATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGc              >  3:3031/1‑149 (MQ=255)

                                                                                                                                       tATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCa                                                                                              >  3:551145/1‑68 (MQ=255)

                                                                                                                                       tATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCa                                                                                              <  4:551145/68‑1 (MQ=255)

                                                                                                                                       tATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGc              >  1:58868/1‑148 (MQ=255)

                                                                                                                                        aTTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCAt            <  4:367525/149‑1 (MQ=255)

                                                                                                                                         ttGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCAt            <  2:255104/148‑1 (MQ=255)

                                                                                                                                          tgtgACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCAtt           >  2:72073/1‑148 (MQ=255)

                                                                                                                                               cGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGt      <  1:441225/148‑1 (MQ=255)

                                                                                                                                                gTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTAt                 >  4:552167/1‑136 (MQ=255)

                                                                                                                                                gTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTAt                 <  3:552167/136‑1 (MQ=255)

                                                                                                                                                  cGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCggg  >  4:114598/1‑149 (MQ=255)

GAGTGACAGCCAGACGCTGGTGGTAAAACTCGGCACCAGTGTGCTAACAGGCGGATCGCGCCGTCTGAACCGTGCCCATATCGTTGAACTTGTTCGCCAGTGCGCGCAGTTACATGCCGCCGGGCATCGGATTGTTATTGTGACGTCGGGC‑GCGATCGCCGCCGGACGTGAGCACCTGGGTTACCCGGAACTGCCAGCGACCATCGCCTCGAAACAACTGCTGGCGGCGGTAGGGCAGAGTCGACTGATTCAACTGTGGGAACAGCTGTTTTCGATTTATGGCATTCACGTCGGG  >  NC_000913/260390‑260684

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