BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,084,593	A→T	5.6%	intergenic (‑67/+5)	PP_t54 ← / ← PP_t53	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,084,593	0	A	T	5.6%	799.3 / 17.8	245	intergenic (‑67/+5)	PP_t54/PP_t53	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base A (142/87); new base T (7/9); total (149/96)
Fisher's exact test for biased strand distribution p-value = 1.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 250 total aligned reads displayed.

AAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATAAAACAACTAC  >  NC_002947/5084447‑5084742

ctgTGGCGTACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCaaaa                                                                                                                                                    <  1:3331839/148‑1 (MQ=255)

    ggCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTgg                                                                                                                                                <  1:2634697/150‑1 (MQ=255)

     gCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGc                                                                                                                                               <  1:1245899/150‑1 (MQ=255)

     gCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGc                                                                                                                                               <  1:6814296/150‑1 (MQ=255)

     gCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGc                                                                                                                                               <  1:4527389/150‑1 (MQ=255)

     gCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGc                                                                                                                                               <  1:1894367/150‑1 (MQ=255)

     gCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTCGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGc                                                                                                                                               <  2:3508823/150‑1 (MQ=255)

        caCTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGa                                                                                                                                            <  1:6403167/149‑1 (MQ=255)

             agcagGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTGGAGTGAAGCCGCCTTGCAACCAAAGCCACTTCAACAGGGGCgggcaca                                                                                                                                        >  1:6659650/1‑145 (MQ=11)

             agcagGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTcag                                                                                                                                       <  1:6221878/150‑1 (MQ=255)

              gcagGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTcagc                                                                                                                                      >  1:2201929/1‑150 (MQ=255)

              gcagGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTcagc                                                                                                                                      <  1:1437543/150‑1 (MQ=255)

              gcagGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTcagc                                                                                                                                      <  1:6961190/150‑1 (MQ=255)

              gcagGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACAAGA‑TACCGCTGGATGATTGAAGCAGCCTTGCAAACAAAGCAACTTCAAAAGTGGCGGACTcagc                                                                                                                                      >  2:7395331/1‑150 (MQ=255)

               cagGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGc                                                                                                                                               <  1:3590260/140‑1 (MQ=255)

                  gATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGa                                                                                                                                  <  1:3769022/150‑1 (MQ=255)

                       gAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGa                                                                                                                             >  2:374087/1‑150 (MQ=255)

                           cTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCt                                                                                                                         >  2:6825068/1‑150 (MQ=255)

                               gACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGa                                                                                                                     <  1:6743059/150‑1 (MQ=255)

                                 ccGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGAcc                                                                                                                   >  1:3284978/1‑150 (MQ=255)

                                  cGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCg                                                                                                                  >  2:6484972/1‑150 (MQ=255)

                                    cgcGGTTCGTAGCCGAGTACTCTAGCCAGCCCAGCTATGAGTCCGCAATTCGCAGTTCTCGCCCAGA‑TCCCGCTGGTTGATTGAAGCAGCCTTGCTAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCt                                                                                                                <  2:4563204/148‑1 (MQ=18)

                                    cTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCCGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGACCCTGCGACCGCt                                                                                                                >  2:3667167/1‑150 (MQ=255)

                                    cTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCt                                                                                                                <  2:4563800/150‑1 (MQ=255)

                                      cGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGCGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGATGCAGCCTTGCAAACAAAGCCACTACAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCg                                                                                                              <  1:4026980/150‑1 (MQ=255)

                                      cGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCg                                                                                                              <  1:5503500/150‑1 (MQ=255)

                                       ggTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCgg                                                                                                             >  1:1183193/1‑150 (MQ=255)

                                           agTAGCCGAGTACACGATCCAGCTGAGCTATGAGTCCGCACTGTGCAGTTTTAGAGCATA‑TACCGGTGTGAGATAGAAGCATCCTCTCTAACAAAGCCACCTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCg                                                                                                         <  2:4864342/149‑1 (MQ=18)

                                            gTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGt                                                                                                        >  1:2974818/1‑150 (MQ=255)

                                             tAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGa                                                                                                                     >  2:6204459/1‑136 (MQ=255)

                                              aGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAg                                                                                                      >  2:5367451/1‑150 (MQ=255)

                                                ccGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGcc                                                                                                    <  1:2121725/150‑1 (MQ=255)

                                                  gAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGa                                                                                                  >  2:165049/1‑150 (MQ=255)

                                                     tACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCATTTTTAGACCTGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTa                                                                                               <  2:2000463/150‑1 (MQ=255)

                                                      aCTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTAc                                                                                              >  2:1339702/1‑150 (MQ=255)

                                                       ctctATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTAct                                                                                             >  1:5171295/1‑150 (MQ=255)

                                                        tctATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTActc                                                                                            >  2:4204312/1‑150 (MQ=255)

                                                        tctATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACCCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTActc                                                                                            >  2:849317/1‑150 (MQ=255)

                                                          tATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                                                          >  1:615621/1‑150 (MQ=255)

                                                          tATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                                                          >  1:523497/1‑150 (MQ=255)

                                                          tATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                                                          >  2:548052/1‑150 (MQ=255)

                                                          tATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                                                          >  1:5676688/1‑150 (MQ=255)

                                                            tCCAGCTGAGCTATGCGTCCGCTAGTGGCAGTATTAGACCAGA‑TACCGCTGGTTGATAGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                                                        <  2:7275321/150‑1 (MQ=18)

                                                            tCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                                                        >  1:1286256/1‑150 (MQ=255)

                                                            tCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                                                        >  2:5556262/1‑150 (MQ=255)

                                                            tCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                                                        <  2:822497/150‑1 (MQ=255)

                                                              cAGCTGAGCTATGAGTCCGCAATTGGCATTTTTAGACCAGA‑TAGCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCa                                                                                      <  2:4404093/150‑1 (MQ=255)

                                                                  tGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTg                                                                                  >  2:6344731/1‑150 (MQ=255)

                                                                   gAGCTAGGAGTCCGCACTTGGCAGTTTTAGACCAGA‑TACCGCTGGGTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                                                 <  2:4072596/150‑1 (MQ=255)

                                                                       tATGAGTCCGCAATTGTCAGTTTTAGACCAGA‑GACCGTTGGTTGATTGACGCAGCCTTGTATAGAAAGCCTCGTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTa                                                                             <  2:4168013/150‑1 (MQ=18)

                                                                        aTGATTCCGCAATTGGCACCTTTAGACCTTT‑TAGCTGTGGTTGATTGAAGGAGCCTTTAAAACAAACCCACTTATAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                                            <  2:1166825/150‑1 (MQ=17)

                                                                          gAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                                          >  1:6992817/1‑150 (MQ=255)

                                                                             tCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAg                                                                                                 >  2:5038679/1‑124 (MQ=255)

                                                                             tCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAg                                                                                                 <  1:5038679/124‑1 (MQ=255)

                                                                             tCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTc                                                                       >  2:2581585/1‑150 (MQ=255)

                                                                                gCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCgt                                                                    >  2:3005438/1‑150 (MQ=255)

                                                                                gCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCgt                                                                    >  1:6510666/1‑150 (MQ=255)

                                                                                       gCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGc                                                             >  1:228240/1‑150 (MQ=255)

                                                                                             ttAGACCAGT‑TACCGATGGTTGATTGTAGGAGCGTTGCTAAGAAAGCCAGTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCa                                                       <  1:2106488/150‑1 (MQ=18)

                                                                                              tAGACCAGA‑GACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCAt                                                      <  2:5123826/150‑1 (MQ=255)

                                                                                                 aCCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATtat                                                   >  1:1380714/1‑150 (MQ=255)

                                                                                                   cAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAg                                                 >  2:3971292/1‑150 (MQ=255)

                                                                                                     gA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGt                                               <  2:2937366/150‑1 (MQ=255)

                                                                                                         aCCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGGGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGc                                            >  1:1028063/1‑150 (MQ=255)

                                                                                                           cGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTg                                          >  2:6653979/1‑150 (MQ=255)

                                                                                                            gCTGGTTGATTTAAGCAGCATTTAAAAGAAAGACACTTAAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATTAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGa                                         <  1:3674008/150‑1 (MQ=21)

                                                                                                            gCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGa                                         >  2:2208099/1‑150 (MQ=255)

                                                                                                              tGGTTGATTGAAGCATCCTCTCAAACAAAGCCACCTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGaaa                                       <  2:3585942/150‑1 (MQ=255)

                                                                                                                   gattgaAGCATCCTTTCAAACAAAGCCTCTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACAttt                                  <  2:7139669/150‑1 (MQ=255)

                                                                                                                   gattgaAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACAttt                                  >  1:4918059/1‑150 (MQ=255)

                                                                                                                     ttgaAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACAttttt                                >  2:4372169/1‑150 (MQ=255)

                                                                                                                      tgaAGCAGCCTTGCAAACAAAGCCACTTCTAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTg                               <  2:1001255/150‑1 (MQ=255)

                                                                                                                         agcagcCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCaa                            >  1:4315586/1‑150 (MQ=255)

                                                                                                                            agcCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCaagaa                         >  1:2785522/1‑150 (MQ=255)

                                                                                                                              ccTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGc                       >  1:7265545/1‑150 (MQ=255)

                                                                                                                                ttGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCtt                     >  1:5726524/1‑150 (MQ=255)

                                                                                                                                ttGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCtt                     >  1:5869626/1‑150 (MQ=255)

                                                                                                                                ttGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCtt                     >  1:7433257/1‑150 (MQ=255)

                                                                                                                                   caaacaaaGCCACGTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTAGAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTc                  <  2:5628977/150‑1 (MQ=255)

                                                                                                                                    caacaacGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCg                 <  2:4807173/149‑1 (MQ=255)

                                                                                                                                      acaaaGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGtt               >  1:588294/1‑150 (MQ=255)

                                                                                                                                       caaaGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTta              <  1:7034357/150‑1 (MQ=255)

                                                                                                                                        tccgCCACTTTTGAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTtat             >  2:5824733/4‑150 (MQ=255)

                                                                                                                                        caaGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTtat             >  2:4968466/2‑150 (MQ=255)

                                                                                                                                          cgCCACTTTTGAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATaa           >  2:3260943/2‑150 (MQ=255)

                                                                                                                                          aGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATaa           <  1:4757789/150‑1 (MQ=255)

                                                                                                                                           gCCACTTTTGAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATaaa          >  2:1336918/1‑150 (MQ=255)

                                                                                                                                           gCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATaaa          >  1:1623704/1‑150 (MQ=255)

                                                                                                                                             cACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACAttttt                                >  1:4869244/1‑126 (MQ=255)

                                                                                                                                             cACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATAAaac        >  1:4204204/1‑150 (MQ=255)

                                                                                                                                             cACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTAAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATAAaac        <  1:739037/150‑1 (MQ=255)

                                                                                                                                                ttCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACAtttt                                 <  1:7172512/122‑1 (MQ=255)

                                                                                                                                                ttCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACAtttt                                 >  2:7172512/1‑122 (MQ=255)

                                                                                                                                                ttCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTtat             <  1:7370850/142‑1 (MQ=255)

                                                                                                                                                ttCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTtat             >  2:7370850/1‑142 (MQ=255)

                                                                                                                                                ttCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATAAaacaag     >  2:5738255/1‑149 (MQ=255)

                                                                                                                                                ttCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATAAaacaac     >  1:6807153/1‑150 (MQ=255)

                                                                                                                                                ttCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATAAaacaac     >  2:1847628/1‑150 (MQ=255)

                                                                                                                                                   aaaaGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATAGAACAactac  >  2:5949188/1‑150 (MQ=255)

                                                                                                                                                   aaaaGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATAAAACAactac  >  1:2583186/1‑150 (MQ=255)

AAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGA‑TACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATAAAACAACTAC  >  NC_002947/5084447‑5084742

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