BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,665,031	A→G	V279A (GTG→GCG)	hcaT ←	HcaT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,665,031	0	A	G	99.6%	833.2 / NA	269	V279A (GTG→GCG)	hcaT	HcaT
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (115/153); minor base T (0/1); total (115/154)
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 284 total aligned reads displayed.

AGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCC  >  NC_000913/2664941‑2665120

aGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCt                                                                                         >  1:7505691/1‑93 (MQ=255)

 ggTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTa                                                                                        >  2:1419934/1‑93 (MQ=255)

  gTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTAc                                                                                       <  1:1194795/93‑1 (MQ=255)

  gTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTAc                                                                                       <  2:1440770/93‑1 (MQ=255)

  gTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTAc                                                                                       >  1:3620159/1‑93 (MQ=255)

    ggCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTAcgc                                                                                     <  1:4619804/93‑1 (MQ=255)

      cAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTAcgccg                                                                                   >  2:53187/1‑93 (MQ=255)

             tGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCg                                                                            >  2:4012370/1‑93 (MQ=255)

             tGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCg                                                                            >  2:4460308/1‑93 (MQ=255)

              gAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGc                                                                           >  2:2869255/1‑93 (MQ=255)

               aaGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGcc                                                                          <  1:1467929/93‑1 (MQ=255)

                aGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCg                                                                         <  1:2283448/93‑1 (MQ=255)

                    gCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGAt                                                                     <  1:4744520/93‑1 (MQ=255)

                     ccGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATc                                                                    >  1:1293758/1‑93 (MQ=255)

                     ccGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATc                                                                    >  1:6018691/1‑93 (MQ=255)

                      cGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCa                                                                   <  2:2026920/93‑1 (MQ=255)

                       gCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCaa                                                                  >  1:193053/1‑93 (MQ=255)

                       gCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCaa                                                                  >  1:795033/1‑93 (MQ=255)

                       gCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCaa                                                                  >  1:4961539/1‑93 (MQ=255)

                        cAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAAc                                                                 >  1:3097973/1‑93 (MQ=255)

                         aaTGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACa                                                                >  2:5708176/1‑93 (MQ=255)

                         aaTGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACa                                                                >  1:6431663/1‑93 (MQ=255)

                         aaTGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACa                                                                >  1:1764499/1‑93 (MQ=255)

                           tGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGc                                                              <  1:7886861/93‑1 (MQ=255)

                           tGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGc                                                              >  1:2004297/1‑93 (MQ=255)

                           tGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGc                                                              >  2:7263296/1‑93 (MQ=255)

                              aGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCata                                                           >  1:3484127/1‑93 (MQ=255)

                               gAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCatat                                                          <  2:4601902/93‑1 (MQ=255)

                                aaTTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATc                                                         >  1:4453977/1‑93 (MQ=255)

                                aaTTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATc                                                         <  2:5169301/93‑1 (MQ=255)

                                aaTTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCACCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATc                                                         <  2:1801305/93‑1 (MQ=255)

                                  tttGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATcgc                                                       >  1:1425075/1‑93 (MQ=255)

                                   ttGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATcgcg                                                      <  1:3276759/93‑1 (MQ=255)

                                     gCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTg                                                    <  2:3269049/93‑1 (MQ=255)

                                      caccacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGc                                                   <  1:2502425/93‑1 (MQ=255)

                                       accacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCa                                                  <  1:2557177/93‑1 (MQ=255)

                                       accacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCa                                                  >  2:7600788/1‑93 (MQ=255)

                                        ccacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCAc                                                 >  1:2829536/1‑93 (MQ=255)

                                           cTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTAc                                              >  2:1905435/1‑93 (MQ=255)

                                           cTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTAc                                              <  1:4234121/93‑1 (MQ=255)

                                            tATCACCCACGGCAACGCCGTAGTTCCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACa                                             <  2:5857526/93‑1 (MQ=255)

                                            tATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACa                                             <  2:7329883/93‑1 (MQ=255)

                                            tATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACa                                             >  2:5345537/1‑93 (MQ=255)

                                             aTCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACaa                                            >  2:3405575/1‑93 (MQ=255)

                                             aTCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACaa                                            >  2:3831527/1‑93 (MQ=255)

                                             aTCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACaa                                            <  1:438249/93‑1 (MQ=255)

                                                aaCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcgg                                         >  2:3741620/1‑93 (MQ=255)

                                                aaCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcgg                                         <  2:6545566/93‑1 (MQ=255)

                                                 aCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggc                                        >  1:4311462/1‑93 (MQ=255)

                                                 aCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggc                                        >  2:6334268/1‑93 (MQ=255)

                                                  ccACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCAGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggcg                                       <  2:6470934/93‑1 (MQ=255)

                                                   cACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAAcggcgg                                      >  1:703720/1‑93 (MQ=255)

                                                     cGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGaa                                    <  1:6525854/93‑1 (MQ=255)

                                                        cAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAg                                 >  2:4997148/1‑93 (MQ=255)

                                                         aaCGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACAGCGGCAAAGt                                >  2:2741533/1‑93 (MQ=255)

                                                          aCGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATGGCGTGCACTACAACGGCGGGAAAGtt                               >  2:2202416/1‑93 (MQ=255)

                                                          aCGCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAggt                               >  2:261851/1‑91 (MQ=255)

                                                            gCCGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTtta                             <  2:7069565/93‑1 (MQ=255)

                                                             ccGTAGTTGCTCCCATAATGCCCGAGCGCGCTACGCCGCAAATCGCCTAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttat                            <  1:1631738/93‑1 (MQ=255)

                                                             ccGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttat                            <  1:5538936/93‑1 (MQ=255)

                                                             ccGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttat                            >  2:2935663/1‑93 (MQ=255)

                                                             ccGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttat                            <  1:3868431/93‑1 (MQ=255)

                                                              cGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatt                           >  2:5083049/1‑93 (MQ=255)

                                                              cGTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatt                           <  1:1011306/93‑1 (MQ=255)

                                                               gTAGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTttatta                          <  2:6039557/93‑1 (MQ=255)

                                                                 aGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACt                        <  2:6222930/93‑1 (MQ=255)

                                                                 aGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACt                        <  1:1966498/93‑1 (MQ=255)

                                                                 aGTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACt                        >  1:523792/1‑93 (MQ=255)

                                                                  gTTGCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTc                       >  1:727611/1‑93 (MQ=255)

                                                                    tGCTCCCATAATGCCCCAGCGCGATACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAg                     <  2:4501845/93‑1 (MQ=255)

                                                                     gCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAgc                    <  2:7875119/93‑1 (MQ=255)

                                                                     gCTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAgc                    <  1:6190057/93‑1 (MQ=255)

                                                                      cTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAgcg                   <  2:2832459/93‑1 (MQ=255)

                                                                      cTCCCATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAgcg                   >  2:2507409/1‑93 (MQ=255)

                                                                        cccATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCa                 <  1:3693608/93‑1 (MQ=255)

                                                                        cccATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCa                 >  1:7262409/1‑93 (MQ=255)

                                                                        cccATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCa                 >  2:1058145/1‑93 (MQ=255)

                                                                         ccATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAgcgcca                >  1:6104253/1‑91 (MQ=255)

                                                                          cATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCaaa               <  2:2456704/93‑1 (MQ=255)

                                                                          cATAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCaaa               <  2:1552150/93‑1 (MQ=255)

                                                                           aTAATGCCCCAGCGCGCTACGCCGCAAATCGGCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAg              <  1:3241472/93‑1 (MQ=255)

                                                                           aTAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAg              <  1:2804337/93‑1 (MQ=255)

                                                                           aTAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAg              <  2:7449269/93‑1 (MQ=255)

                                                                           aTAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAg              <  2:1201680/93‑1 (MQ=255)

                                                                            tAATGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGa             <  2:3101239/93‑1 (MQ=255)

                                                                              aTGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATa           <  1:6910374/93‑1 (MQ=255)

                                                                              aTGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATa           <  1:2758076/93‑1 (MQ=255)

                                                                               tGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa          <  1:5115532/93‑1 (MQ=255)

                                                                               tGCCCCAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa          <  2:6607000/93‑1 (MQ=255)

                                                                               agCCCCAGCGCGCTCCGCCGCAAAGCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATaa          <  2:7692748/92‑1 (MQ=255)

                                                                                 ccccAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATg        <  1:6611682/93‑1 (MQ=255)

                                                                                   ccAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGTACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGAc      <  1:5906354/93‑1 (MQ=255)

                                                                                    cAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACt     <  2:4689602/93‑1 (MQ=255)

                                                                                    cAGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACt     <  2:6790122/93‑1 (MQ=255)

                                                                                     aGCGCGCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACtt    >  2:4131879/1‑93 (MQ=255)

                                                                                      gcgcgcTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTc   >  1:206540/1‑93 (MQ=255)

                                                                                      gcgcgcTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTc   >  2:6398765/1‑93 (MQ=255)

                                                                                      gcgcgcTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTc   <  2:3191075/93‑1 (MQ=255)

                                                                                       cgcgCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTcc  <  2:1611273/93‑1 (MQ=255)

                                                                                       cgcgCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTcc  >  2:3381816/1‑93 (MQ=255)

AGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCC  >  NC_000913/2664941‑2665120

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