BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,122,286	T→G	42.9%	N108H (AAC→CAC)	pyrC ←	dihydro‑orotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,122,286	0	T	G	42.9%	96.4 / inf	224	N108H (AAC→CAC)	pyrC	dihydro‑orotase
Reads supporting (aligned to +/- strand): ref base T (70/58); new base G (41/55); total (111/113)
Fisher's exact test for biased strand distribution p-value = 8.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 229 total aligned reads displayed.

GCGCGCTTCACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAATACGCTGGCGATACG  >  NC_000913/1122140‑1122435

  gcgcTTCACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAGGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCg                                                                                                                                                  >  1:3659317/1‑150 (MQ=255)

  gcgcTTCACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCg                                                                                                                                                  >  1:3252771/1‑150 (MQ=255)

   cgcTTCACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCgg                                                                                                                                                 <  1:109477/150‑1 (MQ=255)

     cTTCACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGt                                                                                                                                               >  1:827289/1‑150 (MQ=255)

      ttCACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTa                                                                                                                                              >  2:269544/1‑150 (MQ=255)

      ttCACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTa                                                                                                                                              >  2:1742957/1‑150 (MQ=255)

      ttCACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTa                                                                                                                                              >  2:331862/1‑150 (MQ=255)

        cACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTaaa                                                                                                                                            <  1:1005350/150‑1 (MQ=255)

        cACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTaaa                                                                                                                                            <  1:1422080/150‑1 (MQ=255)

              cAAAAATGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTg                                                                                                                                      <  1:968637/150‑1 (MQ=255)

               aaaaaTGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGc                                                                                                                                     >  1:508991/1‑150 (MQ=255)

               aaaaaTGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGc                                                                                                                                     >  2:3503557/1‑150 (MQ=255)

               aaaaaTGTCGATATCTGCATGTGTCACTTCACCATGCACCAGAAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGACATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGGTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGc                                                                                                                                     >  2:2340906/1‑150 (MQ=255)

                aaaaTGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCa                                                                                                                                    <  1:984501/150‑1 (MQ=255)

                 aaaTGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAg                                                                                                                                   <  1:3939010/150‑1 (MQ=255)

                   aTGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCg                                                                                                                                 <  1:3567703/150‑1 (MQ=255)

                     gTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGt                                                                                                                               <  1:599991/150‑1 (MQ=255)

                       cGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGa                                                                                                                             >  1:3939418/1‑150 (MQ=255)

                         atatCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAac                                                                                                                           >  2:2034721/1‑150 (MQ=255)

                              tGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGcc                                                                                                                      >  2:3299528/1‑150 (MQ=255)

                                     gtCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTa                                                                                                               >  2:1321339/1‑150 (MQ=255)

                                      tCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTaa                                                                                                              >  1:890873/1‑150 (MQ=255)

                                        aCTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAAt                                                                                                            <  2:2393311/150‑1 (MQ=255)

                                          ttCACCATGCACCAGTAGCGGCATTCCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAATTcc                                                                                                          <  1:3120199/150‑1 (MQ=255)

                                          ttCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATcc                                                                                                          >  2:2525205/1‑150 (MQ=255)

                                           tCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCcg                                                                                                         <  1:1630726/150‑1 (MQ=255)

                                            cACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCcgc                                                                                                        >  1:2108965/1‑150 (MQ=255)

                                            cACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCcgc                                                                                                        >  1:894201/1‑150 (MQ=255)

                                             aCCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCcgcg                                                                                                       <  2:3252771/150‑1 (MQ=255)

                                                  gCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCa                                                                                                  <  1:913907/150‑1 (MQ=255)

                                                   cACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAg                                                                                                 <  1:3620884/150‑1 (MQ=255)

                                                      cAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTc                                                                                              >  1:309151/1‑150 (MQ=255)

                                                      cAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTc                                                                                              <  1:70966/150‑1 (MQ=255)

                                                        gTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCAt                                                                                            >  1:2996307/1‑150 (MQ=255)

                                                          aGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTa                                                                                          <  1:2844673/150‑1 (MQ=255)

                                                          aGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTa                                                                                          <  1:1701149/150‑1 (MQ=255)

                                                            cGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAgg                                                                                        >  2:1613935/1‑150 (MQ=255)

                                                               cATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATc                                                                                     >  2:165254/1‑150 (MQ=255)

                                                                   ccGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGc                                                                                 >  2:3373804/1‑150 (MQ=255)

                                                                    cGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCg                                                                                >  2:1241741/1‑150 (MQ=255)

                                                                             ccATGCGCTCAAGTACCGGCATGATTGCGTCAATGGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTa                                                                       <  1:3056152/150‑1 (MQ=255)

                                                                               aTGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTaaa                                                                     >  1:97601/1‑150 (MQ=255)

                                                                               aTGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTaaa                                                                     <  1:2414909/150‑1 (MQ=255)

                                                                                tGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAAt                                                                    <  1:1949023/150‑1 (MQ=255)

                                                                                 gcgcTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATa                                                                   >  2:2154553/1‑150 (MQ=255)

                                                                                   gcTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAAc                                                                 >  1:854434/1‑150 (MQ=255)

                                                                                   gcTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAAc                                                                 <  1:435392/150‑1 (MQ=255)

                                                                                    cTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACa                                                                >  2:2235170/1‑150 (MQ=255)

                                                                                    cTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACa                                                                <  1:1561235/150‑1 (MQ=255)

                                                                                    cTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACa                                                                <  1:1446361/150‑1 (MQ=255)

                                                                                     tCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAg                                                               <  1:3362131/150‑1 (MQ=255)

                                                                                         gTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGtca                                                           <  1:3299528/150‑1 (MQ=255)

                                                                                             cGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCaa                                                       >  1:1598413/1‑150 (MQ=255)

                                                                                             cGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCaa                                                       >  1:800042/1‑150 (MQ=255)

                                                                                              ggCATGATTGCGTCAAATGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTACGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAAt                                                      >  1:1028310/1‑150 (MQ=255)

                                                                                               gCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATg                                                     >  1:3538887/1‑150 (MQ=255)

                                                                                               gCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATg                                                     <  1:3774418/150‑1 (MQ=255)

                                                                                               gCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATg                                                     <  2:2787920/150‑1 (MQ=255)

                                                                                                cATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATgg                                                    <  2:2956341/150‑1 (MQ=255)

                                                                                                 aTGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATggg                                                   <  1:2766387/150‑1 (MQ=255)

                                                                                                 aTGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATggg                                                   <  1:2191119/150‑1 (MQ=255)

                                                                                                  tGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATgggg                                                  <  1:2333028/150‑1 (MQ=255)

                                                                                                   gATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGt                                                 <  1:1266813/150‑1 (MQ=255)

                                                                                                   gATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGt                                                 >  2:1269916/1‑150 (MQ=255)

                                                                                                    aTTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTg                                                >  2:1886729/1‑150 (MQ=255)

                                                                                                    aTTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTg                                                >  2:377133/1‑150 (MQ=255)

                                                                                                       gCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGaaa                                             >  1:2531382/1‑150 (MQ=255)

                                                                                                         gTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATc                                           >  1:1958368/1‑150 (MQ=255)

                                                                                                         gTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATc                                           >  1:1958425/1‑150 (MQ=255)

                                                                                                            aaTTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGTGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTg                                        <  1:1832957/150‑1 (MQ=255)

                                                                                                            aaTTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGTGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTg                                        <  1:2428565/150‑1 (MQ=255)

                                                                                                            aaTTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTg                                        <  1:1133065/150‑1 (MQ=255)

                                                                                                              ttGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGcc                                      <  1:3139656/150‑1 (MQ=255)

                                                                                                              ttGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAATTAACAGGTCATCAATGGGGTGAAATCGTGcc                                      >  2:2458187/1‑150 (MQ=255)

                                                                                                              ttGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGcc                                      >  2:3062908/1‑150 (MQ=255)

                                                                                                               tGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGccc                                     >  2:726896/1‑150 (MQ=255)

                                                                                                                  cGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGc                                  <  2:157795/150‑1 (MQ=255)

                                                                                                                   gTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCa                                 >  1:3881855/1‑150 (MQ=255)

                                                                                                                   gTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCggta                                 <  1:2495034/150‑3 (MQ=255)

                                                                                                                   gTCACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCa                                 <  1:1742957/150‑1 (MQ=255)

                                                                                                                    tCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAg                                >  2:1151156/1‑150 (MQ=255)

                                                                                                                    tCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAg                                >  2:1150987/1‑150 (MQ=255)

                                                                                                                     cACGCCGTGGCTGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAgg                               >  2:2862796/1‑150 (MQ=255)

                                                                                                                      aCGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGt                              >  2:3034162/1‑150 (MQ=255)

                                                                                                                          cGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACgg                          <  1:2268106/150‑1 (MQ=255)

                                                                                                                              gCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTc                      >  1:177702/1‑150 (MQ=255)

                                                                                                                               cTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCa                     >  2:3659536/1‑150 (MQ=255)

                                                                                                                                tGGAGTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCaa                    <  2:827257/150‑1 (MQ=255)

                                                                                                                                    gTTAGTGGTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAAt                <  1:2568059/150‑1 (MQ=255)

                                                                                                                                     ttAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAATa               <  2:2421039/150‑1 (MQ=255)

                                                                                                                                       aGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAATACg             >  2:1783458/1‑150 (MQ=255)

                                                                                                                                           gTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAATACGCTgg         <  1:2964323/150‑1 (MQ=255)

                                                                                                                                           gTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAATACGCTgg         <  2:3047297/150‑1 (MQ=255)

                                                                                                                                           gTTGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAATACGCTgg         >  1:2608050/1‑150 (MQ=255)

                                                                                                                                            ttGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAATACg             >  2:553040/1‑145 (MQ=255)

                                                                                                                                            ttGCGTGTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAATACg             <  1:553040/145‑1 (MQ=255)

                                                                                                                                                 tttGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAATACGCTGGCGATAc   <  2:139831/150‑1 (MQ=255)

                                                                                                                                                 tttGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAATACGCTGGCGATAc   >  1:3400177/1‑150 (MQ=255)

                                                                                                                                                  ttGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAATACGCTGGCGATACg  >  1:2799399/1‑150 (MQ=255)

                                                                                                                                                  ttGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAATACGCTGGCGATACg  >  1:2799500/1‑150 (MQ=255)

GCGCGCTTCACGATCAAAAATGTCGATATCTGCATGTGTCACTTCACCATGCACCAGTAGCGGCATACCGATTTTTTCCATGCGCTCAAGTACCGGCATGATTGCGTCAATTGACGTCACGCCGTGGCTGGAGTTAGTGGTTGCGTTTGCCGGGTAAAGTTTTGCAGCGGTGAACACGCCTTCGTTAAATCCGCGCTCCAGCTCATTAGGATCCAGCGAATCTGTTAAATAACAGGTCATCAATGGGGTGAAATCGTGCCCGGCAGGTACGGCGTCAAGAATACGCTGGCGATACG  >  NC_000913/1122140‑1122435

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