BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	903,232	Δ1 bp	coding (459/516 nt)	ybjP ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	903,232	0	T	.	98.4%	566.3 / 2.9	128	coding (459/516 nt)	ybjP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base . (59/66); total (60/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 134 total aligned reads displayed.

AATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGA  >  NC_000913/903144‑903319

aaTTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTa                                                                                    <  2:2173868/93‑1 (MQ=255)

 aTTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCC‑AAGTaa                                                                                   <  2:4130600/93‑1 (MQ=255)

 aTTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTaa                                                                                   >  2:2907860/1‑93 (MQ=255)

 aTTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTaa                                                                                   <  2:1713465/93‑1 (MQ=255)

  tttACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAAc                                                                                  <  2:3314794/93‑1 (MQ=255)

  tttACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAAc                                                                                  >  2:1671150/1‑93 (MQ=255)

  tttACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAAc                                                                                  <  2:3805853/93‑1 (MQ=255)

   ttACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACg                                                                                 >  1:3529669/1‑93 (MQ=255)

    tagggggTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGc                                                                                <  1:4491880/90‑1 (MQ=255)

    tACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGc                                                                                <  1:82540/93‑1 (MQ=255)

      cGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGcac                                                                              >  1:3630339/1‑93 (MQ=255)

       ggggTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGcaca                                                                             >  1:2235109/1‑93 (MQ=255)

        gggTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACAt                                                                            >  1:3617043/1‑93 (MQ=255)

         gttttATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACAtc                                                                           <  1:3004533/91‑1 (MQ=255)

         ggTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGCCTGACGGAGCTTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACCCACAtc                                                                           <  1:3004532/93‑1 (MQ=255)

         ggTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACAtc                                                                           <  2:2465051/93‑1 (MQ=255)

          gTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACAtcg                                                                          >  1:4498356/1‑93 (MQ=255)

           tttATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACAtcgt                                                                         >  2:1766635/1‑93 (MQ=255)

             tATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACAtcgtcg                                                                       >  2:3789293/1‑93 (MQ=255)

              aTTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGa                                                                      <  1:2218535/93‑1 (MQ=255)

              aTTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGa                                                                      <  2:4067459/93‑1 (MQ=255)

              aTTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGa                                                                      <  1:3052841/93‑1 (MQ=255)

               ttGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGAt                                                                     >  2:1769668/1‑93 (MQ=255)

                tGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATg                                                                    >  2:3789677/1‑93 (MQ=255)

                tGGTTGATCAAGGAGTTAGCGATTCTCGATGGACTGACGGAGCGGACACGCCGAGGCATGAACGCTACCTCCTAAGCAACGCACATCGTCGAt                                                                     >  1:4186703/1‑93 (MQ=255)

                tGGTTGAGCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATg                                                                    <  1:1050253/93‑1 (MQ=255)

                tGGTTGAACAGGGCGTTAGCGATTCTCGATGGACTGAGGAAGCGTTCCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATg                                                                    <  1:3229520/93‑1 (MQ=255)

                 ggTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATgg                                                                   >  1:2793590/1‑92 (MQ=255)

                 gattGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATgg                                                                   <  1:554774/91‑2 (MQ=255)

                     gATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCa                                                               >  2:3192348/1‑93 (MQ=255)

                       tCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAac                                                             <  2:1435941/93‑1 (MQ=255)

                       tCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAac                                                             >  2:1077639/1‑93 (MQ=255)

                       tCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAac                                                             >  1:2924204/1‑93 (MQ=255)

                         agggCGTTAGCGATTCTCGAGGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAacac                                                           <  2:334514/91‑1 (MQ=255)

                           ggCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTg                                                         >  1:1922961/1‑93 (MQ=255)

                              gTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGcc                                                      <  1:4148272/93‑1 (MQ=255)

                              gTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGcc                                                      >  2:2381942/1‑93 (MQ=255)

                               ttAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCt                                                     >  1:603719/1‑93 (MQ=255)

                               ttAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCt                                                     >  1:1773573/1‑93 (MQ=255)

                               ttAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCt                                                     >  2:4115800/1‑93 (MQ=255)

                               ttAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCt                                                     >  1:1401815/1‑93 (MQ=255)

                                tAGTGATTCTCGCTGGACTGACGGAGCGTCACCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCt                                                     <  1:3691836/93‑1 (MQ=255)

                                 aGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTc                                                   <  1:1422239/93‑1 (MQ=255)

                                  gCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTcc                                                  <  2:4692693/93‑1 (MQ=255)

                                  gCGATCCGCGATGGACTGCCGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTcc                                                  <  1:2304393/93‑1 (MQ=255)

                                     aTTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGa                                               >  2:909061/1‑93 (MQ=255)

                                      ttCTCGATGGACTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGaa                                              <  2:3485709/93‑1 (MQ=255)

                                      ttCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAAGGCACATCGTCGCTGGCCCAACACTGGCCTTCCTGaa                                              >  1:3315502/1‑93 (MQ=255)

                                      ttCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGaa                                              >  2:2063470/1‑93 (MQ=255)

                                      ttCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGaa                                              >  2:332037/1‑93 (MQ=255)

                                      ttCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGaa                                              >  1:3315503/1‑93 (MQ=255)

                                         tcGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCa                                           <  1:2614615/93‑1 (MQ=255)

                                         ccGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCa                                           <  1:2511556/92‑1 (MQ=255)

                                           gATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCAtt                                         >  1:3029836/1‑93 (MQ=255)

                                             tGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTaa                                       >  1:3051604/1‑93 (MQ=255)

                                               gACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAaca                                     >  2:2572056/1‑93 (MQ=255)

                                                aCTGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAacac                                    <  2:3506781/93‑1 (MQ=255)

                                                  tGACGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACtt                                  <  1:4338975/93‑1 (MQ=255)

                                                   gACGGAGCGTCCCCGCCGTGGCATGAACGCTACCG‑CAAAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTc                                 <  1:662201/93‑1 (MQ=255)

                                                   gACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTc                                 >  2:3867926/1‑93 (MQ=255)

                                                    agggAGCGTTCCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCa                                <  2:3569470/91‑1 (MQ=255)

                                                      ggAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATc                              <  2:414555/93‑1 (MQ=255)

                                                      ggAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATAATTAACACTTCATc                              >  2:222978/1‑93 (MQ=255)

                                                        ggCGTCCCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCt                            <  1:1012749/92‑1 (MQ=255)

                                                         gCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTg                           >  2:2175641/1‑93 (MQ=255)

                                                         gCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGACCAACACTGGACTTCCTGAATCATTAACACTTCATACTg                           >  2:456085/1‑93 (MQ=255)

                                                          cGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGc                          <  1:4726595/93‑1 (MQ=255)

                                                           gtccccGCGGTGGCATGAACGCTACCGCC‑AAGTAACGCACGTCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGcc                         <  1:1080618/90‑1 (MQ=255)

                                                           gtccccGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGcc                         <  1:659430/90‑1 (MQ=255)

                                                           gTACCCGCCGTGGCATGAACGCTACCGCC‑ACGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGcc                         <  1:4340016/93‑1 (MQ=255)

                                                           gTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTAATCCTGcc                         >  1:3385175/1‑93 (MQ=255)

                                                            tccccGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCa                        <  2:430622/91‑1 (MQ=255)

                                                            tACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCa                        >  2:910188/1‑93 (MQ=255)

                                                              cccGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAAc                      >  2:1910542/1‑93 (MQ=255)

                                                              cccGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAAc                      <  1:2624495/93‑1 (MQ=255)

                                                              cccGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAAc                      >  1:3076647/1‑93 (MQ=255)

                                                               ccgccgTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAAcc                     <  1:729601/93‑1 (MQ=255)

                                                                cgccgTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAAccc                    <  2:3107384/93‑1 (MQ=255)

                                                                cgccgTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAAccc                    <  2:4366388/93‑1 (MQ=255)

                                                                cgccgTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAAccc                    <  2:1987496/93‑1 (MQ=255)

                                                                 gccgTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCt                   >  1:3509858/1‑93 (MQ=255)

                                                                    gTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGt                >  1:996719/1‑93 (MQ=255)

                                                                     tGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTc               <  2:3139525/93‑1 (MQ=255)

                                                                     tGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTc               >  1:2960982/1‑93 (MQ=255)

                                                                      cgCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCa              <  1:2104536/92‑1 (MQ=255)

                                                                       gCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACGCTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCAc             >  2:1228924/1‑93 (MQ=255)

                                                                       gCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCAc             >  1:4188096/1‑93 (MQ=255)

                                                                        cATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCAc             >  2:2050672/1‑93 (MQ=255)

                                                                        cATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCAcc            >  1:1927369/1‑93 (MQ=255)

                                                                          tGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCt          >  2:4039217/1‑93 (MQ=255)

                                                                          tGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCt          >  1:2230606/1‑93 (MQ=255)

                                                                          agAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCt          <  1:943533/92‑1 (MQ=255)

                                                                            aaCGCTACCGCC‑AAGTACCGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGt        <  2:2824020/93‑1 (MQ=255)

                                                                            aaCGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGGCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGt        <  1:1977361/93‑1 (MQ=255)

                                                                            aaCGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGt        <  1:2117659/93‑1 (MQ=255)

                                                                             aCGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGtt       <  1:4469861/93‑1 (MQ=255)

                                                                             aCGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGtt       <  2:4744016/93‑1 (MQ=255)

                                                                               gCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTg     >  2:2952399/1‑93 (MQ=255)

                                                                                cTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTga    <  2:464964/93‑1 (MQ=255)

                                                                                  aCCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTgaga  <  1:530907/93‑1 (MQ=255)

AATTTACGGGGTTTATTGGTTGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGA  >  NC_000913/903144‑903319

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