BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,820,066	A→G	I72T (ATC→ACC)	osmE ←	OsmE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,820,066	0	A	G	85.6%	614.5 / 86.8	236	I72T (ATC→ACC)	osmE	OsmE
Reads supporting (aligned to +/- strand): ref base A (18/16); new base G (109/93); total (127/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 241 total aligned reads displayed.

AGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCT  >  NC_000913/1819975‑1820158

aGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGAt                                                                                             >  2:4477982/1‑93 (MQ=255)

 gCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTg                                                                                            >  1:4395469/1‑93 (MQ=255)

 gCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTg                                                                                            <  2:4819111/93‑1 (MQ=255)

 gCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATg                                                                                            >  1:2124015/1‑93 (MQ=255)

   acaGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTa                                                                                          <  1:5398183/93‑1 (MQ=255)

   acaGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTa                                                                                          <  1:3645556/93‑1 (MQ=255)

    caGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAg                                                                                         <  2:2004849/93‑1 (MQ=255)

    caGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAg                                                                                         <  2:4844609/93‑1 (MQ=255)

     aGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAgg                                                                                        <  2:5142776/93‑1 (MQ=255)

        tCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCt                                                                                     >  2:522701/1‑93 (MQ=255)

             tAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCaa                                                                                <  2:1705854/93‑1 (MQ=255)

             tAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCaa                                                                                <  2:733981/93‑1 (MQ=255)

              aaCCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAg                                                                               >  1:1116912/1‑93 (MQ=255)

              aaCCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAg                                                                               <  2:2455987/93‑1 (MQ=255)

               aCCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAga                                                                              >  2:584334/1‑92 (MQ=255)

                 cGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTAc                                                                            <  1:5038236/93‑1 (MQ=255)

                  ggAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTAcc                                                                           >  2:3721683/1‑93 (MQ=255)

                   gAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACcg                                                                          >  1:1385702/1‑93 (MQ=255)

                     gTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACcgcg                                                                        >  2:6149555/1‑93 (MQ=255)

                      ttGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGa                                                                       <  1:4694455/93‑1 (MQ=255)

                      ttGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGa                                                                       <  1:3264674/93‑1 (MQ=255)

                        gatgaCATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGc                                                                     <  1:5614891/93‑1 (MQ=255)

                         atgaCATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCa                                                                    <  2:717347/93‑1 (MQ=255)

                         atgaCATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCa                                                                    <  2:85451/93‑1 (MQ=255)

                          tgaCATGTCCGGTATCATCTAACGCGACAAAGTAGGGTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCAt                                                                   >  2:227911/1‑93 (MQ=255)

                           gaCATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATg                                                                  >  2:133440/1‑93 (MQ=255)

                           gaCATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATg                                                                  <  2:2728155/93‑1 (MQ=255)

                            aCATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATgg                                                                 >  1:4213591/1‑93 (MQ=255)

                             cATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGa                                                                >  1:5015799/1‑93 (MQ=255)

                             cATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGa                                                                >  1:352326/1‑93 (MQ=255)

                             cATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGa                                                                <  1:129125/93‑1 (MQ=255)

                              aTGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGAt                                                               >  1:99678/1‑93 (MQ=255)

                              aTGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGAt                                                               <  1:280816/93‑1 (MQ=255)

                                gTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCa                                                             <  1:3007431/93‑1 (MQ=255)

                                gTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCa                                                             <  1:1527417/93‑1 (MQ=255)

                                gTCCGCTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCa                                                             <  2:35289/93‑1 (MQ=255)

                                 tCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCAt                                                            >  1:6244510/1‑93 (MQ=255)

                                   cGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGc                                                          <  1:4043153/93‑1 (MQ=255)

                                   cGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGc                                                          >  2:1171017/1‑93 (MQ=255)

                                    ggTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCt                                                         >  2:2682368/1‑93 (MQ=255)

                                    ggTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCt                                                         >  1:578123/1‑93 (MQ=255)

                                      gatcagcTAACGCGACAAAGGAGGTTTCTGCTTTACCATCGCGGTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCa                                                       <  2:3251611/92‑1 (MQ=255)

                                       atcatcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCAc                                                      >  1:2620309/1‑93 (MQ=255)

                                        tcatcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACt                                                     >  2:1415916/1‑93 (MQ=255)

                                        tcatcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACt                                                     >  2:1019277/1‑93 (MQ=255)

                                        tcatcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACt                                                     >  1:390667/1‑93 (MQ=255)

                                        tcatcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACt                                                     >  2:3350450/1‑93 (MQ=255)

                                         catcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACtt                                                    <  2:1956399/93‑1 (MQ=255)

                                         catcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACtt                                                    >  1:6282496/1‑93 (MQ=255)

                                          atcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTc                                                   >  2:5619932/1‑93 (MQ=255)

                                           tcTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGAAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCa                                                  >  2:1074049/1‑93 (MQ=255)

                                            cTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAg                                                 >  2:4165367/1‑93 (MQ=255)

                                             tAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGa                                                <  1:277361/93‑1 (MQ=255)

                                             tAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGa                                                >  2:408422/1‑93 (MQ=255)

                                              aaCGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGAc                                               <  2:5444719/93‑1 (MQ=255)

                                              aaCGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGAc                                               >  2:1475421/1‑93 (MQ=255)

                                               aCGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACg                                              <  1:5859845/93‑1 (MQ=255)

                                                 gcgACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGaa                                            <  2:3591876/93‑1 (MQ=255)

                                                  cgACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTCCCGCGAGCATGGATCATGCTCACTTCAGACGAAg                                           <  1:1893827/93‑1 (MQ=255)

                                                    aCAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGt                                         >  1:4824957/1‑93 (MQ=255)

                                                     cAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGtt                                        <  1:5688739/93‑1 (MQ=255)

                                                     cAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGtt                                        >  1:256707/1‑93 (MQ=255)

                                                       aaGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTa                                      >  2:3947319/1‑93 (MQ=255)

                                                       aaGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCACGAGCATGTATCATGCTCACTTCAGACGAAGGTTTa                                      >  2:2903738/1‑93 (MQ=255)

                                                        aGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTAc                                     <  2:1548359/93‑1 (MQ=255)

                                                        aGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTAc                                     >  1:1768914/1‑93 (MQ=255)

                                                        aGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTAc                                     >  1:2136599/1‑93 (MQ=255)

                                                         gTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTAcc                                    >  1:2409986/1‑93 (MQ=255)

                                                         gTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTAcc                                    >  1:41277/1‑93 (MQ=255)

                                                         gTAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTAcc                                    <  1:540377/93‑1 (MQ=255)

                                                          tAGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTAccc                                   <  2:5168030/93‑1 (MQ=255)

                                                           aGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCg                                  <  2:1828519/93‑1 (MQ=255)

                                                           aGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCg                                  >  1:5594661/1‑93 (MQ=255)

                                                           aGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCg                                  <  2:4318926/93‑1 (MQ=255)

                                                           aGGTTTCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCg                                  >  1:5258455/1‑93 (MQ=255)

                                                              tttCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCat                               >  1:6154233/1‑92 (MQ=255)

                                                              tttCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCaa                               <  2:5467924/93‑1 (MQ=255)

                                                                tCTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAAtt                             >  1:1944762/1‑93 (MQ=255)

                                                                tCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCTGCAAtt                             >  2:3552469/1‑93 (MQ=255)

                                                                 cTGCTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAAttt                            >  1:2830926/1‑93 (MQ=255)

                                                                    cTTTACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTtgtg                         >  1:3411380/1‑93 (MQ=255)

                                                                       tACCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCaa                      >  1:788415/1‑93 (MQ=255)

                                                                        aCCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAAc                     >  1:5131392/1‑93 (MQ=255)

                                                                        aCCATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAAc                     >  2:5250791/1‑93 (MQ=255)

                                                                         ccATCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAAcc                    <  2:2461578/93‑1 (MQ=255)

                                                                           aTCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTg                  >  1:6610473/1‑93 (MQ=255)

                                                                            tCACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTgc                 <  1:266939/93‑1 (MQ=255)

                                                                             cACGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTgcg                >  1:1554175/1‑93 (MQ=255)

                                                                               cGTTGACCCAGGGTGTAGGTCTGGCACGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGcc              <  2:6644538/93‑1 (MQ=255)

                                                                               cGTTGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGcc              <  1:3593616/93‑1 (MQ=255)

                                                                                 ttGACCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCg            <  2:6580928/93‑1 (MQ=255)

                                                                                    aCCCAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCt         >  1:6041212/1‑93 (MQ=255)

                                                                                    aCCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCt         >  1:6506407/1‑93 (MQ=255)

                                                                                     cccAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTc        >  2:923870/1‑93 (MQ=255)

                                                                                      ccAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCa       >  1:3082139/1‑93 (MQ=255)

                                                                                      ccAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCa       >  2:1104159/1‑93 (MQ=255)

                                                                                       cAGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCAt      >  2:499925/1‑93 (MQ=255)

                                                                                        aGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATg     <  1:2475503/93‑1 (MQ=255)

                                                                                        aGGGTGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATg     <  2:3305924/93‑1 (MQ=255)

                                                                                          gATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGcc   >  2:5315503/1‑93 (MQ=255)

AGCACAGGTCTGATAACCGGAGTTGATGACATGTCCGGTATCATCTAACGCGACAAAGTAGGTTTCTGCTTTACCATCACGTTGACCCAGGATGTAGGTCTGGCAAGTACCGCGAGCATGGATCATGCTCACTTCAGACGAAGGTTTACCCGCAATTTGTGCAACCTGCGCCCGGCTCATGCCT  >  NC_000913/1819975‑1820158

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