BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	477,096	A→G	L251P (CTC→CCC)	ylaB ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	477,096	0	A	G	98.0%	624.5 / 1.3	199	L251P (CTC→CCC)	ylaB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (3/1); new base G (109/86); total (112/87)
Fisher's exact test for biased strand distribution p-value = 6.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

Only 100 of 207 total aligned reads displayed.

ATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGT  >  NC_000913/477006‑477185

aTAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGc                                                                                         >  2:1932511/1‑93 (MQ=255)

 tAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCg                                                                                        >  1:1013067/1‑93 (MQ=255)

  aGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCTCAAAATACGCAGCACAAACAGCGCTGCCAGCAGGCCAATCACCAACCCGGGCgg                                                                                       <  2:3431871/93‑1 (MQ=255)

  aGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCgg                                                                                       <  2:3134841/93‑1 (MQ=255)

  aGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCgg                                                                                       >  2:2163063/1‑93 (MQ=255)

   gCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGt                                                                                      >  2:2800077/1‑93 (MQ=255)

     aTCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTaa                                                                                    >  2:1991979/1‑93 (MQ=255)

      tCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAc                                                                                   >  2:3045160/1‑93 (MQ=255)

      tCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAc                                                                                   >  1:1441446/1‑93 (MQ=255)

       ccTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAcc                                                                                  <  2:348569/93‑1 (MQ=255)

       ccTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAAcc                                                                                  >  1:2560567/1‑93 (MQ=255)

        cTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCa                                                                                 >  1:1554745/1‑93 (MQ=255)

        cTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCa                                                                                 >  2:3857351/1‑93 (MQ=255)

        cTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCa                                                                                 >  1:4568773/1‑93 (MQ=255)

         tGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCaa                                                                                >  1:2625066/1‑93 (MQ=255)

         tGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGACCAACCACCAACCCGGGCGGCAACCaa                                                                                >  2:2502946/1‑93 (MQ=255)

          gCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCaaa                                                                               >  1:4393371/1‑93 (MQ=255)

          gCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCaaa                                                                               <  2:3109515/93‑1 (MQ=255)

             gCCGATGATGCGGTGACTGAATACGGCGCAAAATTCGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATaa                                                                            <  2:3940179/93‑1 (MQ=255)

             gCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATaa                                                                            >  2:3464847/1‑93 (MQ=255)

              ccGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATaaa                                                                           >  1:1194164/1‑93 (MQ=255)

               cGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAg                                                                          >  2:2746255/1‑93 (MQ=255)

               cGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAg                                                                          >  2:1230353/1‑93 (MQ=255)

                gatgatGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGa                                                                         >  1:4148813/1‑93 (MQ=255)

                 atgatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGAc                                                                        <  2:4316140/93‑1 (MQ=255)

                  tgatgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACt                                                                       >  1:1255886/1‑93 (MQ=255)

                     tgCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGc                                                                    >  1:3709306/1‑93 (MQ=255)

                        ggTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCCATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGa                                                                 <  1:4016580/93‑1 (MQ=255)

                        ggTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGa                                                                 <  1:2762591/93‑1 (MQ=255)

                        ggTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGa                                                                 <  1:3150353/93‑1 (MQ=255)

                        ggTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGa                                                                 <  2:1375154/93‑1 (MQ=255)

                        ggTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGa                                                                 <  1:1962714/93‑1 (MQ=255)

                        ggTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGa                                                                 <  1:4003965/93‑1 (MQ=255)

                         gTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGAt                                                                <  2:303690/93‑1 (MQ=255)

                           gACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGc                                                              >  2:799256/1‑93 (MQ=255)

                           gACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGc                                                              <  1:4584436/93‑1 (MQ=255)

                            aCTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGcc                                                             <  2:1071756/93‑1 (MQ=255)

                            aCTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGcc                                                             <  1:2116548/93‑1 (MQ=255)

                              tGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAg                                                           >  1:2824136/1‑93 (MQ=255)

                               gAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGc                                                          >  2:3245832/1‑93 (MQ=255)

                               gAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGc                                                          >  2:2230889/1‑93 (MQ=255)

                               gAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGc                                                          >  1:2504803/1‑93 (MQ=255)

                                 aTACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCt                                                        >  1:424619/1‑93 (MQ=255)

                                   aCGGCGCAAAATACGCAGCACAAACATCGCTG‑CAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGGCGATGCCAGCCTTTc                                                     >  2:1740804/1‑93 (MQ=255)

                                    cGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTc                                                     <  2:2654797/93‑1 (MQ=255)

                                      gcgcAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTg                                                   <  2:4372506/93‑1 (MQ=255)

                                        gcAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCa                                                 >  1:3416828/1‑93 (MQ=255)

                                        gcAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCa                                                 <  1:1046445/93‑1 (MQ=255)

                                          aaaaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCATCCTTTCTGCAAc                                               <  2:2365713/93‑1 (MQ=255)

                                            aaTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAt                                             >  2:1296984/1‑93 (MQ=255)

                                             aTACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAtt                                            >  1:860456/1‑93 (MQ=255)

                                              tACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACAttt                                           >  1:4598491/1‑93 (MQ=255)

                                              tACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCAATGCCAGCCTTTCTGCAACAttt                                           >  2:2588646/1‑93 (MQ=255)

                                                  cagcaCAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGt                                       <  2:2937707/93‑1 (MQ=255)

                                                  cagcaCAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGt                                       <  1:3981387/93‑1 (MQ=255)

                                                    gcaTAAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTg                                     >  1:56602/1‑93 (MQ=255)

                                                     cacaAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGa                                    >  2:1140075/1‑93 (MQ=255)

                                                     cacaAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGa                                    <  1:1184781/93‑1 (MQ=255)

                                                      acaAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGaa                                   <  2:3286954/93‑1 (MQ=255)

                                                      acaAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGaa                                   >  2:1641309/1‑93 (MQ=255)

                                                      acaAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGaa                                   >  2:4406959/1‑93 (MQ=255)

                                                      acaAACATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGaa                                   <  1:4476856/93‑1 (MQ=255)

                                                        aaaCATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGc                                 >  2:3605587/1‑93 (MQ=255)

                                                         aaCATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGcc                                >  2:507053/1‑93 (MQ=255)

                                                         aaCATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGcc                                <  2:4041259/93‑1 (MQ=255)

                                                           cATCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCa                              >  1:2002362/1‑93 (MQ=255)

                                                            aTCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAt                             >  2:4500089/1‑93 (MQ=255)

                                                             tCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtg                            <  1:2853376/93‑1 (MQ=255)

                                                             tCGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtg                            >  2:3936975/1‑93 (MQ=255)

                                                              cGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtgt                           <  1:141953/93‑1 (MQ=255)

                                                              cGCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtgt                           <  1:3040310/93‑1 (MQ=255)

                                                               gCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtgtg                          >  2:3681885/1‑93 (MQ=255)

                                                               gCTGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCAtgtg                          >  1:2917096/1‑93 (MQ=255)

                                                                 tGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgat                        <  1:436676/93‑1 (MQ=255)

                                                                 tGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgat                        >  2:2135657/1‑93 (MQ=255)

                                                                 tGCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgat                        <  1:3651679/93‑1 (MQ=255)

                                                                  gCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgatg                       >  1:3702066/1‑93 (MQ=255)

                                                                  gCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgatg                       <  1:3259281/93‑1 (MQ=255)

                                                                  gCCAGCAGGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGtgatg                       >  1:1659184/1‑93 (MQ=255)

                                                                    cagcagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTgatgat                     >  1:1511217/1‑93 (MQ=255)

                                                                     agcagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATc                    >  2:4151723/1‑93 (MQ=255)

                                                                     agcagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATc                    <  2:4620327/93‑1 (MQ=255)

                                                                       cagGCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGa                  <  1:3226233/93‑1 (MQ=255)

                                                                         ggCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGaaa                >  1:752094/1‑93 (MQ=255)

                                                                          gCCAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAat               >  2:778522/1‑93 (MQ=255)

                                                                            cAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAatat             >  1:2275065/1‑93 (MQ=255)

                                                                            cAATCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAatat             >  1:3879599/1‑93 (MQ=255)

                                                                             aaTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATAtt            <  1:1716713/93‑1 (MQ=255)

                                                                             aaTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATAtt            >  2:4177159/1‑93 (MQ=255)

                                                                              aTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTc           <  1:917347/93‑1 (MQ=255)

                                                                              aTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATGTTCGTTGAAGCCCATGTGATGATCGAAATATTc           <  1:1045954/93‑1 (MQ=255)

                                                                              aTCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCAAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGTAATATTc           >  1:4689534/1‑93 (MQ=255)

                                                                               tCACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCCAAATATTCa          <  2:1329908/93‑1 (MQ=255)

                                                                                cACCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAt         <  1:2213179/93‑1 (MQ=255)

                                                                                 aCCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtc        <  1:2377739/93‑1 (MQ=255)

                                                                                 aCCAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtc        <  2:1713137/93‑1 (MQ=255)

                                                                                  ccAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtct       >  1:4548077/1‑93 (MQ=255)

                                                                                   cAACCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCAtctc      >  1:4412210/1‑93 (MQ=255)

                                                                                    aaCCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTcc     >  2:3199296/1‑93 (MQ=255)

                                                                                    aaCCCGGGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTcc     <  2:4486632/93‑1 (MQ=255)

ATAGCATCCTGCAGCCGATGATGCGGTGACTGAATACGGCGCAAAATACGCAGCACAAACATCGCTGCCAGCAGGCCAATCACCAACCCGAGCGGTAACCAAATAAAGACTTGCCGATGCCAGCCTTTCTGCAACATTTTCGTTGAAGCCCATGTGATGATCGAAATATTCATCTCCGGT  >  NC_000913/477006‑477185

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