BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,357,063	G→A	5.7%	A38A (GCC→GCT)	PP_2950 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,357,063	0	G	A	5.7%	991.0 / 31.6	314	A38A (GCC→GCT)	PP_2950	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (167/129); new base A (15/3); total (182/132)
Fisher's exact test for biased strand distribution p-value = 2.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 321 total aligned reads displayed.

ATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTTGGG  >  NC_002947/3356914‑3357212

      ttGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATc                                                                                                                                                 >  1:585772/1‑150 (MQ=255)

          cGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCATACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCtttt                                                                                                                                             <  2:2835404/150‑1 (MQ=255)

           ggTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTc                                                                                                                                            <  2:4255422/150‑1 (MQ=255)

           ggTGAAGGCCCAGAGCACATTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTc                                                                                                                                            <  1:4771186/150‑1 (MQ=255)

            gTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCg                                                                                                                                           >  1:4968927/1‑150 (MQ=255)

             tGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCgg                                                                                                                                          <  2:5523737/150‑1 (MQ=255)

                    ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                                                   >  1:101337/1‑150 (MQ=255)

                    ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                                                   >  2:420494/1‑150 (MQ=255)

                    ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                                                   >  2:4892600/1‑150 (MQ=255)

                    ccAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTcgcg                                                                                                                                   >  2:1461457/1‑150 (MQ=255)

                     cAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCgg                                                                                                                                  <  2:8742533/150‑1 (MQ=255)

                         gCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTAt                                                                                                                              >  1:7413657/1‑150 (MQ=255)

                           tctTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTc                                                                                                                            <  2:7648900/150‑1 (MQ=255)

                               tGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGc                                                                                                                        <  1:7580482/150‑1 (MQ=255)

                                  tCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCAt                                                                                                                     >  2:6299684/1‑150 (MQ=255)

                                    aaCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATcc                                                                                                                   >  2:8329773/1‑150 (MQ=255)

                                    aaCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATcc                                                                                                                   >  1:2875830/1‑150 (MQ=255)

                                    aaCCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATcc                                                                                                                   >  1:6370372/1‑150 (MQ=255)

                                      ccGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTc                                                                                                                 <  1:6461282/150‑1 (MQ=255)

                                       cGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCa                                                                                                                >  2:6830836/1‑150 (MQ=255)

                                        gATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAg                                                                                                               >  2:8904025/1‑150 (MQ=255)

                                        caTAGTCCAGTTCTTTCTTCCACTTGCGAGTTACGAAGTCCCGAGCTTCTACACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAg                                                                                                               <  2:8271133/149‑1 (MQ=255)

                                         aTAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGc                                                                                                              >  1:914336/1‑150 (MQ=255)

                                          tactCCAGTTCCTTATTCCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGca                                                                                                             >  1:4254847/4‑149 (MQ=255)

                                            gTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAg                                                                                                           >  2:4383047/1‑150 (MQ=255)

                                             tCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGt                                                                                                          <  1:7800667/150‑1 (MQ=255)

                                             tCCAGTTCCTTATTCCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAgg                                                                                                          >  1:534974/1‑149 (MQ=255)

                                             tCCAGTTCCTTATTCCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAgg                                                                                                          >  1:5070848/1‑149 (MQ=255)

                                                     tttcttCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtt                                                                                                  >  2:1190035/1‑150 (MQ=255)

                                                     tttcttCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtt                                                                                                  >  2:8635875/1‑150 (MQ=255)

                                                     tttcttCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGACAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtt                                                                                                  >  2:533153/1‑150 (MQ=255)

                                                     cttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCtt                                                                                                  >  1:464019/5‑150 (MQ=255)

                                                      ttcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCttt                                                                                                 <  2:3440458/150‑1 (MQ=255)

                                                      ttcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCttt                                                                                                 <  1:2664431/150‑1 (MQ=255)

                                                      ttattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCttt                                                                                                 >  2:3987359/4‑150 (MQ=255)

                                                       tcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtttt                                                                                                <  2:159913/150‑1 (MQ=255)

                                                       tcttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtttt                                                                                                <  2:7090428/150‑1 (MQ=255)

                                                       tattcCACTTGCGGGTCAGGAAGGCTTCGGCTTCGGCACGGTTTGTAACGAATTCAACGCGGTTGCGCCACCGCGTGTCCAAGGTATAGGCGAGAGCGACTTTCGCGGGATCGCGGGTATTCCAGCCATCTTCCGCAAGGCGGACCtttt                                                                                                <  2:77904/148‑1 (MQ=255)

                                                         ttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCa                                                                                              >  2:9088814/1‑150 (MQ=255)

                                                         ttcCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCa                                                                                              >  2:3652628/1‑150 (MQ=255)

                                                         ttcCACTTGCGAGTTAGGAAGTCCAGAGCTTCTTCACGGTTGGTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGAATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGACATCCTCAGCGAGTCGGACCTTTTCa                                                                                              >  1:8147628/1‑150 (MQ=255)

                                                            cACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATc                                                                                           >  2:2577482/1‑150 (MQ=255)

                                                            cACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATc                                                                                           <  2:139986/150‑1 (MQ=255)

                                                                  cGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGag                                                                                     >  1:2437606/1‑149 (MQ=255)

                                                                   gAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGAtt                                                                                    >  2:2354668/1‑150 (MQ=255)

                                                                      ttAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTcgc                                                                                 >  2:5006492/1‑150 (MQ=255)

                                                                      ttAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTcgc                                                                                 <  2:914336/150‑1 (MQ=255)

                                                                        aGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGa                                                                               <  1:6214588/150‑1 (MQ=255)

                                                                          gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtttt                                                                                                <  2:6636007/131‑1 (MQ=255)

                                                                          gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCtttt                                                                                                >  1:6636007/1‑131 (MQ=255)

                                                                          gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGt                                                                             >  1:4595796/1‑150 (MQ=255)

                                                                          gAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGt                                                                             <  1:5171175/150‑1 (MQ=255)

                                                                           aaGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTg                                                                            <  1:5624508/150‑1 (MQ=255)

                                                                              tCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGaaa                                                                         >  2:2047827/1‑150 (MQ=255)

                                                                              tCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGaaa                                                                         <  2:4560067/150‑1 (MQ=255)

                                                                               cccGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAag                                                                        >  2:6421660/1‑150 (MQ=255)

                                                                               cccGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAag                                                                        <  2:2875830/150‑1 (MQ=255)

                                                                                ccGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTAt                                                                                                                              >  2:182504/1‑95 (MQ=255)

                                                                                ccGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAagg                                                                       >  2:7271746/1‑150 (MQ=255)

                                                                                 cGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAagga                                                                      <  1:8494567/150‑1 (MQ=255)

                                                                                        cttcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTg                                                               <  1:8837307/150‑1 (MQ=255)

                                                                                         ttcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTgg                                                              >  1:719614/1‑150 (MQ=255)

                                                                                         ttcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTgg                                                              >  2:3784077/1‑150 (MQ=255)

                                                                                          tcACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGc                                                             <  2:2821487/150‑1 (MQ=255)

                                                                                           cACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCg                                                            <  2:4599515/150‑1 (MQ=255)

                                                                                                ttgttTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGa                                                       <  2:3082489/150‑1 (MQ=255)

                                                                                                ttgttTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGa                                                       <  1:7234567/150‑1 (MQ=255)

                                                                                                  gttTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAgt                                                     <  1:1234662/150‑1 (MQ=255)

                                                                                                  gttTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAgt                                                     <  1:1800280/150‑1 (MQ=255)

                                                                                                    ttGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAgtgt                                                   <  1:8764398/150‑1 (MQ=255)

                                                                                                        aaaCTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTc                                               <  1:8012039/150‑1 (MQ=255)

                                                                                                        aaaCTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTc                                               >  1:2926160/1‑150 (MQ=255)

                                                                                                        aaaCTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTc                                               >  1:8997268/1‑150 (MQ=255)

                                                                                                            tCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTg                                           >  1:7679237/1‑150 (MQ=255)

                                                                                                                 aCGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACAt                                      >  2:3297239/1‑150 (MQ=255)

                                                                                                                    gTTACGCCAGTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACAtggt                                   >  1:2114349/1‑150 (MQ=255)

                                                                                                                     ttACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACAtggtg                                  <  2:2659993/150‑1 (MQ=255)

                                                                                                                       aCGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATggtggt                                >  1:3218386/1‑150 (MQ=255)

                                                                                                                            aTTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTc                           <  2:7552385/150‑1 (MQ=255)

                                                                                                                                    tCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCttt                   >  2:5626916/1‑150 (MQ=255)

                                                                                                                                    tCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCttt                   <  1:99371/150‑1 (MQ=255)

                                                                                                                                    tCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCttt                   >  2:1804345/1‑150 (MQ=255)

                                                                                                                                    tCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCttt                   >  2:4791107/1‑150 (MQ=255)

                                                                                                                                    tCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCttt                   >  1:3450560/1‑150 (MQ=255)

                                                                                                                                     ccACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGTCCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCtttt                  <  2:7977398/150‑1 (MQ=255)

                                                                                                                                     ccACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCtttt                  <  2:2256072/150‑1 (MQ=255)

                                                                                                                                      cACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAggagga                                                                   >  2:5969025/1‑100 (MQ=255)

                                                                                                                                      cACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTg                 >  2:7814518/1‑150 (MQ=255)

                                                                                                                                      cACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTg                 <  1:1565074/150‑1 (MQ=255)

                                                                                                                                      cACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTg                 >  1:487193/1‑150 (MQ=255)

                                                                                                                                        ccGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAg                                                                                                               <  2:6890400/54‑1 (MQ=255)

                                                                                                                                        ccGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAg                                                                                                               >  1:6890400/1‑54 (MQ=255)

                                                                                                                                         cGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGAt              <  2:1906464/150‑1 (MQ=255)

                                                                                                                                          gTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATg             <  2:5893172/150‑1 (MQ=255)

                                                                                                                                           tataGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATgg            >  1:8332086/1‑150 (MQ=255)

                                                                                                                                            ataGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGt           >  1:5430471/1‑150 (MQ=255)

                                                                                                                                             taGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGtt          <  1:3884055/150‑1 (MQ=255)

                                                                                                                                               ggCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCa        <  2:5255898/150‑1 (MQ=255)

                                                                                                                                                  cAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGtt     >  1:4032287/1‑150 (MQ=255)

                                                                                                                                                  cAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGtt     <  1:1841223/150‑1 (MQ=255)

ATGCGGTTGTCGGTGAAGGCCCAGAGCTCTTTGATCAACCGATAGTCCAGTTCTTTCTTCCACTTGCGAGTTAGGAAGTCCCGAGCTTCTTCACGGTTGTTTGCAAACTCGGCACGGTTACGCCATTGGGTGTCCACCGTATAGGCCAGGGCAATCTTTTCGGGTTCGCGGGTATTCCAGCCATCCTCAGCGAGTCGGACCTTTTCAATCGCCGATTCGCGAGTGAAAGGAGGAAGTGGCGGACGAGTGTGTTCGTTGGACATGGTGGTGCCTCCAATCTTTTGGATGGTTCAGTTGGG  >  NC_002947/3356914‑3357212

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