BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	525,741	G→C	8.3%	E827D (GAG→GAC) ‡	rhsD →	protein RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	525,741	0	G	C	8.3%	352.4 / 18.9	108	E827D (GAG→GAC) ‡	rhsD	protein RhsD
Reads supporting (aligned to +/- strand): ref base G (54/45); new base C (4/5); total (58/50)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

Only 100 of 112 total aligned reads displayed.

ATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGT  >  NC_000913/525594‑525890

 tACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACg                                                                                                                                                    >  1:795512/1‑150 (MQ=255)

   cAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGt                                                                                                                                                  >  1:87783/1‑150 (MQ=255)

    aaCGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTg                                                                                                                                                 <  2:1498106/150‑1 (MQ=255)

     aCGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTCTCACGGTgc                                                                                                                                                <  1:863573/150‑1 (MQ=255)

      cGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTgcg                                                                                                                                               <  2:112804/150‑1 (MQ=255)

       gAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTCTCACGGTgcag                                                                                                                                              <  2:446003/150‑3 (MQ=255)

       gAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGGGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTCTCACGGTgcag                                                                                                                                              <  1:1322920/150‑3 (MQ=255)

        aGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCa                                                                                                                                             <  1:554171/150‑1 (MQ=255)

        aGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTCTCACGGTgcagc                                                                                                                                             <  1:489740/150‑4 (MQ=255)

        aGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTCTCACGGTgcagc                                                                                                                                             <  2:654098/150‑4 (MQ=255)

         gCAGGTGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGGAGCGGTTACCTGCAGGGAATCAAGCAGGGCGGGACGACGCTGGTCGAGTATACGCTGGACAGGATGCACCGGGAGCCGGTGCGCAg                                                                                                                                            >  2:366976/1‑150 (MQ=18)

                 gggCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGgca                                                                                                                                    <  2:443172/149‑1 (MQ=255)

                   gCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGgcagc                                                                                                                                  <  1:1662264/150‑1 (MQ=255)

                    cAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGgcagca                                                                                                                                 <  2:1401473/150‑1 (MQ=255)

                      aaCCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATg                                                                                                                               <  2:253117/150‑1 (MQ=255)

                        ccGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGTAGCATggc                                                                                                                             >  2:451484/1‑150 (MQ=255)

                        ccGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggc                                                                                                                             >  1:114257/1‑150 (MQ=255)

                        ccGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggc                                                                                                                             >  1:823538/1‑150 (MQ=255)

                         cgcgTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggca                                                                                                                            >  1:1034833/1‑150 (MQ=255)

                          gcgTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcag                                                                                                                           <  2:538277/150‑1 (MQ=255)

                            gTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggc                                                                                                                         <  1:664505/150‑1 (MQ=255)

                            gTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggc                                                                                                                         >  1:773967/1‑150 (MQ=255)

                            gTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggc                                                                                                                         <  2:286914/150‑1 (MQ=255)

                               aCGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGt                                                                                                                      <  1:534681/150‑1 (MQ=255)

                                cGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTgcgc                                                                                                                                              <  2:1229951/125‑1 (MQ=255)

                                cGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTgcgc                                                                                                                                              >  1:1229951/1‑125 (MQ=255)

                                cGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTa                                                                                                                     >  1:838396/1‑150 (MQ=255)

                                cGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTACACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTa                                                                                                                     >  1:775154/1‑150 (MQ=255)

                                       cAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGc                                                                                                              <  2:742380/150‑1 (MQ=255)

                                         gCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCAt                                                                                                            >  1:1107333/1‑150 (MQ=255)

                                          ccTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGGGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATa                                                                                                           <  1:163886/150‑1 (MQ=255)

                                           cTGCCGCCGTTGGAGTGGATAACGAATGGCAGCGGTTACCTGGCGGGAATAAACATGGCCGGGACGCCGATGGACGAGTATACGCGGGACAGGATCCACCGTGAGACGGTGCGCAGATTCGGAAGCATGGCAGGAAGTAATGCCGCATAc                                                                                                          >  2:1655664/1‑150 (MQ=18)

                                            tGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACg                                                                                                         >  2:311642/1‑150 (MQ=255)

                                              ccgccgGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGaa                                                                                                       <  1:850379/150‑1 (MQ=255)

                                               cgccgGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAAc                                                                                                      <  1:1251252/150‑1 (MQ=255)

                                               cgccgGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAAc                                                                                                      <  1:311642/150‑1 (MQ=255)

                                                gccgGTGGAGTGGCTGACCTATAGCAGCGGTTACCTCGCCCGAATTAAGCTGGGCGGGACGCCGATGGTCGAGTATACGCGGGACAGCCAGCACCGTTAGACGGTGCGCAGCTTCGGCATCATGGCAGGCAGTAATCCCGCATACGAAca                                                                                                     >  2:1567852/1‑149 (MQ=17)

                                                   ggTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGAAGCATGGCAGGCAGTAATGCCGCATACGAACTGAc                                                                                                  >  1:674807/1‑150 (MQ=255)

                                                    gTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGAcc                                                                                                 >  2:864422/1‑150 (MQ=255)

                                                     tGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCa                                                                                                <  2:823538/150‑1 (MQ=255)

                                                     tGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCa                                                                                                >  1:869365/1‑150 (MQ=255)

                                                           ggCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACAt                                                                                          >  2:1143084/1‑150 (MQ=255)

                                                             cTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATac                                                                                        >  1:1486320/1‑150 (MQ=255)

                                                                aCGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACAcc                                                                                     >  2:545484/1‑150 (MQ=255)

                                                                 cGTATGGCAGCGGTTACCTGGAGGGAATCAAGCTGGGCGGGACGACGCTGGTCGAGTATACGCGCGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACAccc                                                                                    >  2:781727/1‑150 (MQ=255)

                                                                  gTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACAcccc                                                                                   <  2:804097/150‑1 (MQ=255)

                                                                   tATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACAccccc                                                                                  >  1:896778/1‑150 (MQ=255)

                                                                      ggCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCa                                                                               >  1:77812/1‑150 (MQ=255)

                                                                      ggCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCa                                                                               >  1:787354/1‑150 (MQ=255)

                                                                      ggCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCa                                                                               <  1:1143084/150‑1 (MQ=255)

                                                                      ggCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCa                                                                               <  2:712106/150‑1 (MQ=255)

                                                                        cAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAgg                                                                             <  2:1034833/150‑1 (MQ=255)

                                                                          gCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGcc                                                                           >  1:888665/1‑150 (MQ=255)

                                                                            ggTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAg                                                                         >  2:1247827/1‑150 (MQ=255)

                                                                            ggTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAg                                                                         <  2:324156/150‑1 (MQ=255)

                                                                               tACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTa                                                                      <  1:1087012/150‑1 (MQ=255)

                                                                                 ccTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACa                                                                    >  1:1011452/1‑150 (MQ=255)

                                                                                 ccTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACa                                                                    >  2:794315/1‑150 (MQ=255)

                                                                                  cTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACag                                                                   >  1:1262669/1‑150 (MQ=255)

                                                                                  cTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGACTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACag                                                                   <  2:42140/150‑1 (MQ=255)

                                                                                   tGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACagc                                                                  >  2:416874/1‑149 (MQ=255)

                                                                                    ggCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACCTACACCCCCGCAGGCCAGTTACagag                                                                 >  2:1379331/1‑150 (MQ=255)

                                                                                      cGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGcc                                                               <  1:283971/150‑1 (MQ=255)

                                                                                        ggAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAg                                                             >  1:370765/1‑150 (MQ=255)

                                                                                         gAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCTGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGc                                                            >  1:279079/1‑150 (MQ=255)

                                                                                         gAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGc                                                            <  1:1447135/150‑1 (MQ=255)

                                                                                               aGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGa                                                      <  2:795512/150‑1 (MQ=255)

                                                                                               aGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGa                                                      >  1:1257637/1‑150 (MQ=255)

                                                                                                gCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGaa                                                     >  1:187500/1‑150 (MQ=255)

                                                                                                    ggCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGc                                                 >  1:1283414/1‑150 (MQ=255)

                                                                                                        ggACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTgg                                             <  1:1535187/150‑1 (MQ=255)

                                                                                                        ggACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTgg                                             >  1:1416547/1‑150 (MQ=255)

                                                                                                           cgccgcTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGtat                                          >  1:268001/1‑150 (MQ=255)

                                                                                                           cgccgcTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGtat                                          <  2:238393/150‑1 (MQ=255)

                                                                                                              cgcTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATg                                       <  2:77812/150‑1 (MQ=255)

                                                                                                                cagggcGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGAc                                     <  1:1563247/145‑1 (MQ=255)

                                                                                                                  ggTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCg                                   >  2:1430294/1‑150 (MQ=255)

                                                                                                                    tCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTg                                 <  2:1011452/150‑1 (MQ=255)

                                                                                                                      gAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGTAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGAc                               <  1:451484/150‑1 (MQ=255)

                                                                                                                      gAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGAc                               >  2:1419508/1‑150 (MQ=255)

                                                                                                                           taCGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACgg                          >  1:57071/1‑150 (MQ=255)

                                                                                                                            tcgcgGGACAGGCTGCACCGTGTGTCGGTGCGCAGGTTCGGCAGCATGGCAGGGAGTAATGCCGCATACGAACTGTCCAGGACATACACCCCCGCAGGCGAGTTACAGACCCAGCACCTTAACAGCCTTGTATATGACCGTGACTACggg                         <  1:1567852/149‑1 (MQ=255)

                                                                                                                            aCGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACggg                         >  2:1578110/1‑150 (MQ=255)

                                                                                                                            aCGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACggg                         >  1:208186/1‑150 (MQ=255)

                                                                                                                             cgcgGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGt                        <  1:1519305/150‑1 (MQ=255)

                                                                                                                                 ggACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAg                    >  1:43231/1‑150 (MQ=255)

                                                                                                                                  gACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGt                   >  2:1663106/1‑150 (MQ=255)

                                                                                                                                  gACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGt                   <  2:254888/150‑1 (MQ=255)

                                                                                                                                   aCAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTg                  <  1:1663106/150‑1 (MQ=255)

                                                                                                                                   aCAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTg                  <  2:674807/150‑1 (MQ=255)

                                                                                                                                     gctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGAc                >  1:1049806/6‑150 (MQ=255)

                                                                                                                                     gctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGAc                >  1:1089402/6‑150 (MQ=255)

                                                                                                                                     aGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGAc                >  1:1171925/1‑150 (MQ=255)

                                                                                                                                      ggCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACa               >  1:1243716/1‑150 (MQ=40)

                                                                                                                                      ggCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACa               >  1:525285/1‑150 (MQ=40)

                                                                                                                                       tgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACaa              >  1:861147/4‑150 (MQ=255)

                                                                                                                                          gCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACgg           >  1:283744/1‑150 (MQ=255)

                                                                                                                                               gttagaCGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCGAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGAcc      <  2:1162609/147‑1 (MQ=255)

                                                                                                                                                  agaCGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTgg   <  1:610764/150‑1 (MQ=39)

                                                                                                                                                   gaCGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGt  >  1:1011329/1‑150 (MQ=39)

ATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTGGT  >  NC_000913/525594‑525890

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