BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,788,401	G→C	5.7%	R78R (CGC→CGG)	PP_3350 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,788,401	0	G	C	5.7%	333.3 / 7.9	105	R78R (CGC→CGG)	PP_3350	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (46/53); new base C (2/4); total (48/57)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.73e-01

Only 100 of 105 total aligned reads displayed.

GGGTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGC  >  NC_002947/3788256‑3788542

gggTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCgg                                                                                                                                            <  7:294675/149‑1 (MQ=255)

gggTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCgg                                                                                                                                            <  2:204413/149‑1 (MQ=255)

    aTTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCt                                                                                                                                        >  7:7348/1‑149 (MQ=255)

      tctcGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCt                                                                                                                                        >  6:373253/1‑147 (MQ=255)

       ctcGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCa                                                                                                                                      <  7:49153/148‑1 (MQ=255)

       ctcGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGCTACCGGGccgat                                                                                                                                      <  8:74372/148‑4 (MQ=255)

        tcGACCTGGAAGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAga                                                                                                                                    >  7:171869/1‑148 (MQ=255)

          gACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGCTACCGGGccgatg                                                                                                                                     <  8:315065/146‑5 (MQ=255)

          gACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGCTACCGGGccgatg                                                                                                                                     <  1:265270/146‑5 (MQ=255)

             cTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGcccc                                                                                                                                 >  5:380786/1‑147 (MQ=255)

                gACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGtt                                                                                                                            <  6:113828/149‑1 (MQ=255)

                gACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGt                                                                                                                             <  3:17146/148‑1 (MQ=255)

                  cGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGc                                                                                                                          <  2:61466/149‑1 (MQ=255)

                   gAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTg                                                                                                                           >  7:68597/1‑147 (MQ=255)

                   gAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCa                                                                                                                         <  4:35030/149‑1 (MQ=255)

                     aGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAgt                                                                                                                       >  5:138752/1‑148 (MQ=255)

                      gggaCGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAgtc                                                                                                                      <  3:139016/146‑3 (MQ=255)

                           gCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGcc                                                                                                                 >  8:194596/1‑149 (MQ=255)

                             aCGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCg                                                                                                                >  5:303395/1‑148 (MQ=255)

                             aCGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGa                                                                                                               >  3:272919/1‑149 (MQ=255)

                                 cATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGActct                                                                                                           <  1:280815/149‑1 (MQ=255)

                                  aTCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTg                                                                                                          <  8:7348/149‑1 (MQ=255)

                                     gtcgAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCTTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACa                                                                                                       <  6:273751/149‑1 (MQ=255)

                                     gtcgAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACa                                                                                                       <  5:373253/149‑1 (MQ=255)

                                     gtcgAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACa                                                                                                       >  7:288616/1‑149 (MQ=255)

                                      tcgAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTCCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAg                                                                                                      >  8:211140/1‑149 (MQ=255)

                                        gAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGcc                                                                                                    <  6:125248/149‑1 (MQ=255)

                                        gAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGCTACCGGGccgatgc                                                                                                                                    <  3:151376/117‑6 (MQ=255)

                                        gAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGCTACCGGGccgatgc                                                                                                                                    >  4:151376/1‑112 (MQ=255)

                                         aaCAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGcc                                                                                                    >  5:170434/1‑148 (MQ=255)

                                          aCAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGAc                                                                                                              <  3:282129/137‑1 (MQ=255)

                                          aCAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGAc                                                                                                              >  4:282129/1‑137 (MQ=255)

                                           cAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGAt                                                                                                 >  7:184308/1‑149 (MQ=255)

                                            aGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATg                                                                                                <  8:68597/149‑1 (MQ=255)

                                             ggTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATgg                                                                                               >  7:79895/1‑149 (MQ=255)

                                               tGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTa                                                                                             <  6:117030/149‑1 (MQ=255)

                                                   ggTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAgcgc                                                                                         <  1:41003/149‑1 (MQ=255)

                                                    gTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAgcg                                                                                          >  5:157627/1‑147 (MQ=255)

                                                      cGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAgg                                                                                      <  6:157627/149‑1 (MQ=255)

                                                        gCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGACGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAgga                                                                                     >  6:127690/1‑147 (MQ=255)

                                                         cATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCGAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATg                                                                                                <  8:308978/136‑1 (MQ=255)

                                                         cATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCGAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATg                                                                                                >  7:308978/1‑136 (MQ=255)

                                                            gggCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGa                                                                                                  >  7:148687/1‑131 (MQ=255)

                                                            gggCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGa                                                                                                  <  8:148687/131‑1 (MQ=255)

                                                                cTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGCCAGCCGATGGTAGCGCAGGCCAGCGAAAg                                                                            >  4:245675/1‑149 (MQ=255)

                                                                cTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAg                                                                            >  1:291394/1‑149 (MQ=255)

                                                                   cAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAggg                                                                          <  3:28109/148‑1 (MQ=255)

                                                                    aGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGc                                                                        <  2:291394/149‑1 (MQ=255)

                                                                         ttGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTgg                                                                   >  1:329757/1‑149 (MQ=255)

                                                                         ttGCTGCAGATAGGAATCGTCGTCCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTgg                                                                   <  7:211140/149‑1 (MQ=255)

                                                                             tgcAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCt                                                                                                                                        <  2:223129/76‑1 (MQ=255)

                                                                             tgcAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCt                                                                                                                                        >  1:223129/1‑76 (MQ=255)

                                                                             tgcAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCtt                                                               >  3:49292/1‑149 (MQ=255)

                                                                              gcAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTg                                                              >  1:121586/1‑149 (MQ=255)

                                                                                aGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCAt                                                                                  <  7:319482/127‑1 (MQ=255)

                                                                                aGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCAt                                                                                  >  8:319482/1‑127 (MQ=255)

                                                                                 gATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCg                                                           <  1:225504/149‑1 (MQ=255)

                                                                                 gATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCAGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCg                                                           >  5:320579/1‑149 (MQ=255)

                                                                                  aTAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGGCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCgg                                                          >  8:244528/1‑149 (MQ=255)

                                                                                  aTAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCgg                                                          >  4:12139/1‑149 (MQ=255)

                                                                                    aGGATTCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCTACTCTTACAGCCGATGGTAGCGCAGGCCATc                                                                                 <  1:40248/124‑1 (MQ=255)

                                                                                    aGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCgggg                                                        >  3:138920/1‑149 (MQ=255)

                                                                                    aGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAACGGTTAGTCCGAAGCATCGGCCCGGTAGCGGGACTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGAAAGACGATGGTATCGCAGGCCATc                                                                                 >  2:40248/1‑124 (MQ=255)

                                                                                     ggAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCgggg                                                        <  8:288616/148‑1 (MQ=255)

                                                                                     ggAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCgggg                                                        <  7:194596/148‑1 (MQ=255)

                                                                                      gAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAgg                                                                           <  7:316329/128‑1 (MQ=255)

                                                                                      gAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAgg                                                                           >  8:316329/1‑128 (MQ=255)

                                                                                      gAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTg                                                                    <  5:340818/135‑1 (MQ=255)

                                                                                      gAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTg                                                                    >  6:340818/1‑135 (MQ=255)

                                                                                         tcgtcgTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtcgtc                                                   >  1:233649/1‑149 (MQ=255)

                                                                                           gtcgtTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtcgtc                                                   >  8:10167/1‑147 (MQ=255)

                                                                                                  aGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGATGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCTGTCAGGCCGACTCTGACAGGCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCg                                          <  2:213636/149‑1 (MQ=255)

                                                                                                   gCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCgg                                         >  5:10523/1‑149 (MQ=255)

                                                                                                   gCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCgg                                         >  2:201925/1‑149 (MQ=255)

                                                                                                    cgcgcgCATGCCAAAGAAGGGTTTGTCCGATGCATCGGCCGGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCgg                                         <  5:106104/146‑1 (MQ=255)

                                                                                                    cccGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCgg                                         >  1:69891/1‑148 (MQ=255)

                                                                                                    cccGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGt                                        <  4:272919/149‑1 (MQ=255)

                                                                                                     ccGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTa                                       <  8:143796/149‑1 (MQ=255)

                                                                                                     ccGCGCAGGCCAAAGAAGGGTTAGACAGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCAACGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCgg                                         >  7:235527/1‑147 (MQ=255)

                                                                                                        cgcAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTc                                    <  7:278634/149‑1 (MQ=255)

                                                                                                        cgcAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTc                                    <  5:127690/149‑1 (MQ=255)

                                                                                                         gcAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTcc                                   <  2:69891/149‑1 (MQ=255)

                                                                                                           aGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTc                                 <  8:79895/149‑1 (MQ=255)

                                                                                                            ggCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCg                                <  4:138920/149‑1 (MQ=255)

                                                                                                             gCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCgaga                                               >  7:299452/1‑133 (MQ=255)

                                                                                                             gCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCgaga                                               <  8:299452/133‑1 (MQ=255)

                                                                                                               cAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGa                               <  1:342728/147‑1 (MQ=255)

                                                                                                               cAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGa                               >  2:342728/1‑147 (MQ=255)

                                                                                                                       ggTTTGTCCGAAGCATGGGCCCGGTAGCGGGCCTCAGCCCCGTTTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAg                      <  8:235527/148‑1 (MQ=255)

                                                                                                                         tttGTCCGAAGCATCGGCCCGGTAGCGGGCCGCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAgggg                   <  8:184308/149‑1 (MQ=255)

                                                                                                                          ttGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGc                  <  7:244528/149‑1 (MQ=255)

                                                                                                                               cGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCt              <  6:138752/148‑1 (MQ=255)

                                                                                                                               cGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTg             >  5:339085/1‑149 (MQ=255)

                                                                                                                               cGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTg             >  8:41525/1‑149 (MQ=255)

                                                                                                                                   gCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCt         >  1:207986/1‑149 (MQ=255)

                                                                                                                                     tgaggCCCGCTACCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTcc       >  8:235440/4‑149 (MQ=255)

                                                                                                                                     aTCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTcc       <  3:12139/149‑1 (MQ=255)

                                                                                                                                      tCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTccc      <  6:170434/149‑1 (MQ=255)

                                                                                                                                         gCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCccgg   >  4:207329/1‑149 (MQ=255)

                                                                                                                                          cccGGTAGCGGGCCTCAGCCCCGAGTTGCAGTCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGGGGGGTCGTCGTGAATGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCccggc  <  2:329757/149‑1 (MQ=255)

GGGTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGC  >  NC_002947/3788256‑3788542

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 21 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`