BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,071,901	G→A	A197T (GCC→ACC)	yihV →	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,071,901	0	G	A	99.4%	1267.9 / ‑4.2	356	A197T (GCC→ACC)	yihV	putative kinase
Reads supporting (aligned to +/- strand): ref base G (0/2); new base A (184/170); total (184/172)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 361 total aligned reads displayed.

CCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGGCCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAATGGT  >  NC_000913/4071811‑4071990

ccgcAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGAcc                                                                                         >  2:3354265/1‑93 (MQ=255)

ccgcAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGAcc                                                                                         <  2:5297533/93‑1 (MQ=255)

 cgcAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCa                                                                                        <  1:1147210/93‑1 (MQ=255)

 cgcAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCa                                                                                        <  2:5397470/93‑1 (MQ=255)

 cgcAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCa                                                                                        <  2:8651179/93‑1 (MQ=255)

  gcAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAg                                                                                       <  2:1552752/93‑1 (MQ=255)

  gcAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAg                                                                                       <  2:2702557/93‑1 (MQ=255)

   cAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGt                                                                                      >  2:3242165/1‑93 (MQ=255)

    aGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTg                                                                                     <  1:5723780/93‑1 (MQ=255)

    aGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTg                                                                                     >  1:5369890/1‑93 (MQ=255)

    aGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTGTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTg                                                                                     <  1:7325414/93‑1 (MQ=255)

     ggATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTgc                                                                                    <  1:2881137/93‑1 (MQ=255)

      gATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTgcg                                                                                   >  1:3054189/1‑93 (MQ=255)

      gATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTgcg                                                                                   <  2:9162994/93‑1 (MQ=255)

       atatCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTgcgc                                                                                  >  2:8171758/1‑93 (MQ=255)

        tatCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCt                                                                                 >  1:4927318/1‑93 (MQ=255)

        tatCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCt                                                                                 >  2:9526367/1‑93 (MQ=255)

         atCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTa                                                                                <  2:1230257/93‑1 (MQ=255)

         atCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTa                                                                                >  2:4976714/1‑93 (MQ=255)

            aGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTaaaa                                                                             <  1:7957883/93‑1 (MQ=255)

             gTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAAc                                                                            >  2:7871469/1‑93 (MQ=255)

             gTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAAc                                                                            <  2:3701466/93‑1 (MQ=255)

             gTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTCAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAAc                                                                            <  1:1802053/93‑1 (MQ=255)

               gAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAg                                                                          >  1:7243478/1‑93 (MQ=255)

                aGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAgg                                                                         <  2:9377701/93‑1 (MQ=255)

                 gCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGc                                                                        >  2:5835520/1‑93 (MQ=255)

                   tggtggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGcac                                                                      <  2:4865673/93‑1 (MQ=255)

                   tggtggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGcac                                                                      >  2:8943524/1‑93 (MQ=255)

                    ggtggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGcaca                                                                     <  2:7443984/93‑1 (MQ=255)

                    ggtggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGGGCGCTAAAACAGGcaca                                                                     >  2:2496641/1‑93 (MQ=255)

                      tggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGCAAGAGATGAACAGTGCGCTAACACAGGCACaaa                                                                   >  2:7062843/1‑93 (MQ=255)

                      tggCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACaaa                                                                   <  2:9319788/93‑1 (MQ=255)

                       ggCATTAAGCGATCACGCGGCCTTTTCAGAACCAGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAAc                                                                  <  1:7138649/93‑1 (MQ=255)

                            tAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTca                                                             >  1:6064347/1‑93 (MQ=255)

                             aaGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTcac                                                            <  1:7744169/93‑1 (MQ=255)

                             aaGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTcac                                                            >  2:2238266/1‑93 (MQ=255)

                             aaGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTcac                                                            <  2:6321004/93‑1 (MQ=255)

                                cGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACaaa                                                         <  2:4293504/93‑1 (MQ=255)

                                cGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACaaa                                                         <  1:1587002/93‑1 (MQ=255)

                                cGATCACGCGGCCTTGTCAGAACCGGGTCTGGCGCGCTTACCGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACaaa                                                         <  1:6440312/93‑1 (MQ=255)

                                 gATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAAt                                                        >  1:6656338/1‑93 (MQ=255)

                                  aTCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAGGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATg                                                       >  2:3540557/1‑93 (MQ=255)

                                  aTCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATg                                                       <  2:7657517/93‑1 (MQ=255)

                                  aTCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATg                                                       >  1:9109474/1‑93 (MQ=255)

                                       gcgGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCAt                                                  >  1:6771923/1‑93 (MQ=255)

                                        cgGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATg                                                 >  1:9672007/1‑93 (MQ=255)

                                         ggCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGt                                                >  1:8656947/1‑93 (MQ=255)

                                          gCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTc                                               >  2:6057307/1‑93 (MQ=255)

                                          gCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTc                                               >  2:189916/1‑93 (MQ=255)

                                           ccTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCt                                              >  1:3718270/1‑93 (MQ=255)

                                            cTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTa                                             <  1:463408/93‑1 (MQ=255)

                                             ttttCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCGGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAt                                            >  1:5273152/1‑93 (MQ=255)

                                             ttttCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAt                                            <  2:9771762/93‑1 (MQ=255)

                                             ttttCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAt                                            >  2:8021576/1‑93 (MQ=255)

                                             ttttCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAt                                            >  2:8425541/1‑93 (MQ=255)

                                              tttCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAtg                                           >  2:1305190/1‑93 (MQ=255)

                                              tttCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTAtg                                           <  1:7561427/93‑1 (MQ=255)

                                                 cAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGa                                        <  2:8063724/93‑1 (MQ=255)

                                                    aaCCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGAccc                                     <  2:4140524/93‑1 (MQ=255)

                                                     aCCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGAcccg                                    >  1:8427524/1‑92 (MQ=255)

                                                      ccGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAg                                   <  1:8024941/93‑1 (MQ=255)

                                                       cGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAgg                                  <  2:4029344/93‑1 (MQ=255)

                                                       cGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAgg                                  <  2:83169/93‑1 (MQ=255)

                                                        gggTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAggg                                 <  1:2037235/93‑1 (MQ=255)

                                                        gggTCTGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAggg                                 <  2:6907532/93‑1 (MQ=255)

                                                          gTCTGGCGCGCTTAACGGGCGTAAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAgggaa                               <  1:8946140/93‑3 (MQ=255)

                                                             tGGCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGGAgcg                            >  1:2243005/1‑93 (MQ=255)

                                                              ggCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAgcgc                           <  1:3314636/93‑1 (MQ=255)

                                                              ggCGCGCTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAgcgc                           <  2:6967585/93‑1 (MQ=255)

                                                               gcgcgcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCa                          <  2:6472174/93‑1 (MQ=255)

                                                                 gcgcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAgg                        >  1:9651003/1‑93 (MQ=255)

                                                                  cgcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGc                       <  2:8649394/93‑1 (MQ=255)

                                                                  cgcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGc                       <  2:9639082/93‑1 (MQ=255)

                                                                   gcTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCt                      >  1:3853419/1‑93 (MQ=255)

                                                                    cTTAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTg                     <  2:980908/93‑1 (MQ=255)

                                                                     ttAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGc                    <  2:7579273/93‑1 (MQ=255)

                                                                      tAACGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCg                   <  1:766590/93‑1 (MQ=255)

                                                                       aaCGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGa                  <  1:8728404/93‑1 (MQ=255)

                                                                        aCGGGCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGAc                 <  2:1400613/93‑1 (MQ=255)

                                                                          gggCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActg               <  1:354359/93‑1 (MQ=255)

                                                                            gCGTGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggc             >  2:8715397/1‑93 (MQ=255)

                                                                               tGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggctgg          >  1:6222695/1‑93 (MQ=255)

                                                                               tGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggctgg          >  1:10022314/1‑93 (MQ=255)

                                                                               tGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggctgg          >  2:2352994/1‑93 (MQ=255)

                                                                               tGAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGActggctgg          <  1:6805047/93‑1 (MQ=255)

                                                                                gAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGa         >  1:4462845/1‑93 (MQ=255)

                                                                                gAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGa         >  2:6781880/1‑93 (MQ=255)

                                                                                gAAAGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGa         >  2:3906546/1‑93 (MQ=255)

                                                                                 aaaGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGaa        >  1:4968351/1‑93 (MQ=255)

                                                                                 aaaGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGaa        >  2:4564679/1‑93 (MQ=255)

                                                                                 aaaGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGaa        >  2:2413576/1‑93 (MQ=255)

                                                                                  aaGAGATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGaaa       >  1:4145011/1‑93 (MQ=255)

                                                                                   agagATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGaaaa      <  1:5903033/93‑1 (MQ=255)

                                                                                   agagATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGaaaa      <  1:9369053/93‑1 (MQ=255)

                                                                                   agagATGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGaaaa      >  1:1434515/1‑93 (MQ=255)

                                                                                    gagaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAt     <  2:4325454/93‑1 (MQ=255)

                                                                                    gagaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAt     >  2:7313583/1‑93 (MQ=255)

                                                                                      gaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGACAAtgg   >  1:6270540/1‑93 (MQ=255)

                                                                                      gaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAtgg   <  1:2451679/93‑1 (MQ=255)

                                                                                      gaTGACCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAAtgg   >  2:2127737/1‑93 (MQ=255)

CCGCAGGATATCAGTGAGCTGGTGGCATTAAGCGATCACGCGGCCTTTTCAGAACCGGGTCTGGCGCGCTTAACGGGCGTGAAAGAGATGGCCAGTGCGCTAAAACAGGCACAAACGCTCACAAATGGTCATGTCTATGTGACCCAGGGTAGCGCAGGCTGCGACTGGCTGGAAAATGGT  >  NC_000913/4071811‑4071990

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