BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	291,960	T→C	R71R (AGA→AGG)	yagK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	291,960	0	T	C	99.3%	455.3 / ‑4.9	144	R71R (AGA→AGG)	yagK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (68/75); total (68/76)
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.05e-01

Only 100 of 149 total aligned reads displayed.

ATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCTCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAACTATTTTGGGATA  >  NC_000913/291870‑292049

aTACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGccct                                                                                         >  2:4047471/1‑93 (MQ=255)

aTACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGcact                                                                                         >  1:4082150/1‑93 (MQ=255)

 tcctctTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTg                                                                                        <  2:1127042/91‑1 (MQ=255)

  aCTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGa                                                                                       <  1:4561196/93‑1 (MQ=255)

  aCTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGa                                                                                       <  1:2855124/93‑1 (MQ=255)

   ctctTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAg                                                                                      <  2:82414/93‑1 (MQ=255)

   ctctTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAg                                                                                      >  2:3751855/1‑93 (MQ=255)

   ctctTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAg                                                                                      <  1:3335096/93‑1 (MQ=255)

    tctTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAgg                                                                                     >  1:1246314/1‑93 (MQ=255)

    tctTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAgg                                                                                     <  2:3708443/93‑1 (MQ=255)

    tctTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAgg                                                                                     <  2:2223678/93‑1 (MQ=255)

    tctTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAgg                                                                                     <  1:2471969/93‑1 (MQ=255)

     ctTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAggg                                                                                    >  2:840448/1‑93 (MQ=255)

      tttCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGa                                                                                   >  1:3205103/1‑93 (MQ=255)

      tttCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGa                                                                                   <  1:29219/93‑1 (MQ=255)

      tttCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGa                                                                                   >  1:1983364/1‑93 (MQ=255)

          gCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCa                                                                               >  2:460332/1‑93 (MQ=255)

           cccATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCAc                                                                              <  2:3118931/93‑1 (MQ=255)

            ccATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACg                                                                             >  1:749124/1‑93 (MQ=255)

              atatAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCa                                                                           >  1:1895903/1‑93 (MQ=255)

              atatAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCa                                                                           <  2:1664691/93‑1 (MQ=255)

               tataATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCAt                                                                          <  1:1287563/93‑1 (MQ=255)

               tataATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCAt                                                                          <  2:4276733/93‑1 (MQ=255)

               tataATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCAt                                                                          >  1:4471287/1‑93 (MQ=255)

               tataATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCAt                                                                          >  1:4424071/1‑93 (MQ=255)

                ataataAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATc                                                                         >  1:537521/1‑93 (MQ=255)

                ataataAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATc                                                                         <  1:4669297/93‑1 (MQ=255)

                ataataAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATc                                                                         <  1:583270/93‑1 (MQ=255)

                   ataaataaAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCgg                                                                      >  2:110304/1‑93 (MQ=255)

                   ataaataaAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGACCTGAGGGATTCACGCATCCgg                                                                      >  2:771884/1‑93 (MQ=255)

                    taaataaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCggg                                                                     >  2:1677372/1‑93 (MQ=255)

                    taaataaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCggg                                                                     <  2:4424891/93‑1 (MQ=255)

                    taaataaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCggg                                                                     >  2:665597/1‑93 (MQ=255)

                     aaataaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGa                                                                    <  2:2471540/93‑1 (MQ=255)

                     aaataaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGa                                                                    <  1:1911280/93‑1 (MQ=255)

                     aaataaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGa                                                                    <  1:2787823/93‑1 (MQ=255)

                      aataaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGat                                                                   >  2:3397410/1‑93 (MQ=255)

                      aataaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGat                                                                   >  1:1114072/1‑93 (MQ=255)

                       ataaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACTGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGata                                                                  <  1:2881274/93‑1 (MQ=255)

                       ataaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGata                                                                  <  1:3269151/93‑1 (MQ=255)

                        taaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGatat                                                                 <  1:2347491/93‑1 (MQ=255)

                        taaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGatat                                                                 <  1:2728832/93‑1 (MQ=255)

                        taaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGatat                                                                 >  2:3827075/1‑93 (MQ=255)

                        taaaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCAGGatat                                                                 <  1:2728831/93‑1 (MQ=255)

                          aaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCa                                                               >  2:4591800/1‑93 (MQ=255)

                          aaGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCa                                                               >  2:1829505/1‑93 (MQ=255)

                            gggCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCAcc                                                             >  2:1221132/1‑93 (MQ=255)

                              gCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCTCTGAGGGATTCACGCATCCGGGATATCAcccc                                                           <  1:398296/93‑1 (MQ=255)

                              gCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCAcccc                                                           <  1:2773375/93‑1 (MQ=255)

                                aaCGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCccag                                                         <  1:2846208/93‑1 (MQ=255)

                                aaCGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCccag                                                         >  1:2366962/1‑93 (MQ=255)

                                  cGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCccaggt                                                       <  2:2160715/93‑1 (MQ=255)

                                  cGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCccaggt                                                       >  2:4539098/1‑93 (MQ=255)

                                  cGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCccaggt                                                       <  2:3488444/93‑1 (MQ=255)

                                   ggTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCNTCCGGGATATCACCccaggtt                                                      >  1:2217807/1‑93 (MQ=255)

                                      aaaTGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCccaggttcca                                                   <  2:2234118/93‑1 (MQ=255)

                                       aaTGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCccaggttccag                                                  >  1:940031/1‑93 (MQ=255)

                                       aaTGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCccaggttccag                                                  <  2:4031141/93‑1 (MQ=255)

                                       aaTGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCccaggttccag                                                  <  1:3197406/93‑1 (MQ=255)

                                           cGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTa                                              <  1:1597085/93‑1 (MQ=255)

                                             tctATCCTCACGCCCCTTACGCCAACGGTCCGCCTCCAGCCTAGCTCTGAGCGATTCACGCATCCGGCATATCGCCCCAGGTTGCAGGTTAgg                                            <  2:3344667/91‑1 (MQ=255)

                                               tATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGa                                          <  2:3605798/93‑1 (MQ=255)

                                                aTCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGaa                                         <  1:3203915/93‑1 (MQ=255)

                                                aTCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGaa                                         <  2:1238796/93‑1 (MQ=255)

                                                aTCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGaa                                         >  2:2850892/1‑93 (MQ=255)

                                                  ccTCACGCCCCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGc                                       <  2:2937309/93‑1 (MQ=255)

                                                  ccTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGc                                       <  1:3854499/93‑1 (MQ=255)

                                                   cTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCa                                      >  1:4333687/1‑93 (MQ=255)

                                                   cTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCa                                      >  1:897191/1‑93 (MQ=255)

                                                    tCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCaa                                     <  1:1266682/93‑1 (MQ=255)

                                                    tCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCaa                                     >  2:1038856/1‑93 (MQ=255)

                                                       cGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAg                                  >  1:818754/1‑93 (MQ=255)

                                                        gCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGa                                 >  1:4560706/1‑93 (MQ=255)

                                                        gCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGa                                 <  1:1250774/93‑1 (MQ=255)

                                                         cACCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGat                                >  1:3179387/1‑93 (MQ=255)

                                                          aCCTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGata                               >  2:4465493/1‑93 (MQ=255)

                                                            cTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATAtt                             <  1:4092976/93‑1 (MQ=255)

                                                            cTTACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATAtt                             >  2:3651223/1‑93 (MQ=255)

                                                              tACGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGt                           >  2:4069949/1‑93 (MQ=255)

                                                                cGGGTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCa                         <  1:3217252/93‑1 (MQ=255)

                                                                 gggTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCAc                        <  2:2315460/93‑1 (MQ=255)

                                                                  ggTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCAcc                       >  2:2651174/1‑93 (MQ=255)

                                                                  ggTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCAcc                       <  2:1370801/93‑1 (MQ=255)

                                                                   gTACGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCg                      <  1:4042170/93‑1 (MQ=255)

                                                                      cGGTCAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTa                   <  1:1747504/93‑1 (MQ=255)

                                                                          cAGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGGTCCAGGTTAGGGAAGCGACAGATATTGTCACCGTTATCaa               >  1:4363890/1‑93 (MQ=255)

                                                                           aGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAAc              <  1:3799094/93‑1 (MQ=255)

                                                                           aGCCTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAAc              <  1:3773322/93‑1 (MQ=255)

                                                                           aGCCTCCAGCTTAGCCCTGAGCGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAAc              <  2:1433685/93‑1 (MQ=255)

                                                                             ccTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAACTa            <  1:3063385/93‑1 (MQ=255)

                                                                              cTCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAACTAt           <  1:1738968/93‑1 (MQ=255)

                                                                               tCCAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAACTAtt          >  1:2473189/1‑93 (MQ=255)

                                                                                ccAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCCACTAttt         >  2:1969463/1‑93 (MQ=255)

                                                                                ccAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAACTAttt         >  1:3898893/1‑93 (MQ=255)

                                                                                 cAGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAACTAtttt        <  2:1022129/93‑1 (MQ=255)

                                                                                  aGCTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTACGGAAGCAACAGATATTGTCACCGTTATCAACTATTTTg       >  2:807796/1‑93 (MQ=255)

                                                                                    cTTAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAACTATTTTggg     <  2:511817/93‑1 (MQ=255)

                                                                                     ttAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAACTATTTTGGGa    <  1:1554953/93‑1 (MQ=255)

                                                                                     ttAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAACTATTTTGGGa    <  2:2021385/93‑1 (MQ=255)

                                                                                      tAGCCCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAACTATTTTGGGAt   <  2:2581759/93‑1 (MQ=255)

ATACTCTTTCGCCCATATAATAAATAAAGGGCAACGGTAAATGCGTTTATCCTCACGCACCTTACGGGTACGGTCAGCCTCCAGCTTAGCTCTGAGGGATTCACGCATCCGGGATATCACCCCAGGTTCCAGGTTAGGGAAGCAACAGATATTGTCACCGTTATCAACTATTTTGGGATA  >  NC_000913/291870‑292049

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