BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,418,843	A→G	K46E (AAG→GAG)	ydaT →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,418,843	0	A	G	96.4%	396.0 / 3.5	137	K46E (AAG→GAG)	ydaT	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (4/1); new base G (47/85); total (51/86)
Fisher's exact test for biased strand distribution p-value = 6.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

Only 100 of 142 total aligned reads displayed.

GCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCA  >  NC_000913/1418753‑1418932

gcagcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAg                                                                                         <  2:2506631/93‑1 (MQ=255)

gcagcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAg                                                                                         <  2:113291/93‑1 (MQ=255)

  agcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGACGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGa                                                                                       >  1:2519202/1‑93 (MQ=255)

   gcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGAc                                                                                      <  2:1027975/93‑1 (MQ=255)

   gcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGAc                                                                                      <  2:3152379/93‑1 (MQ=255)

   gcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGAc                                                                                      >  1:3565074/1‑93 (MQ=255)

   gcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGAc                                                                                      <  2:3784378/93‑1 (MQ=255)

   gcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGAc                                                                                      <  1:3472862/93‑1 (MQ=255)

   gcCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCACGCGAGGAc                                                                                      <  2:4581269/93‑1 (MQ=255)

    ccGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACt                                                                                     <  1:3244506/93‑1 (MQ=255)

     cGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTg                                                                                    <  2:920503/93‑1 (MQ=255)

      gAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTgg                                                                                   >  1:2374881/1‑93 (MQ=255)

      gAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTgg                                                                                   >  1:2684887/1‑93 (MQ=255)

      gAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTgg                                                                                   >  2:2801013/1‑93 (MQ=255)

           aGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTGCAGGCGAGGACTGGAACaa                                                                              >  2:1438387/1‑93 (MQ=255)

           aGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACaa                                                                              >  1:2366717/1‑93 (MQ=255)

             ggcaggcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATc                                                                            <  1:3737566/93‑1 (MQ=255)

             ggcaggcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATc                                                                            <  2:3770791/93‑1 (MQ=255)

              gcaggcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATca                                                                           <  1:3762925/93‑1 (MQ=255)

              gcaggcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATca                                                                           <  1:3052098/93‑1 (MQ=255)

                aggcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATcagc                                                                         <  1:3443666/93‑1 (MQ=255)

                aggcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATcagc                                                                         <  2:3398102/93‑1 (MQ=255)

                 tgcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCACGCGAGGACTGGAACAATcagca                                                                        >  2:1173142/2‑93 (MQ=255)

                 ggcaTGGGTAGCCAATGGAATTACTGAAGAATATCTGCGCCAGGGGGGCGACGAATTGGCCCTGGTTCCAGGCGAGGACTGGAACAATcagca                                                                        <  2:2991334/93‑1 (MQ=255)

                 ggcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATcagca                                                                        >  2:4297623/1‑93 (MQ=255)

                 ggcaTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATcagca                                                                        >  1:1739380/1‑93 (MQ=255)

                   caTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGa                                                                      <  1:3021236/93‑1 (MQ=255)

                    aTGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGaa                                                                     <  1:2422005/93‑1 (MQ=255)

                     tGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAat                                                                    <  1:1578919/93‑1 (MQ=255)

                     tGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTAGAACAATCAGCAGAat                                                                    >  2:3949648/1‑93 (MQ=255)

                      gggTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAata                                                                   <  2:828929/93‑1 (MQ=255)

                       ggTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAatat                                                                  <  1:3861272/93‑1 (MQ=255)

                       ggTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAatat                                                                  >  2:133088/1‑93 (MQ=255)

                        gTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGCATTGCCTCTGGCTTCAGGCGAGGACCGGAACAAAAAGGCGAatata                                                                 >  2:1448942/1‑92 (MQ=255)

                        gTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAAtatc                                                                 <  2:3961693/93‑1 (MQ=255)

                           gCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAAtatctat                                                              <  1:729566/93‑1 (MQ=255)

                            ccAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAAtatctatc                                                             <  2:2242726/93‑1 (MQ=255)

                            ccAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAAtatctatc                                                             <  2:3598457/93‑1 (MQ=255)

                            ccAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAAtatctatc                                                             <  1:2670074/93‑1 (MQ=255)

                             cAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCa                                                            >  1:3356294/1‑93 (MQ=255)

                              aaTGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGGATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCAc                                                           >  2:3826565/1‑93 (MQ=255)

                              aaTGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCAc                                                           >  2:3264463/1‑93 (MQ=255)

                                  cAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCgtt                                                       <  2:1147491/93‑1 (MQ=255)

                                   aaTTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCgttg                                                      >  2:4336472/1‑93 (MQ=255)

                                     ttACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCgttggt                                                    >  1:1831273/1‑93 (MQ=255)

                                        cTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGa                                                 >  2:1630826/1‑93 (MQ=255)

                                        cTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGa                                                 >  2:416639/1‑93 (MQ=255)

                                         tGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGaa                                                >  2:472451/1‑93 (MQ=255)

                                         tGAAGAACATCTGCCCCAGGGTGGCGGCAAATTGCCCCTGGTGCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGaa                                                <  1:2997355/93‑1 (MQ=255)

                                          gaagaaTATCTGTGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCATCGTTGGTTGaaa                                               >  2:2961671/1‑93 (MQ=255)

                                          gaagaaTATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGaaa                                               <  1:3764216/93‑1 (MQ=255)

                                           aagaaTATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAg                                              >  2:2252542/1‑93 (MQ=255)

                                            agaaTATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAgg                                             >  2:3282475/1‑93 (MQ=255)

                                              aaTATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTg                                           >  1:4382888/1‑93 (MQ=255)

                                               atatCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGa                                          >  1:2790358/1‑93 (MQ=255)

                                                 atCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGaaa                                        >  1:4678811/1‑93 (MQ=255)

                                                 atCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGaaa                                        <  1:3033585/93‑1 (MQ=255)

                                                 agcTGCGCCAGGGGGGCGGCGAATTGCCCGTGGTTCCAGGCGAGGACTGTAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGaaa                                        <  1:1507270/91‑1 (MQ=255)

                                                  tCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAAc                                       >  2:1663282/1‑93 (MQ=255)

                                                  tCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAAc                                       <  1:2236419/93‑1 (MQ=255)

                                                  tCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAAc                                       <  1:3196293/93‑1 (MQ=255)

                                                      cgcCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGaaa                                   <  1:1492244/93‑1 (MQ=255)

                                                       gcCAGGGGGGCGTCGAATTGCCCCTGTTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGaaaa                                  <  1:416113/93‑1 (MQ=255)

                                                         cAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACg                                >  2:3168125/1‑93 (MQ=255)

                                                         cAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACg                                <  1:3689506/93‑1 (MQ=255)

                                                           ggggggCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCa                              >  1:57209/1‑93 (MQ=255)

                                                           ggggggCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCa                              >  1:1681090/1‑93 (MQ=255)

                                                           ggggggCGGCGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCa                              <  2:1098311/93‑1 (MQ=255)

                                                               ggcggcGAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAaga                          <  1:3720523/93‑1 (MQ=255)

                                                               ggcggcGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAaga                          <  1:4078211/93‑1 (MQ=255)

                                                                gcggcgAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAagag                         <  2:72599/93‑1 (MQ=255)

                                                                gcggcgAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAagag                         <  1:3603759/93‑1 (MQ=255)

                                                                gcgacgAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTTAAAGGTGAAACGAAAACGCAAagag                         <  2:3488800/93‑1 (MQ=255)

                                                                 cggcgAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAgaga                        <  1:2780480/93‑1 (MQ=255)

                                                                 cggcgAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAgaga                        <  1:4135628/93‑1 (MQ=255)

                                                                  ggcgAATTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGaa                       >  1:350515/1‑93 (MQ=255)

                                                                    cgAATTGCCCCGGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGaaaa                     <  2:4529748/93‑1 (MQ=255)

                                                                      aaTTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGaaaaaa                   <  2:2405798/93‑1 (MQ=255)

                                                                      aaTTGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGaaaaaa                   <  1:3869666/93‑1 (MQ=255)

                                                                        ttGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAAtt                 <  2:298777/93‑1 (MQ=255)

                                                                         tGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTc                <  1:1403581/93‑1 (MQ=255)

                                                                         tGCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTc                <  2:1614560/93‑1 (MQ=255)

                                                                          gCCCCTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCa               <  2:4002638/93‑1 (MQ=255)

                                                                            cccTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGa             <  2:1737615/93‑1 (MQ=255)

                                                                            cccTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGa             <  1:1157056/93‑1 (MQ=255)

                                                                            cccTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGa             <  1:1435427/93‑1 (MQ=255)

                                                                            cccTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGa             >  2:16761/1‑93 (MQ=255)

                                                                            cccTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGa             <  1:2015269/93‑1 (MQ=255)

                                                                             ccTGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGaa            >  1:3339320/1‑93 (MQ=255)

                                                                               tGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGc          <  1:2416824/93‑1 (MQ=255)

                                                                               tGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGc          >  2:229558/1‑93 (MQ=255)

                                                                               tGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGc          >  1:714981/1‑93 (MQ=255)

                                                                               tGGTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGc          >  1:4646212/1‑93 (MQ=255)

                                                                                ggTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCt         <  1:4374439/93‑1 (MQ=255)

                                                                                 gTTCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTg        <  1:1797254/93‑1 (MQ=255)

                                                                                  ttCCAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGa       <  1:195507/93‑1 (MQ=255)

                                                                                    ccAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATc     <  2:2873331/93‑1 (MQ=255)

                                                                                     cAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATcc    >  2:1525702/1‑93 (MQ=255)

                                                                                     cAGGCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATcc    <  2:723355/93‑1 (MQ=255)

                                                                                       ggCGAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCa  <  1:3014324/93‑1 (MQ=255)

GCAGCCGAAAAAGGGCAGGCATGGGTAGCCAATGCAATTACTGAAGAATATCTGCGCCAGGGGGGCGGCGAATTGCCCCTGGTTCCAGGCAAGGACTGGAACAATCAGCAGAATATCTATCACCGTTGGTTGAAAGGTGAAACGAAAACGCAAAGAGAAAAAATTCAGAAGCTGATCCCA  >  NC_000913/1418753‑1418932

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