BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	325,767	A→G	V494A (GTA→GCA)	betA ←	choline dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	325,767	0	A	G	79.3%	598.3 / 150.6	251	V494A (GTA→GCA)	betA	choline dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (28/24); new base G (105/94); total (133/118)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 256 total aligned reads displayed.

TTTTCTCGCCAATCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCA  >  NC_000913/325669‑325865

ttttCTCGCCAATCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGcg                                                                                                  >  1:434962/1‑101 (MQ=255)

  ttCTCGCCAATCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGcgcg                                                                                                >  2:434966/1‑101 (MQ=255)

   tctcGCCAATCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCgg                                                                                               <  1:434597/101‑1 (MQ=255)

   tctcGCCAATCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCgg                                                                                               >  1:434967/1‑101 (MQ=255)

   tctcGCCAATCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCgg                                                                                               >  1:434970/1‑101 (MQ=255)

    ctcGCCAATCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGc                                                                                              >  1:434973/1‑101 (MQ=255)

        ccAATCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTc                                                                                          >  2:434977/1‑101 (MQ=255)

        ccAATCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTc                                                                                          <  1:434588/101‑1 (MQ=255)

         cAATCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCg                                                                                         <  1:434724/101‑1 (MQ=255)

          aaTCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGc                                                                                        >  2:434981/1‑101 (MQ=255)

          aaTCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGc                                                                                        <  1:434570/101‑1 (MQ=255)

          aaTCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGc                                                                                        >  1:434980/1‑101 (MQ=255)

           aTCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGcc                                                                                       <  2:434644/101‑1 (MQ=255)

             cATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGt                                                                                     <  2:434449/101‑1 (MQ=255)

              aTAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTc                                                                                    <  2:434238/101‑1 (MQ=255)

              aTAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTc                                                                                    >  2:434987/1‑101 (MQ=255)

              aTAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTc                                                                                    >  1:434986/1‑101 (MQ=255)

                aaTTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCaa                                                                                  <  1:434611/101‑1 (MQ=255)

                aaTTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCa                                                                                   >  2:434990/1‑100 (MQ=255)

                  ttGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAAcc                                                                                >  1:434994/1‑101 (MQ=255)

                     tCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACg                                                                             >  1:434996/1‑101 (MQ=255)

                      cGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACgg                                                                            >  2:434997/1‑101 (MQ=255)

                       gTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCa                                                                                   <  2:434582/93‑1 (MQ=255)

                           cGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACAtc                                                                       <  2:434624/101‑1 (MQ=255)

                           cGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACAtc                                                                       >  2:435001/1‑101 (MQ=255)

                            gTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACAtct                                                                      >  2:435003/1‑101 (MQ=255)

                             ttCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACg                                                                             <  2:434551/93‑1 (MQ=255)

                             ttCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACAtctc                                                                     >  1:435004/1‑101 (MQ=255)

                                aaaTTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCtcgtc                                                                  >  1:435010/1‑101 (MQ=255)

                                aaaTTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCtcgtc                                                                  <  2:434647/101‑1 (MQ=255)

                                aaaTTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCtcgtc                                                                  >  2:435008/1‑101 (MQ=255)

                                  aTTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTcgtcgt                                                                >  2:435012/1‑101 (MQ=255)

                                   ttCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTa                                                               <  1:434564/101‑1 (MQ=255)

                                      ccGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAAcc                                                            >  1:435015/1‑101 (MQ=255)

                                        ggTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCa                                                          >  1:435019/1‑101 (MQ=255)

                                        ggTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCa                                                          <  1:434771/101‑1 (MQ=255)

                                        ggTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCa                                                          <  2:434612/101‑1 (MQ=255)

                                         gTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCAt                                                         >  2:435021/1‑101 (MQ=255)

                                           gATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCAttt                                                       >  1:435022/1‑101 (MQ=255)

                                             tAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTg                                                     <  2:434513/101‑1 (MQ=255)

                                             tAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTg                                                     <  2:434749/101‑1 (MQ=255)

                                              aaTCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGc                                                    >  2:435025/1‑101 (MQ=255)

                                               aTCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTGGTAACCCATTTTGCa                                                   <  2:434415/101‑1 (MQ=255)

                                                  tGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGt                                                >  1:435029/1‑101 (MQ=255)

                                                   gCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTa                                               >  2:435031/1‑101 (MQ=255)

                                                   gCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTa                                               >  1:435032/1‑101 (MQ=255)

                                                    cGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTAc                                              >  1:435037/1‑101 (MQ=255)

                                                    cGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTAc                                              >  1:435035/1‑101 (MQ=255)

                                                     ggCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTAcc                                             <  2:434766/101‑1 (MQ=255)

                                                      gCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCg                                            <  2:434740/101‑1 (MQ=255)

                                                       cATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGc                                           >  1:435041/1‑101 (MQ=255)

                                                       cATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGc                                           <  1:434747/101‑1 (MQ=255)

                                                        aTAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCa                                          >  1:435043/1‑101 (MQ=255)

                                                        aTAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCa                                          >  2:435046/1‑101 (MQ=255)

                                                        aTAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCa                                          <  1:434815/101‑1 (MQ=255)

                                                         tAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCAc                                         >  1:435048/1‑101 (MQ=255)

                                                          aaTCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACg                                        <  2:434361/101‑1 (MQ=255)

                                                          aaTCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACg                                        >  2:435049/1‑101 (MQ=255)

                                                           aTCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACgg                                       >  2:435053/1‑101 (MQ=255)

                                                           aTCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCa                                                                                   <  1:434581/57‑1 (MQ=255)

                                                            tCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGa                                      <  2:434620/101‑1 (MQ=255)

                                                            tCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGa                                      <  2:434569/101‑1 (MQ=255)

                                                              gACGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATg                                    >  1:435055/1‑101 (MQ=255)

                                                                cGCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGa                                  >  1:435058/1‑101 (MQ=255)

                                                                cGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGa                                  >  2:435059/1‑101 (MQ=255)

                                                                 gCATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGaa                                 >  1:435063/1‑101 (MQ=255)

                                                                 gCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGaa                                 >  2:435065/1‑101 (MQ=255)

                                                                  cATCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAg                                >  1:435067/1‑101 (MQ=255)

                                                                   aTCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAgg                               >  1:435070/1‑101 (MQ=255)

                                                                   aTCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAgg                               <  1:434670/101‑1 (MQ=255)

                                                                    tCCACCACACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGc                              >  2:435072/1‑101 (MQ=255)

                                                                       accacACGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGt                           <  1:434693/101‑1 (MQ=255)

                                                                       accacACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGt                           <  1:434468/101‑1 (MQ=255)

                                                                       accacACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGt                           <  2:434720/101‑1 (MQ=255)

                                                                          acacGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTc                        >  1:435075/1‑101 (MQ=255)

                                                                          acacGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTc                        <  1:434327/101‑1 (MQ=255)

                                                                           cacGCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCg                       <  1:434508/101‑1 (MQ=255)

                                                                              gCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCg                    <  2:434696/101‑1 (MQ=255)

                                                                              gCAGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCg                    >  1:435079/1‑101 (MQ=255)

                                                                              gCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCg                    <  1:434608/101‑1 (MQ=255)

                                                                               cAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGt                   <  1:434796/101‑1 (MQ=255)

                                                                               cAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGt                   >  2:435082/1‑101 (MQ=255)

                                                                                aGGCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTg                  >  1:435083/1‑101 (MQ=255)

                                                                                 ggCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTgg                 <  2:434814/101‑1 (MQ=255)

                                                                                  gCCTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGt                >  1:435085/1‑101 (MQ=255)

                                                                                   ccTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGtt               <  1:434638/101‑1 (MQ=255)

                                                                                   ccTTCTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGtt               <  1:434497/101‑1 (MQ=255)

                                                                                     ttCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTAc             <  2:434784/101‑1 (MQ=255)

                                                                                       cTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCg                       <  2:434672/89‑1 (MQ=255)

                                                                                       cTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGc           <  1:434702/101‑1 (MQ=255)

                                                                                       cTAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGc           <  1:434667/101‑1 (MQ=255)

                                                                                        tAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCa          >  1:435089/1‑101 (MQ=255)

                                                                                        tAACCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCa          >  2:435090/1‑101 (MQ=255)

                                                                                         aaCCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAAcc                                                            <  2:434725/50‑1 (MQ=255)

                                                                                          aCCCGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACg        <  1:434622/101‑1 (MQ=255)

                                                                                           cccGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGa       <  2:434874/101‑1 (MQ=255)

                                                                                            ccGTGTGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGaa      >  2:435094/1‑101 (MQ=255)

                                                                                            ccGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGaa      >  1:435095/1‑101 (MQ=255)

                                                                                              gtgtGCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACt    <  1:434891/101‑1 (MQ=255)

                                                                                               tgtgCGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTc   <  2:434687/101‑1 (MQ=255)

TTTTCTCGCCAATCATAATTGTCGTGGCGTTCAAATTCCCGGTGATAATCTGCGGCATAATCGACGCATCCACCACACGCAGGCCTTCTAACCCGTGTACGCGGCCTTCGCCGTCAACCACGGACATCTCGTCGTAACCCATTTTGCAGGTACCGCACGGATGGAAGGCGGTTTCGGCGTGGTTACGCACGAACTCA  >  NC_000913/325669‑325865

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