BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	537,464	T→A	6.7%	noncoding (57/66 nt)	PP_mr07 →	P26 ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	537,464	0	T	A	6.7%	792.0 / 30.7	237	noncoding (57/66 nt)	PP_mr07	P26 ncRNA
Reads supporting (aligned to +/- strand): ref base T (98/123); new base A (10/6); total (108/129)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 251 total aligned reads displayed.

TCAGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGAC  >  NC_002947/537315‑537608

tcAGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTAt                                                                                                                                                  >  1:608398/1‑150 (MQ=255)

tcAGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTAt                                                                                                                                                  >  2:618736/1‑150 (MQ=255)

tcAGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTAt                                                                                                                                                  >  2:4712886/1‑150 (MQ=255)

tcAGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTAt                                                                                                                                                  >  2:4162845/1‑150 (MQ=255)

 cAGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATg                                                                                                                                                 <  1:3074687/150‑1 (MQ=255)

    gTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCa                                                                                                                                              >  2:1769289/1‑150 (MQ=255)

    gTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCa                                                                                                                                              <  1:2237347/150‑1 (MQ=255)

        tAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCgg                                                                                                                                          >  2:1060468/1‑150 (MQ=255)

         aGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGc                                                                                                                                         >  1:2635290/1‑150 (MQ=255)

          gCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCACCAGcc                                                                                                                                        <  2:1126108/150‑1 (MQ=255)

           cTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCt                                                                                                                                       <  2:3219850/150‑1 (MQ=255)

            tGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCtt                                                                                                                                      >  2:3248777/1‑150 (MQ=255)

            tGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCCGGTGGCATAAGCCACCAGccat                                                                                                                                      <  2:1458616/150‑3 (MQ=255)

             gAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCttt                                                                                                                                     >  1:3666251/1‑150 (MQ=255)

                ggCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTc                                                                                                                                  <  1:1370355/150‑1 (MQ=255)

                 gCGTTTCGAAAGATGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTcc                                                                                                                                 <  1:5438715/150‑1 (MQ=255)

                   gTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGt                                                                                                                               >  1:98352/1‑150 (MQ=255)

                     ttCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTa                                                                                                                             >  2:5319596/1‑150 (MQ=255)

                      tCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTAt                                                                                                                            <  1:50638/150‑1 (MQ=255)

                        gAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtg                                                                                                                          >  1:4131846/1‑150 (MQ=255)

                         aaaGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtgg                                                                                                                         >  1:5605788/1‑150 (MQ=255)

                         aaaGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtgg                                                                                                                         >  1:426973/1‑150 (MQ=255)

                          aagaagCCGCTGAAGACGCTAAGAAGAAGCCGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtggt                                                                                                                        <  1:4509093/150‑1 (MQ=255)

                            gaagCCGCTGCAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtggtgg                                                                                                                      <  1:3917515/150‑1 (MQ=255)

                            gaagCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATtggtgg                                                                                                                      <  1:905879/150‑1 (MQ=255)

                               gCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCg                                                                                                                   <  1:432230/150‑1 (MQ=255)

                                 cGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGt                                                                                                                 >  1:2799779/1‑150 (MQ=255)

                                 cGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGt                                                                                                                 >  1:1825863/1‑150 (MQ=255)

                                  gCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTc                                                                                                                >  2:2454220/1‑150 (MQ=255)

                                      taGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGc                                                                                                            <  1:5525453/149‑1 (MQ=255)

                                      aaGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGc                                                                                                            <  2:4887340/150‑1 (MQ=255)

                                        gACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCg                                                                                                          <  2:4010128/150‑1 (MQ=255)

                                           gCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTc                                                                                                       >  1:1858420/1‑150 (MQ=255)

                                            cTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTcc                                                                                                      <  1:3414809/150‑1 (MQ=255)

                                             tAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTAtt                                                                                                                           >  2:4613837/1‑128 (MQ=255)

                                             tAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTAtt                                                                                                                           <  1:4613837/128‑1 (MQ=255)

                                             tAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTccc                                                                                                     >  1:2931133/1‑150 (MQ=255)

                                             tAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTccc                                                                                                     >  1:484385/1‑150 (MQ=255)

                                                 aagaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGaga                                                                                                 >  1:479164/1‑148 (MQ=255)

                                                 aagaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATa                                                                                                 >  1:130330/1‑150 (MQ=255)

                                                  agaagCTGGAAGATGCTGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATaa                                                                                                <  2:5402080/150‑1 (MQ=255)

                                                  agaagCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATaa                                                                                                <  1:291428/150‑1 (MQ=255)

                                                      gCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAg                                                                                            >  2:2780989/1‑150 (MQ=255)

                                                        tGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCt                                                                                                                                       <  2:3296439/105‑1 (MQ=255)

                                                        tGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCt                                                                                          >  2:722897/1‑150 (MQ=255)

                                                          gaagaaGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGc                                                                                        <  1:70967/150‑1 (MQ=255)

                                                            agaagCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCaa                                                                                      <  2:920058/150‑1 (MQ=255)

                                                            agaagCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCaa                                                                                      <  1:343345/150‑1 (MQ=255)

                                                                 cAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGacac                                                                                 <  2:2662291/150‑1 (MQ=255)

                                                                 cAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGacac                                                                                 <  2:2038426/150‑1 (MQ=255)

                                                                      gcTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAtg                                                                            >  2:373210/1‑150 (MQ=255)

                                                                      gcTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTCGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAtg                                                                            >  1:3809469/1‑150 (MQ=255)

                                                                       cTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAtgt                                                                           <  1:1060468/150‑1 (MQ=255)

                                                                       cTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAtgt                                                                           >  2:3800248/1‑150 (MQ=255)

                                                                         aaaGTCGAGCTGAAGTAATCTCGACCGTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATgtgg                                                                         <  2:4505345/150‑1 (MQ=255)

                                                                              cGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGgc                                                                    <  1:1913729/150‑1 (MQ=255)

                                                                                aGCTGACGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTGTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACTAGCTGCAAGACACCCATGTGGGTGgcgc                                                                  <  1:5003637/150‑1 (MQ=255)

                                                                                       gTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAATCCaa                                                           <  2:58066/150‑1 (MQ=255)

                                                                                        tAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAg                                                          <  2:3449870/150‑1 (MQ=255)

                                                                                         aaTCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGTCTGATGGCTGGTGGCTTATGCCACCGTCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTACCGATAACGAGGTGCAAGACACCCATGTGGGTGGCGCAAACCAAgg                                                         <  2:1859313/150‑1 (MQ=255)

                                                                                         aaTCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAgg                                                         >  2:2576920/1‑150 (MQ=255)

                                                                                          aTCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAggg                                                        <  1:4651389/150‑1 (MQ=255)

                                                                                           tctcGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAgggg                                                       >  1:802859/1‑150 (MQ=255)

                                                                                            ctcGACGTTGCGATGACAGCCTGAGCGTTGAGCGATAGGCTGATGGCTGGTGTCTTATGCCACCGGACTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGTTGGCGCAAACATAGGGGt                                                      <  2:123175/150‑1 (MQ=255)

                                                                                              cGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGttt                                                    <  1:2783864/150‑1 (MQ=255)

                                                                                                 ccTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCa                                                 <  1:2446641/150‑1 (MQ=255)

                                                                                                  cTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCAc                                                >  2:2906279/1‑150 (MQ=255)

                                                                                                  cTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTGTGCAc                                                >  1:4045382/1‑150 (MQ=255)

                                                                                                  cTTCCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCAc                                                >  1:2391870/1‑150 (MQ=255)

                                                                                                   ttGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACg                                               >  1:826848/1‑150 (MQ=255)

                                                                                                   ttGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACg                                               >  2:3473574/1‑150 (MQ=255)

                                                                                                           ccAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCATTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTg                                       <  2:3956191/150‑1 (MQ=255)

                                                                                                            cAGCCTGTGCGTTGTGCGAAAGGGTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTgg                                      <  1:1188943/150‑1 (MQ=255)

                                                                                                              gCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCTAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGgct                                    <  2:2591579/150‑1 (MQ=255)

                                                                                                                      gTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCccgc                            <  1:4095918/150‑1 (MQ=255)

                                                                                                                        tGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCccgccg                          <  1:1606029/150‑1 (MQ=255)

                                                                                                                         gcgCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCGGGCTGCTCCCGCCgg                         <  1:4069306/148‑1 (MQ=255)

                                                                                                                           gCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTAAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGga                       <  1:3011876/150‑1 (MQ=255)

                                                                                                                              aaaGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                    <  1:5552048/150‑1 (MQ=255)

                                                                                                                              aaaGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                    <  2:4414777/150‑1 (MQ=255)

                                                                                                                                 ggCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGcc                 <  2:2218627/150‑1 (MQ=255)

                                                                                                                                 ggCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGcc                 >  1:4256016/1‑150 (MQ=255)

                                                                                                                                 ggCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATGGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGcc                 <  2:2228139/150‑1 (MQ=255)

                                                                                                                                  gCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCa                <  1:3884419/150‑1 (MQ=255)

                                                                                                                                    tGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCaaa              >  2:802141/1‑150 (MQ=255)

                                                                                                                                    ggATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCaaa              <  2:2984406/149‑1 (MQ=255)

                                                                                                                                     aaaggCCGGTGGCATAAGCCACCAGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAAc             >  1:2871081/4‑150 (MQ=255)

                                                                                                                                      aTGGCTGGTGGCTTATGCCACCGGCCTGGTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTGTCTGGCTGCTCCCGCCGGGAGAAGCCAAACa            >  2:3762521/1‑150 (MQ=255)

                                                                                                                                       aggCCGGTGGCATAAGCCACCAGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAt                                                                             <  1:5113538/83‑1 (MQ=38)

                                                                                                                                       aggCCGGTGGCATAAGCCACCAGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCAt                                                                             >  2:5113538/2‑84 (MQ=38)

                                                                                                                                        ggCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg          <  1:5319596/150‑1 (MQ=255)

                                                                                                                                        ggCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg          <  1:2196730/150‑1 (MQ=255)

                                                                                                                                        ggCCGGTGGCATAAGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg          >  2:1294866/1‑150 (MQ=255)

                                                                                                                                        ggCCGGTGGCATAAGCCACCAGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg          >  1:1987941/1‑150 (MQ=255)

                                                                                                                                          cTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACAAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCa        <  1:3271005/150‑1 (MQ=255)

                                                                                                                                           tggtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAg       >  1:2574163/1‑150 (MQ=255)

                                                                                                                                           tggtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAg       >  1:5226018/1‑150 (MQ=255)

                                                                                                                                             gtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGt     <  1:4712886/150‑1 (MQ=255)

                                                                                                                                               ggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGa   <  1:160473/150‑1 (MQ=255)

                                                                                                                                               ggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGa   <  1:490448/150‑1 (MQ=255)

TCAGGTTGTAGCTGAAGGCGTTTCGAAAGAAGCCGCTGAAGACGCTAAGAAGAAGCTGGAAGAAGCAGGCGCTAAAGTCGAGCTGAAGTAATCTCGACCTTGCGATGCCAGCCTGAGCGTTGAGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGAC  >  NC_002947/537315‑537608

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