BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,147,430	A→G	W139R (TGG→CGG)	yijP ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,147,430	0	A	G	84.8%	690.8 / 115.6	303	W139R (TGG→CGG)	yijP	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (22/24); new base G (118/139); total (140/163)
Fisher's exact test for biased strand distribution p-value = 8.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

Only 100 of 308 total aligned reads displayed.

CGATCGGATGCAGAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCACAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCTGGCTTCGTTGGTG  >  NC_000913/4147340‑4147521

 gATCGGATGCAGAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGcag                                                                                          >  2:1087981/1‑93 (MQ=255)

 gATCGGATGCAGAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCAcag                                                                                          >  2:1542026/1‑93 (MQ=255)

   tCGGATGCAGAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCgcagca                                                                                        >  1:6676955/1‑93 (MQ=255)

     ggATGCAGAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGa                                                                                      >  1:1697476/1‑93 (MQ=255)

      gATGCAGAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGaa                                                                                     <  2:6787580/93‑1 (MQ=255)

      gATGCAGAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGaa                                                                                     >  2:262087/1‑93 (MQ=255)

      gATGCAGAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGaa                                                                                     <  2:3076744/93‑1 (MQ=255)

      gATGCAGAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCACAGCAGaa                                                                                     <  2:7403791/93‑1 (MQ=255)

          cAGAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGc                                                                                 >  2:1699549/1‑93 (MQ=255)

            gAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCa                                                                               >  1:2281191/1‑93 (MQ=255)

            gAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCa                                                                               >  1:5889790/1‑93 (MQ=255)

            gAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCa                                                                               >  2:6217084/1‑93 (MQ=255)

             aaTCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCAc                                                                              <  2:4710494/93‑1 (MQ=255)

             aaTCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCAc                                                                              >  2:2837171/1‑93 (MQ=255)

              aTCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCACAGCAGAACTGCCAcc                                                                             <  1:3228622/93‑1 (MQ=255)

               tCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCAccg                                                                            >  2:1894394/1‑93 (MQ=255)

                 aaGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCAccgcc                                                                          >  1:5197731/1‑93 (MQ=255)

                    ccATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGta                                                                       >  2:511290/1‑93 (MQ=255)

                    ccATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCACAGCAGAACTGCCACCGCCGta                                                                       >  1:7192165/1‑93 (MQ=255)

                      aTAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGtata                                                                     <  2:7314100/93‑1 (MQ=255)

                      aTAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGtata                                                                     <  2:6996752/93‑1 (MQ=255)

                      aTAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCACAGCAGAACTGCCACCGCCGtata                                                                     >  1:3039461/1‑93 (MQ=255)

                       tAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAg                                                                    >  1:6435808/1‑93 (MQ=255)

                        aaaGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAgg                                                                   >  1:2438903/1‑93 (MQ=255)

                          aGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGcc                                                                 >  1:3258053/1‑93 (MQ=255)

                           gCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCa                                                                >  1:5568636/1‑93 (MQ=255)

                              gggCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAgcg                                                             >  1:284913/1‑93 (MQ=255)

                               ggCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAgcgc                                                            <  2:4357435/93‑1 (MQ=255)

                                  aaaaGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGAt                                                         <  1:3819423/93‑1 (MQ=255)

                                  aaaaGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGAt                                                         >  1:6057442/1‑93 (MQ=255)

                                    aaGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATaa                                                       <  1:4680703/93‑1 (MQ=255)

                                    aaGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATaa                                                       <  1:2703990/93‑1 (MQ=255)

                                      gagaCAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCACAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGc                                                     >  1:2202120/1‑93 (MQ=255)

                                         acaacaTAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGcaca                                                  <  2:7211883/93‑1 (MQ=255)

                                          caacaTAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACaa                                                 >  1:108524/1‑93 (MQ=255)

                                           aacaTAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAAt                                                <  1:673197/93‑1 (MQ=255)

                                           aacaTAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAAt                                                <  1:117996/93‑1 (MQ=255)

                                           aacaTAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAAt                                                <  2:6592231/93‑1 (MQ=255)

                                            acaTAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCACAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAAtt                                               >  1:250319/1‑93 (MQ=255)

                                             caTAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAAttt                                              >  1:3757861/1‑93 (MQ=255)

                                             caTAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAAttt                                              >  1:1774721/1‑93 (MQ=255)

                                              aTAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAAtttt                                             >  1:1335820/1‑93 (MQ=255)

                                               tAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTc                                            >  1:2381256/1‑93 (MQ=255)

                                               tAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTc                                            >  2:6363389/1‑93 (MQ=255)

                                                aaCGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCa                                           <  2:2740387/93‑1 (MQ=255)

                                                aaCGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCACAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCa                                           <  1:1438809/93‑1 (MQ=255)

                                                 aCGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAg                                          <  1:7529154/93‑1 (MQ=255)

                                                   gCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGc                                        >  2:3771170/1‑93 (MQ=255)

                                                    ccACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCt                                       <  2:3025436/93‑1 (MQ=255)

                                                      aCGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGa                                     <  1:3434322/93‑1 (MQ=255)

                                                      aCGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGa                                     >  1:2936561/1‑93 (MQ=255)

                                                        ggCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGaaa                                   >  2:7514873/1‑93 (MQ=255)

                                                         gCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAAt                                  <  2:3654255/93‑1 (MQ=255)

                                                         gCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAAt                                  >  2:3280958/1‑93 (MQ=255)

                                                         gCTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAAt                                  >  1:4497422/1‑93 (MQ=255)

                                                          cTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATa                                 <  1:4155882/93‑1 (MQ=255)

                                                          cTTTGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATa                                 <  1:5334677/93‑1 (MQ=255)

                                                            ttGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACt                               <  2:807931/93‑1 (MQ=255)

                                                             tGGAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTg                              <  1:578188/93‑1 (MQ=255)

                                                              ggAATATAGACCGGGCGCAGGCGTGTCCACAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTgg                             <  1:5836957/93‑1 (MQ=255)

                                                               gAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGc                            <  2:2458876/93‑1 (MQ=255)

                                                               gAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGc                            <  2:133265/93‑1 (MQ=255)

                                                               gAATATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGAAAAGCACAATTTTCAGGCTGAAATACTGGc                            >  2:3453012/1‑93 (MQ=255)

                                                                aaTATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCt                           >  2:2933062/1‑93 (MQ=255)

                                                                aaTATAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCt                           <  2:2256982/93‑1 (MQ=255)

                                                                aaTATAGACCGGGCGCAGGCGTGTCCACAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCt                           <  1:6559969/93‑1 (MQ=255)

                                                                 atatAGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCtt                          >  2:6194059/1‑93 (MQ=255)

                                                                 atatAGACCGGGCGCAGGCGTGTCCACAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCtt                          >  2:2979723/1‑93 (MQ=255)

                                                                  tataGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCttt                         <  2:7247577/93‑2 (MQ=255)

                                                                  tataGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTa                         <  2:498860/93‑1 (MQ=255)

                                                                  tataGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTa                         >  2:5276590/1‑93 (MQ=255)

                                                                  tataGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTa                         >  2:2211134/1‑93 (MQ=255)

                                                                  tataGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTa                         >  2:3504457/1‑93 (MQ=255)

                                                                    taGCCCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTaaa                       <  2:7144562/93‑1 (MQ=255)

                                                                    taGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTaaa                       <  1:6966157/93‑1 (MQ=255)

                                                                    taGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTaaa                       <  1:2348803/93‑1 (MQ=255)

                                                                    taGACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTaaa                       <  1:869491/93‑1 (MQ=255)

                                                                      gACCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATa                     <  2:1628197/93‑1 (MQ=255)

                                                                       aCCGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATAc                    <  2:4210638/93‑1 (MQ=255)

                                                                        ccGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACt                   <  1:1878543/93‑1 (MQ=255)

                                                                        ccGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACt                   <  2:5362218/93‑1 (MQ=255)

                                                                         cGGGCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTc                  <  1:3184437/93‑1 (MQ=255)

                                                                           ggCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGc                <  1:6231357/93‑1 (MQ=255)

                                                                           ggCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGc                >  1:3630372/1‑93 (MQ=255)

                                                                           ggCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGc                >  2:2135/1‑93 (MQ=255)

                                                                           ggCGCAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGc                >  2:4115769/1‑93 (MQ=255)

                                                                            gcgcAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCt               <  2:4897977/93‑1 (MQ=255)

                                                                            gcgcAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCt               <  2:6479092/93‑1 (MQ=255)

                                                                            gcgcAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCt               <  1:5457427/93‑1 (MQ=255)

                                                                             cgcAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCTg              <  2:3842132/93‑1 (MQ=255)

                                                                             cgcAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCTg              <  2:5574437/93‑1 (MQ=255)

                                                                             cgcAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCTg              <  1:6838918/93‑1 (MQ=255)

                                                                               cAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCTGGc            >  1:730684/1‑93 (MQ=255)

                                                                               cAGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCTGGc            <  2:5058813/93‑1 (MQ=255)

                                                                               cAGGCGTGTCCACAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCTGGc            <  1:3929716/93‑1 (MQ=255)

                                                                                aGGCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCTGGCt           <  1:4362731/93‑1 (MQ=255)

                                                                                 ggCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCTGGCtt          >  1:4896151/1‑93 (MQ=255)

                                                                                 ggCGTGTCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCTGGCtt          >  2:3373623/1‑93 (MQ=255)

                                                                                     tgtCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCTGGCTTCgtt      >  2:1418494/1‑93 (MQ=255)

                                                                                       tCCGCAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCTGGCTTCgttgg    >  2:5906863/1‑93 (MQ=255)

CGATCGGATGCAGAATCAAGCCATAAAGCAGGGCAAAAGAGACAACATAACGCCACGGCTTTGGAATATAGACCGGGCGCAGGCGTGTCCACAGCAGAACTGCCACCGCCGTATAGGCCAGCGCGATAAGCACAATTTTCAGGCTGAAATACTGGCTTAAATACTCGCTGGCTTCGTTGGTG  >  NC_000913/4147340‑4147521

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