BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	903,254	A→G	V146A (GTC→GCC)	ybjP ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	903,254	0	A	G	97.7%	397.3 / ‑0.7	131	V146A (GTC→GCC)	ybjP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (53/75); total (54/77)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

Only 100 of 137 total aligned reads displayed.

TGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTACGGGCA  >  NC_000913/903164‑903346

tGACCAAGGCGTTAGCGTTTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGccc                                                                                           <  2:3652373/93‑1 (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:1877251/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)

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

         cGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCCACCCTCGc                                                                                  >  2:3201536/1‑93 (MQ=255)

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

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

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

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

           ttAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCt                                                                                >  1:603719/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:4105837/1‑93 (MQ=255)

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

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

               cGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCt                                                                            <  1:3080177/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                                                                         >  1:3315503/1‑93 (MQ=255)

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

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

                       gATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCAtt                                                                    >  1:3029837/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                                                            <  1:662202/93‑1 (MQ=255)

                                 cGGAGCGTCCCCGCCGTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCAt                                                          <  2:381392/93‑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‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCtc                                                      >  2:3448546/1‑92 (MQ=255)

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

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

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

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

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

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

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

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

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

                                             gccgTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCt                                              <  1:987110/93‑1 (MQ=255)

                                              ccgTGGCATGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTg                                             <  2:566212/93‑1 (MQ=255)

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

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

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

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

                                                      tGAACGCTACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCt                                     >  2:4039217/1‑93 (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)

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

                                                             tACCGCC‑AAGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTgag                              >  2:2987685/1‑93 (MQ=255)

                                                                    caaGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCAc                        <  2:1972016/92‑1 (MQ=255)

                                                                    caaGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCAc                        <  1:735867/92‑1 (MQ=255)

                                                                    caaGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCAc                        <  1:4204559/92‑1 (MQ=255)

                                                                    caaGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAAAACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCAc                        <  2:2937450/92‑1 (MQ=255)

                                                                     aaGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCAcg                       >  1:2697020/1‑93 (MQ=255)

                                                                      aGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCAcgc                      <  2:1501080/93‑1 (MQ=255)

                                                                      aGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCAcgc                      <  2:440821/93‑1 (MQ=255)

                                                                      aGTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCAcgc                      <  2:1501081/93‑1 (MQ=255)

                                                                       gTAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCAcgcg                     >  1:1943918/1‑93 (MQ=255)

                                                                        tAACGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACgcgc                    <  1:4425003/93‑1 (MQ=255)

                                                                         aaCGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCa                   <  1:867402/93‑1 (MQ=255)

                                                                         aaCGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCa                   <  1:4176746/93‑1 (MQ=255)

                                                                          aCGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCaa                  <  2:4515261/93‑1 (MQ=255)

                                                                          aCGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCaa                  >  1:4443547/1‑93 (MQ=255)

                                                                           cGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAAt                 <  2:3320467/93‑1 (MQ=255)

                                                                           cGCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAAt                 <  1:4727750/93‑1 (MQ=255)

                                                                            gCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATg                >  1:859989/1‑93 (MQ=255)

                                                                            gCACATCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATg                <  2:4042967/93‑1 (MQ=255)

                                                                            gCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATg                <  2:4524682/93‑1 (MQ=255)

                                                                            ccacaTCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATg                <  1:2390520/92‑1 (MQ=255)

                                                                               caTCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGaa             >  2:201754/1‑93 (MQ=255)

                                                                                aTCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAAt            >  1:648639/1‑93 (MQ=255)

                                                                                aTCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAAt            >  1:222406/1‑93 (MQ=255)

                                                                                aTCGTCGATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAAt            <  2:2968855/93‑1 (MQ=255)

                                                                                 tcgtcgATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCACCCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATg           <  1:2359281/93‑1 (MQ=255)

                                                                                 tcgtcgATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATg           >  2:3336289/1‑93 (MQ=255)

                                                                                   gtcgATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGtt         <  1:391796/93‑1 (MQ=255)

                                                                                   gtcgATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGtt         <  2:4332252/93‑1 (MQ=255)

                                                                                    tcgATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATACTGCCAACCCTGGTCACTCTGTTTGAGATCCACGCGCAATGGAATGTTa        <  2:2218850/93‑1 (MQ=255)

                                                                                     cgATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTAc       <  2:2696724/93‑1 (MQ=255)

                                                                                     cgATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTAc       <  1:1689534/93‑1 (MQ=255)

                                                                                      gATGGCCCAACGCTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTACg      <  2:91172/93‑1 (MQ=255)

                                                                                      gATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTACg      >  2:3674233/1‑93 (MQ=255)

                                                                                      gATGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTACg      <  1:1863416/93‑1 (MQ=255)

                                                                                       aTGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTACgg     <  2:4733720/93‑1 (MQ=255)

                                                                                       aTGGCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTACgg     <  1:1270625/93‑1 (MQ=255)

                                                                                          aCCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTACGGGCa  >  2:744995/1‑93 (MQ=255)

TGATCAAGGCGTTAGCGATTCTCGATGGACTGACGGAGCGTACCCGCCGTGGCATGAACGCTACCGCCTAAGTAACGCACATCGTCGATGACCCAACACTGGCCTTCCTGAATCATTAACACTTCATCCTGCCAACCCTGGTCACCCTGTTTGAGATCCACGCGCAATGGAATGTTACGGGCA  >  NC_000913/903164‑903346

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