BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,856,333	T→G	15.6%	E199A (GAG→GCG) ‡	ydjH ←	putative kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,856,333	0	T	G	15.6%	341.0 / 51.3	136	E199A (GAG→GCG)	ydjH	putative kinase
Reads supporting (aligned to +/- strand): ref base T (68/46); new base G (9/13); total (77/59)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

CGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAA  >  NC_000913/1856235‑1856432

cGTCTTTACCCGTTTTAATCACCACCGTTTTTTTTTTTTTCTCTACAAAGCAGCCAGCAATTTCATCCATTTTTTCTTTCCCGCTTAGTAATTTTGCCt                                                                                                     >  2:3356804/1‑99 (MQ=255)

cGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTc                                                                                                    >  2:3198630/1‑100 (MQ=255)

 gTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCg                                                                                                   >  1:2046455/1‑100 (MQ=255)

  tCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCgg                                                                                                  <  1:1414857/100‑1 (MQ=255)

  tCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCgg                                                                                                  <  2:2792672/100‑1 (MQ=255)

   cTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGc                                                                                                 <  1:3545097/100‑1 (MQ=255)

   cTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGc                                                                                                 <  1:1211253/100‑1 (MQ=255)

   cTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGc                                                                                                 >  1:3498309/1‑100 (MQ=255)

      tACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCaaa                                                                                              >  2:3831545/1‑100 (MQ=255)

        cccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAAt                                                                                            <  2:2243769/100‑1 (MQ=255)

        cccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAAt                                                                                            >  2:1012256/1‑100 (MQ=255)

        cccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAAt                                                                                            <  1:2216060/100‑1 (MQ=255)

        cccGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCATTCCCGGTGGGTAATTTTGCCTCGGCaaaa                                                                                             >  2:2827017/1‑99 (MQ=255)

        cccGTTTTAATCACCACCGTTTTTACACCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAAt                                                                                            >  1:2552612/1‑100 (MQ=255)

          cGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTa                                                                                          <  1:1544836/100‑1 (MQ=255)

          cGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAAtt                                                                                           <  2:1945218/99‑1 (MQ=255)

           gTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAg                                                                                         >  1:1655321/1‑100 (MQ=255)

           gTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTac                                                                                         <  1:1772333/100‑2 (MQ=255)

           gTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTac                                                                                         <  2:1956931/100‑2 (MQ=255)

           gTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATGTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTac                                                                                         <  2:1426584/100‑2 (MQ=255)

             tttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGa                                                                                       <  1:1849911/100‑1 (MQ=255)

             tttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGa                                                                                       <  1:1260818/100‑1 (MQ=255)

              ttAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGaa                                                                                      <  1:2168501/100‑1 (MQ=255)

               tAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCGAGGCAAAATTactcac                                                                                     <  1:2856677/100‑6 (MQ=39)

                aaTCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAAc                                                                                    >  1:3285102/1‑100 (MQ=255)

                aaTCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAAc                                                                                    <  1:3014698/100‑1 (MQ=255)

                aaTCACCACCATTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAAc                                                                                    <  1:4110810/100‑1 (MQ=255)

                 aTCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACa                                                                                   <  1:3831570/100‑1 (MQ=255)

                 aTCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACa                                                                                   <  1:1031867/100‑1 (MQ=255)

                  ncaccacCGTTTTTACGCCGGACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAg                                                                                  >  2:2557519/2‑100 (MQ=255)

                   caccacCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGa                                                                                 >  2:43614/1‑100 (MQ=255)

                   caccacCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGa                                                                                 <  2:2885600/100‑1 (MQ=255)

                   caccacCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGa                                                                                 >  2:2671651/1‑100 (MQ=255)

                   caccacCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGa                                                                                 >  2:284585/1‑100 (MQ=255)

                   caccacCGATTTTACGCCGCACGCAAGAAAGCAGTCAACAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAACATTAGGAAACAGa                                                                                 >  2:2248821/1‑100 (MQ=255)

                    accaccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGAt                                                                                >  1:4186038/1‑100 (MQ=255)

                    accaccGATTTCACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGAt                                                                                >  2:2206178/1‑100 (MQ=255)

                     ccaccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATa                                                                               >  1:671927/1‑100 (MQ=255)

                     ccaccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATa                                                                               >  2:1123681/1‑100 (MQ=255)

                     ccaccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATa                                                                               >  1:987206/1‑100 (MQ=255)

                      caccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATaa                                                                              >  1:2321126/1‑100 (MQ=255)

                      caccGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATaa                                                                              >  1:3654441/1‑100 (MQ=255)

                         cGTTTTTACGCTGCATGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATc                                                                            >  2:2288603/1‑99 (MQ=255)

                         cGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCg                                                                           <  1:2385622/100‑1 (MQ=255)

                         cGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAACTTCGCCTCGGCAAAATTAGGAAACAGATAATCg                                                                           <  2:2751052/100‑1 (MQ=255)

                          gTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGa                                                                          >  1:3278705/1‑100 (MQ=255)

                          gTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGa                                                                          <  1:4132027/100‑1 (MQ=255)

                                gcgccgcACGCAAGAAAGCAGTCAGCAATGTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAAc                                                                    <  1:3068161/99‑1 (MQ=255)

                                aCGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAAc                                                                    <  1:2206183/100‑1 (MQ=255)

                                 cgccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACt                                                                   >  1:3528298/1‑100 (MQ=255)

                                  gccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTc                                                                  >  1:114257/1‑100 (MQ=255)

                                  gccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTc                                                                  >  2:58771/1‑100 (MQ=255)

                                   ccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCa                                                                 <  1:1596858/100‑1 (MQ=255)

                                   ccgcACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCa                                                                 >  2:293442/1‑100 (MQ=255)

                                     gcacgcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTGGGCAAAATTAGGAAACAGATAATCGACATAACTCAAt                                                               <  2:2249145/100‑1 (MQ=255)

                                     gcacgcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAAt                                                               >  1:1485853/1‑100 (MQ=255)

                                       atgctttaaaGCTTTCAGCTATTTAATCCAGTGTCCCTTTCCCGGTGGGTAATTTTGCCGCGGAAAGTTTAGGTAACAGATAATCGACATAACTCTAtt                                                              <  2:3921240/92‑2 (MQ=255)

                                        cgcaAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCt                                                            >  1:2725425/1‑100 (MQ=255)

                                            aGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGc                                                        <  2:3708119/100‑1 (MQ=255)

                                             gAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCTCa                                                       <  2:2151194/100‑1 (MQ=255)

                                             gAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCa                                                       <  2:1418169/100‑1 (MQ=255)

                                               aaGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCGAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCaaa                                                     >  2:2688407/1‑100 (MQ=255)

                                               aaGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCaaa                                                     >  1:2746839/1‑100 (MQ=255)

                                                    gTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatca                                                >  1:39650/1‑100 (MQ=255)

                                                    gTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCACAAATatca                                                >  2:2903223/1‑100 (MQ=255)

                                                      cAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGCAACAGATAATCGACATAACTCAATGCTTCGCAAATatcatc                                              >  1:647914/1‑100 (MQ=255)

                                                      cAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATatcatc                                              >  1:3100947/1‑100 (MQ=255)

                                                      cAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTACCTCGGCAAAATTAGGAAAGAGATACTCGACATAACTCAATGCTTCGCAAATatcat                                               >  2:2289893/1‑99 (MQ=255)

                                                       aGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATcc                                             >  2:1107610/1‑100 (MQ=255)

                                                        gCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCa                                            >  2:964650/1‑100 (MQ=255)

                                                           aTTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCg                                         <  1:1286475/100‑1 (MQ=255)

                                                           aTTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCAGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCg                                         >  1:3653147/1‑100 (MQ=255)

                                                              tCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGttt                                      >  1:2167519/1‑100 (MQ=255)

                                                               cATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttc                                     <  2:39650/100‑1 (MQ=255)

                                                                 tCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcat                                   <  1:1720417/100‑1 (MQ=255)

                                                                   cAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTttcattc                                 >  2:1814/1‑100 (MQ=255)

                                                                          tctTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACcgcg                          >  2:909319/1‑100 (MQ=255)

                                                                            tttCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGt                        >  2:1325640/1‑100 (MQ=255)

                                                                            tttCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGt                        <  1:2900257/100‑1 (MQ=255)

                                                                              tCCCGGTGAGTAATTTTGCCTCGTCAAAATTAGGAAACAGATAATCGTCATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGttt                      <  2:697238/100‑1 (MQ=255)

                                                                                ccGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGa                    >  2:2592658/1‑100 (MQ=255)

                                                                                ccGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGa                    >  2:1696198/1‑100 (MQ=255)

                                                                                 cGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGAt                   <  1:2868097/100‑1 (MQ=255)

                                                                                  ggTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATc                  >  1:893590/1‑100 (MQ=255)

                                                                                  ggTGAGTAATCTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATTTCATCCAGCGTTTCATTCAACCGCGGTTTGATc                  <  2:2746832/100‑1 (MQ=255)

                                                                                  atctcctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATc                  <  1:1891547/94‑1 (MQ=255)

                                                                                     tcctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCata               >  1:2525725/4‑100 (MQ=255)

                                                                                      cctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCatat              >  1:1587223/3‑100 (MQ=255)

                                                                                      cctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCata               >  1:671318/3‑99 (MQ=255)

                                                                                      cctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCagat              >  2:1891544/3‑97 (MQ=255)

                                                                                      cctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGAAATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACTTCGGTGTGATCatat              >  2:3190021/3‑100 (MQ=255)

                                                                                       ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATc             >  1:1067188/2‑100 (MQ=255)

                                                                                       ctAATTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATACAGCGTTTCATTCAACCGCGGTTTAATCATATc             >  2:620392/2‑100 (MQ=255)

                                                                                         aaTTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCgg           <  1:2376472/100‑1 (MQ=255)

                                                                                         aaTTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCgg           >  1:4154029/1‑100 (MQ=255)

                                                                                         aaTTTTGCCGAGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCgg           >  1:2688380/1‑100 (MQ=255)

                                                                                                ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGAt    >  1:3485661/1‑100 (MQ=255)

                                                                                                ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGAt    <  1:4136194/100‑1 (MQ=255)

                                                                                                ccTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGAt    >  2:2141258/1‑100 (MQ=255)

                                                                                                  tCGGCAAAATTAGAAAACAGATAATCGACATAACTCAAGGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATaa  >  2:2590819/1‑100 (MQ=255)

CGTCTTTACCCGTTTTAATCACCACCGTTTTTACGCCGCACGCAAGAAAGCAGTCAGCAATTTCATCCAGTGTCTCTTTCCCGGTGAGTAATTTTGCCTCGGCAAAATTAGGAAACAGATAATCGACATAACTCAATGCTTCGCAAATATCATCCAGCGTTTCATTCAACCGCGGTTTGATCATATCGGCACAGATAA  >  NC_000913/1856235‑1856432

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