BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,615,947	G→T	88.4%	E73D (GAG→GAT)	radA →	DNA repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,615,947	0	G	T	88.4%	761.7 / 100.1	300	E73D (GAG→GAT)	radA	DNA repair protein
Reads supporting (aligned to +/- strand): ref base G (19/16); new base T (128/137); total (147/153)
Fisher's exact test for biased strand distribution p-value = 5.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 302 total aligned reads displayed.

TGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGA  >  CP009273/4615801‑4616095

    gtgcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCcg                                                                                                                                                >  2:1139723/1‑149 (MQ=255)

    gtgcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCcg                                                                                                                                                >  1:1773243/1‑149 (MQ=255)

    gtgcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCcg                                                                                                                                                <  2:3469409/149‑1 (MQ=255)

      gcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCcgcg                                                                                                                                              <  2:2865882/149‑1 (MQ=255)

        cTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGtt                                                                                                                                            >  2:2805631/1‑149 (MQ=255)

         tGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGttt                                                                                                                                           >  1:3179566/1‑149 (MQ=255)

           tCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGttttt                                                                                                                                         >  2:4736910/1‑149 (MQ=255)

           tCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGttttt                                                                                                                                         >  1:746124/1‑149 (MQ=255)

               gCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCAc                                                                                                                                     >  1:1394825/1‑149 (MQ=255)

                 cTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCg                                                                                                                                   >  2:2021569/1‑149 (MQ=255)

                  tGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCgg                                                                                                                                  >  1:4373299/1‑149 (MQ=255)

                   ggAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGa                                                                                                                                 >  2:2192239/1‑149 (MQ=255)

                   ggAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGa                                                                                                                                 >  1:1825161/1‑149 (MQ=255)

                   ggAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGa                                                                                                                                 >  1:824097/1‑149 (MQ=255)

                    gAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGAt                                                                                                                                <  2:1571841/149‑1 (MQ=255)

                      acacCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGAttt                                                                                                                              >  2:1801662/1‑149 (MQ=255)

                        acCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTaa                                                                                                                            <  1:669607/149‑1 (MQ=255)

                            tCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAagag                                                                                                                        >  1:4468094/1‑149 (MQ=255)

                            tCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAagag                                                                                                                        <  2:609837/149‑1 (MQ=255)

                             cACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGt                                                                                                                       <  1:1720436/149‑1 (MQ=255)

                             cACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGt                                                                                                                       <  1:4323557/149‑1 (MQ=255)

                                 gAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGa                                                                                                                   >  2:1270230/1‑149 (MQ=255)

                                  aGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGAc                                                                                                                  >  1:2816149/1‑149 (MQ=255)

                                  aGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGAc                                                                                                                  <  2:422465/149‑1 (MQ=255)

                                   ggTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGAcc                                                                                                                 >  1:2208975/1‑149 (MQ=255)

                                     tGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACcgc                                                                                                               >  2:4108740/1‑149 (MQ=255)

                                       cGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGt                                                                                                             <  1:3981585/149‑1 (MQ=255)

                                         tCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTAc                                                                                                           >  2:743540/1‑149 (MQ=255)

                                          cTTGCTGCGTCGCCATTGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACt                                                                                                          <  1:3254240/149‑1 (MQ=255)

                                          cTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACt                                                                                                          <  2:3696169/149‑1 (MQ=255)

                                           ttGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCCCCGGATTTAAAGAGTTCGACCGCGTACTa                                                                                                         >  1:2169140/1‑149 (MQ=255)

                                             gctgcGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAgg                                                                                                       >  2:2454022/1‑149 (MQ=255)

                                               tgcGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcg                                                                                                     >  2:1550728/1‑149 (MQ=255)

                                                  gTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcg                                                                                                  >  1:828171/1‑149 (MQ=255)

                                                  gTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcg                                                                                                  >  1:4035901/1‑149 (MQ=255)

                                                     gCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGgcggcggcg                                                                                               >  1:1199194/1‑149 (MQ=255)

                                                      ccAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCgt                                                                                              <  1:2643664/149‑1 (MQ=255)

                                                        aaTGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCgtgg                                                                                            >  2:4159871/1‑149 (MQ=255)

                                                         aTGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCgtggt                                                                                           >  2:2695009/1‑149 (MQ=255)

                                                         aTGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCgtggt                                                                                           >  1:1948252/1‑149 (MQ=255)

                                                           ggtggCGCTTAACTAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGc                                                                                         <  2:2638298/149‑1 (MQ=255)

                                                           ggtggCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGc                                                                                         <  2:584193/149‑1 (MQ=255)

                                                              ggCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAg                                                                                      >  1:2051637/1‑149 (MQ=255)

                                                               gcgcGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAgg                                                                                     >  1:1502636/1‑149 (MQ=255)

                                                                  cgTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAg                                                                                  <  1:3902736/149‑1 (MQ=255)

                                                                   gTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGCTCTGCCGCGTTTGTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGt                                                                                 <  1:384725/149‑1 (MQ=255)

                                                                     aaCGCGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGc                                                                               <  1:2376417/149‑1 (MQ=255)

                                                                     aaCGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGc                                                                               <  1:1379953/149‑1 (MQ=255)

                                                                       cGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCa                                                                             <  1:4272521/149‑1 (MQ=255)

                                                                       cGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCa                                                                             <  2:336336/149‑1 (MQ=255)

                                                                        gAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCAt                                                                            <  2:4289814/149‑1 (MQ=255)

                                                                        gAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCAt                                                                            >  1:3529844/1‑149 (MQ=255)

                                                                         aGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCAtt                                                                           >  1:2522075/1‑149 (MQ=255)

                                                                            gTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTg                                                                        >  2:2313966/1‑149 (MQ=255)

                                                                             tctcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                       >  2:3317225/1‑149 (MQ=255)

                                                                             tctcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                       >  1:3388998/1‑149 (MQ=255)

                                                                             tctcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                       >  2:2964323/1‑149 (MQ=255)

                                                                             tctcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                       >  1:3644098/1‑149 (MQ=255)

                                                                                cAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTg                                                                    <  1:2914351/149‑1 (MQ=255)

                                                                                cAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTg                                                                    >  1:1013632/1‑149 (MQ=255)

                                                                                     gCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGt                                                               <  1:1569974/149‑1 (MQ=255)

                                                                                         tGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAcc                                                           >  1:4298052/1‑149 (MQ=255)

                                                                                           ccGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCt                                                         <  1:965904/149‑1 (MQ=255)

                                                                                             ggTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTgg                                                       <  1:3442119/149‑1 (MQ=255)

                                                                                                 gcgcCGGGGTGGCAAACGTCCAGAAACTCTCCGATATCAGCCTTGCAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCg                                                   <  2:3853398/149‑1 (MQ=255)

                                                                                                       gggTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGaaa                                             <  2:1816547/149‑1 (MQ=255)

                                                                                                        ggTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAAt                                            <  2:3223034/149‑1 (MQ=255)

                                                                                                         gTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATc                                           <  1:642791/149‑1 (MQ=255)

                                                                                                           ggCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCa                                         <  1:3604013/149‑1 (MQ=255)

                                                                                                              aaaaGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgc                                      <  1:1866467/149‑1 (MQ=255)

                                                                                                                  gTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctgct                                  <  2:486686/149‑1 (MQ=255)

                                                                                                                  gTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTGAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctgct                                  <  2:2927197/149‑1 (MQ=255)

                                                                                                                    ccAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTAc                                <  1:2992683/149‑1 (MQ=255)

                                                                                                                       gAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGc                             >  1:4199064/1‑149 (MQ=255)

                                                                                                                          aCTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGACATCCACGCTGCTACTGCaaa                          >  1:2504378/1‑149 (MQ=255)

                                                                                                                          aCTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCaaa                          <  2:4437110/149‑1 (MQ=255)

                                                                                                                          aCTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATGCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCaaa                          >  2:3103109/1‑149 (MQ=255)

                                                                                                                            tctcCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACg                        >  2:2354891/1‑149 (MQ=255)

                                                                                                                            tctcCGATATCAGCCTTGAAGATCTGCCGCGTTTGTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACg                        <  2:4523218/149‑1 (MQ=255)

                                                                                                                             ctcCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGc                       <  2:1769820/149‑1 (MQ=255)

                                                                                                                             ctcCGATATCAGCCATGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGc                       >  1:2424681/1‑149 (MQ=255)

                                                                                                                             ctcCGATATCAGCCATGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGc                       >  1:133425/1‑149 (MQ=255)

                                                                                                                                 gATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCtgtg                   >  2:671841/1‑149 (MQ=255)

                                                                                                                                 gATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCtgtg                   <  2:3793496/149‑1 (MQ=255)

                                                                                                                                  atatCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGc                  >  2:3229604/1‑149 (MQ=255)

                                                                                                                                      cAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTTTGCAAAc              <  2:2528416/149‑1 (MQ=255)

                                                                                                                                      cAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAAc              >  2:2742067/1‑149 (MQ=255)

                                                                                                                                      cAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAAc              >  2:350649/1‑149 (MQ=255)

                                                                                                                                      cAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAAc              <  1:1487166/149‑1 (MQ=255)

                                                                                                                                       aGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACt             <  1:4762009/149‑1 (MQ=255)

                                                                                                                                        gCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTg            <  1:1677939/149‑1 (MQ=255)

                                                                                                                                         ccTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTgg           <  1:925355/149‑1 (MQ=255)

                                                                                                                                          cTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGc          >  1:605814/1‑149 (MQ=255)

                                                                                                                                           ttGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGcc         <  2:4008545/149‑1 (MQ=255)

                                                                                                                                           ttGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGcc         >  2:2289363/1‑149 (MQ=255)

                                                                                                                                           ttGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGcc         >  2:1840839/1‑149 (MQ=255)

                                                                                                                                           ttGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGGTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGcc         <  1:458338/149‑1 (MQ=255)

                                                                                                                                            tGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGccc        <  1:4682011/149‑1 (MQ=255)

                                                                                                                                             gAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCca       >  2:949287/1‑149 (MQ=255)

                                                                                                                                             gAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCca       >  1:1058627/1‑149 (MQ=255)

TGCAGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGA  >  CP009273/4615801‑4616095

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: `‑`