BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,945,896	T→C	V63A (GTA→GCA)	yifE →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,945,896	0	T	C	99.3%	875.9 / ‑3.2	274	V63A (GTA→GCA)	yifE	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/2); new base C (136/135); total (137/137)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-01

Only 100 of 286 total aligned reads displayed.

AACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACC  >  NC_000913/3945806‑3945987

aaCTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCAt                                                                                           >  2:619273/1‑93 (MQ=255)

aaCTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCAt                                                                                           <  1:3264615/93‑1 (MQ=255)

  cTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGc                                                                                         >  1:2397824/1‑93 (MQ=255)

     cTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGt                                                                                      <  1:3583771/93‑1 (MQ=255)

      ttGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTg                                                                                     <  2:4013963/93‑1 (MQ=255)

       tGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTgg                                                                                    >  1:1148980/1‑93 (MQ=255)

        gAGGGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGc                                                                                   >  1:7423/1‑93 (MQ=255)

        gAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGc                                                                                   >  1:2909056/1‑93 (MQ=255)

          gCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGa                                                                                 >  1:2420542/1‑93 (MQ=255)

          gCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGa                                                                                 <  2:4155351/93‑1 (MQ=255)

            gTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAAc                                                                               >  1:3626139/1‑93 (MQ=255)

             tCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACg                                                                              >  1:3563709/1‑93 (MQ=255)

             tCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACg                                                                              >  2:1704343/1‑93 (MQ=255)

               aTGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTg                                                                            >  2:2980551/1‑93 (MQ=255)

                tGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGa                                                                           <  2:975802/93‑1 (MQ=255)

                 ggTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAg                                                                          <  2:1825284/93‑1 (MQ=255)

                   ttATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGcc                                                                        >  1:4642038/1‑93 (MQ=255)

                     aTGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAg                                                                      >  1:4244139/1‑93 (MQ=255)

                     aTGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAg                                                                      <  2:4411626/93‑1 (MQ=255)

                       gCCTTCAATGAGTTGGATCTTGGCAAGCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTg                                                                    <  1:4405869/93‑1 (MQ=255)

                       gCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTg                                                                    >  2:517243/1‑93 (MQ=255)

                        ccTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTgc                                                                   >  2:1119676/1‑92 (MQ=255)

                         cTTCAATGAGTTGGATCTGGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGAc                                                                  <  1:2838828/93‑1 (MQ=255)

                          ttCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACa                                                                 >  2:4429323/1‑93 (MQ=255)

                            cAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGa                                                               <  2:1555792/93‑1 (MQ=255)

                            cAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGa                                                               <  1:2133815/93‑1 (MQ=255)

                            cAATGAGTGGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGa                                                               <  2:4613780/93‑1 (MQ=255)

                             aaTCAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGaa                                                              <  1:165402/93‑1 (MQ=255)

                                 aGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAg                                                          >  2:1924241/1‑93 (MQ=255)

                                  gTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGGCGTGGCGAACGTGAGCCAGTGACAGAAGCAGa                                                         >  2:2082287/1‑93 (MQ=255)

                                  gTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGa                                                         <  2:590120/93‑1 (MQ=255)

                                   ttGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGaa                                                        <  1:679785/93‑1 (MQ=255)

                                    tGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAc                                                       >  1:4520484/1‑93 (MQ=255)

                                    tGGATCATGGAAACCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAc                                                       <  1:1905460/93‑1 (MQ=255)

                                     ggATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcg                                                      >  1:1074290/1‑93 (MQ=255)

                                     ggATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcg                                                      >  2:2040497/1‑93 (MQ=255)

                                     ggATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcg                                                      >  1:3416831/1‑93 (MQ=255)

                                       aTCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcgcg                                                    >  2:3215892/1‑93 (MQ=255)

                                       aTCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAAcgcg                                                    <  1:3832445/93‑1 (MQ=255)

                                        tCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgt                                                   >  2:3343171/1‑93 (MQ=255)

                                        tCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgt                                                   <  1:2647665/93‑1 (MQ=255)

                                          ttGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgtgt                                                 <  2:3679814/93‑1 (MQ=255)

                                          ttGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgtgt                                                 >  1:3526612/1‑93 (MQ=255)

                                          ttGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCgtgt                                                 >  1:825929/1‑93 (MQ=255)

                                           tGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGtgtg                                                >  2:1322707/1‑93 (MQ=255)

                                           tGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGtgtg                                                <  2:2242948/93‑1 (MQ=255)

                                             gCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGt                                              <  2:47164/93‑1 (MQ=255)

                                              cAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTc                                             >  1:4302381/1‑93 (MQ=255)

                                               aaaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTcc                                            <  2:3006858/93‑1 (MQ=255)

                                               aaaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTcc                                            >  2:265884/1‑93 (MQ=255)

                                                aaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCa                                           >  2:480306/1‑93 (MQ=255)

                                                aaCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGCTCCa                                           >  1:3174671/1‑93 (MQ=255)

                                                 aCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCaa                                          >  1:502933/1‑93 (MQ=255)

                                                 aCGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCaa                                          <  2:3318312/93‑1 (MQ=255)

                                                   gcgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCTAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGt                                        <  2:2773624/93‑1 (MQ=255)

                                                   gcgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGt                                        >  1:388933/1‑93 (MQ=255)

                                                    cgAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGta                                       <  2:2393090/93‑1 (MQ=255)

                                                     gAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGtat                                      >  1:632534/1‑93 (MQ=255)

                                                       gCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTatat                                    >  1:261737/1‑93 (MQ=255)

                                                        ccGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATg                                   <  2:3909894/93‑1 (MQ=255)

                                                           gTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcg                                <  1:3345840/93‑1 (MQ=255)

                                                           gTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcg                                >  2:832399/1‑93 (MQ=255)

                                                           gTTACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcg                                >  2:4600004/1‑93 (MQ=255)

                                                            ttgccGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcgc                               <  2:350172/90‑1 (MQ=255)

                                                             tACCGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGAcgcg                              >  1:4090289/1‑93 (MQ=255)

                                                               ccGAGGAAGAGAACCGCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTa                            <  1:2009190/93‑1 (MQ=255)

                                                               ccGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTa                            >  2:1783021/1‑93 (MQ=255)

                                                               ccGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTa                            >  1:3704141/1‑93 (MQ=255)

                                                                cGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTAt                           <  2:3974458/93‑1 (MQ=255)

                                                                cGAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTAt                           <  2:3799346/93‑1 (MQ=255)

                                                                 gAGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTAtt                          <  1:3120456/93‑1 (MQ=255)

                                                                  aGGAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTa                         <  1:988770/93‑1 (MQ=255)

                                                                   ggAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGCCCAAGTATATGACGCGTATTag                        >  2:3535180/1‑92 (MQ=255)

                                                                   ggAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGACCAAGTATATGACGAGTATTaa                        >  1:1821558/1‑93 (MQ=255)

                                                                    gAAGAGAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAg                       >  1:1010142/1‑93 (MQ=255)

                                                                      agagAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCg                     >  1:1437352/1‑93 (MQ=255)

                                                                      agagAAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCg                     <  2:1300508/93‑1 (MQ=255)

                                                                        agaAACTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTc                   >  2:2738352/1‑93 (MQ=255)

                                                                          aaaCTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCa                 <  1:2580466/93‑1 (MQ=255)

                                                                          aaaCTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCa                 >  1:1531789/1‑93 (MQ=255)

                                                                           aaCTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaa                >  1:4542405/1‑93 (MQ=255)

                                                                            aCTCTTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaaa               <  2:673008/93‑1 (MQ=255)

                                                                             ctctTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaaaa              <  1:3218981/93‑1 (MQ=255)

                                                                             ctctTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCaaaa              <  2:4359168/93‑1 (MQ=255)

                                                                              tctTCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAAc             >  2:3433588/1‑93 (MQ=255)

                                                                                ttCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGt           <  1:4590798/93‑1 (MQ=255)

                                                                                ttCGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGt           <  1:2095461/93‑1 (MQ=255)

                                                                                  cGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGttt         >  1:4002296/1‑93 (MQ=255)

                                                                                  cGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGttt         <  1:1604483/93‑1 (MQ=255)

                                                                                  cGTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGttt         >  1:87141/1‑93 (MQ=255)

                                                                                   gTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGtttt        <  2:4062331/93‑1 (MQ=255)

                                                                                   gTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGtttt        >  2:1517647/1‑93 (MQ=255)

                                                                                   gTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGtttt        <  1:2234933/93‑1 (MQ=255)

                                                                                   gTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGtttt        >  1:2174757/1‑93 (MQ=255)

                                                                                   gTAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGtttt        >  2:3167285/1‑93 (MQ=255)

                                                                                    tAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTc       <  2:1424623/93‑1 (MQ=255)

                                                                                    tAGCAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTc       <  1:1333186/93‑1 (MQ=255)

                                                                                     agcagcATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTca      <  1:1693743/93‑1 (MQ=255)

                                                                                      gcagcaTGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGAGTCCAAAACGTTTTcac     >  2:1477633/1‑93 (MQ=255)

                                                                                       cAGCATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCCAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTcaca    >  1:2271989/1‑93 (MQ=255)

AACTGCTTGAGCGTCATGGTTATGCCTTCAATGAGTTGGATCTTGGCAAACGCGAGCCGGTTACCGAGGAAGAGAAACTCTTCGTAGCAGTATGCCGTGGCGAACGTGAGCCAGTGACAGAAGCAGAACGCGTGTGGTCCAAGTATATGACGCGTATTAAGCGTCCAAAACGTTTTCACACC  >  NC_000913/3945806‑3945987

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