BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,682,619	A→G	R199R (CGT→CGC)	yhjK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,682,619	0	A	G	82.5%	672.6 / 146.1	319	R199R (CGT→CGC)	yhjK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (27/29); new base G (130/133); total (157/162)
Fisher's exact test for biased strand distribution p-value = 8.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

Only 100 of 323 total aligned reads displayed.

TTTTCGTTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCAAAATCAGGC  >  NC_000913/3682527‑3682710

ttttCGTTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGgaga                                                                                             <  1:6967963/93‑1 (MQ=255)

  ttCGTTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGcg                                                                                           <  1:4786054/93‑1 (MQ=255)

    cGTTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGcgcg                                                                                         <  1:3382126/93‑1 (MQ=255)

    cGTTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGcgcg                                                                                         >  1:6885873/1‑93 (MQ=255)

      ttctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGt                                                                                       >  2:854502/1‑93 (MQ=255)

        ctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTaa                                                                                     <  2:7132035/93‑1 (MQ=255)

        ctgttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTaa                                                                                     <  1:5282224/93‑1 (MQ=255)

          gttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATg                                                                                   >  2:1233517/1‑93 (MQ=255)

          gttctTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATg                                                                                   >  2:4497123/1‑93 (MQ=255)

           ttcttcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGc                                                                                  <  2:6865848/93‑1 (MQ=255)

           ttcttcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGc                                                                                  >  1:4607600/1‑93 (MQ=255)

            tcttcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGcc                                                                                 >  1:7232422/1‑93 (MQ=255)

            tcttcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGcc                                                                                 <  1:828944/93‑1 (MQ=255)

             cttcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCa                                                                                >  2:3618450/1‑93 (MQ=255)

             cttcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCa                                                                                <  1:4005796/93‑1 (MQ=255)

             cttcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCa                                                                                >  1:5497333/1‑93 (MQ=255)

              ttcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAg                                                                               >  1:5197940/1‑93 (MQ=255)

              ttcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAg                                                                               >  1:2411027/1‑93 (MQ=255)

               tcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGt                                                                              <  2:455236/93‑1 (MQ=255)

               tcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGt                                                                              >  2:4705466/1‑93 (MQ=255)

               tcATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGt                                                                              <  2:636303/93‑1 (MQ=255)

                 aTAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAACGCCAGTTg                                                                            <  1:2567743/93‑1 (MQ=255)

                  tAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTgg                                                                           >  1:3017032/1‑93 (MQ=255)

                    aTGGCGCTGCAGCATTTGCTGGTGAAGGTTGAAACTGCGCACCAAAATACCGTTTTCACCGTCCGGAGGCAGGCGCGGTAATGCCAGTTGGTg                                                                         <  1:6154859/93‑1 (MQ=255)

                       gcgcTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGAcc                                                                      >  2:1970142/1‑93 (MQ=255)

                         gcTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACcaa                                                                    <  2:1542606/93‑1 (MQ=255)

                          cTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACcaac                                                                   <  1:3245266/93‑1 (MQ=255)

                          cTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACcaac                                                                   >  1:2851707/1‑93 (MQ=255)

                          cTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACcaac                                                                   >  1:5406405/1‑93 (MQ=255)

                           tGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACcaaca                                                                  <  1:3742273/93‑1 (MQ=255)

                            gcagcaATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGATAATGCCAGTTGGTGACcaacaa                                                                 >  2:6971745/1‑93 (MQ=255)

                              agcaATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGc                                                               <  1:2713922/93‑1 (MQ=255)

                                caATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTc                                                             <  1:1249706/93‑1 (MQ=255)

                                 aaTTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTcc                                                            <  2:3454000/93‑1 (MQ=255)

                                 aaTTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTcc                                                            >  2:2254352/1‑93 (MQ=255)

                                 aaTTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTcc                                                            >  2:1143902/1‑93 (MQ=255)

                                  aTTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCt                                                           <  1:2146016/93‑1 (MQ=255)

                                  aTTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCt                                                           >  1:5920957/1‑93 (MQ=255)

                                  aTTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCt                                                           <  2:4120815/93‑1 (MQ=255)

                                  aTTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCt                                                           <  2:2380455/93‑1 (MQ=255)

                                   ttGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCtt                                                          >  1:4365638/1‑93 (MQ=255)

                                    tGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTg                                                         >  2:4445442/1‑93 (MQ=255)

                                     gCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTgg                                                        >  2:7055748/1‑93 (MQ=255)

                                     gCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTgg                                                        >  2:4954860/1‑93 (MQ=255)

                                     gCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTgg                                                        <  1:1947684/93‑1 (MQ=255)

                                       tGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCt                                                      >  2:2970654/1‑93 (MQ=255)

                                         gTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTgg                                                    <  1:4561702/93‑1 (MQ=255)

                                         gTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTgg                                                    >  2:1377250/1‑93 (MQ=255)

                                          ttGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTggg                                                   >  1:3686792/1‑93 (MQ=255)

                                           tGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGa                                                  <  2:3079019/93‑1 (MQ=255)

                                           tGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGa                                                  <  2:624379/93‑1 (MQ=255)

                                            gAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGAt                                                 <  2:5580169/93‑1 (MQ=255)

                                            gAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGAt                                                 >  1:3077518/1‑93 (MQ=255)

                                             aGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATg                                                <  2:5117746/93‑1 (MQ=255)

                                             aGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATg                                                <  1:6582485/93‑1 (MQ=255)

                                             aGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATg                                                <  1:173962/93‑1 (MQ=255)

                                               gTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGc                                              <  2:504185/93‑1 (MQ=255)

                                               gTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGc                                              <  2:2525938/93‑1 (MQ=255)

                                                ttGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCg                                             >  1:6698807/1‑93 (MQ=255)

                                                ttGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCg                                             >  2:42782/1‑93 (MQ=255)

                                                ttGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCg                                             <  1:1322671/93‑1 (MQ=255)

                                                 tGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCGTTGGCTGGGATGGCGt                                            <  2:5802989/93‑1 (MQ=255)

                                                   tAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTa                                          <  2:5648713/93‑1 (MQ=255)

                                                   tAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGCGACCAACAAGCTCCTTGGCTGGGATGGCGTTa                                          <  1:115150/93‑1 (MQ=255)

                                                    aaCTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTaa                                         <  2:5385827/93‑1 (MQ=255)

                                                     aCTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAg                                        <  1:1390306/93‑1 (MQ=255)

                                                       tGCGCACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGtt                                      <  1:3527058/93‑1 (MQ=255)

                                                          gcACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTcgc                                   <  2:2744069/93‑1 (MQ=255)

                                                          gcACCAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTcgc                                   >  1:5670755/1‑93 (MQ=255)

                                                          gcACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTcgc                                   <  2:725081/93‑1 (MQ=255)

                                                             ccAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAg                                >  2:5967774/1‑93 (MQ=255)

                                                             ccAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAg                                <  1:3870660/93‑1 (MQ=255)

                                                              cAACATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGc                               >  2:6474517/1‑93 (MQ=255)

                                                              cAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGc                               >  1:3581965/1‑93 (MQ=255)

                                                               aaCATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCa                              <  2:6767124/93‑1 (MQ=255)

                                                                aCATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCaa                             >  1:2251757/1‑93 (MQ=255)

                                                                aCATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCaa                             <  2:5754292/93‑1 (MQ=255)

                                                                 cATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAat                            >  2:328476/1‑93 (MQ=255)

                                                                 cATACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAat                            <  1:6065572/93‑1 (MQ=255)

                                                                  aTACCGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAata                           <  1:5815943/93‑1 (MQ=255)

                                                                   tACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAatat                          >  1:6423604/1‑93 (MQ=255)

                                                                     ccGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTg                        >  2:5486868/1‑93 (MQ=255)

                                                                      cGATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGc                       >  1:874854/1‑93 (MQ=255)

                                                                       gATTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCg                      >  1:5098191/1‑93 (MQ=255)

                                                                        aTTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGt                     >  2:566471/1‑93 (MQ=255)

                                                                        aTTTCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGt                     >  2:7475485/1‑93 (MQ=255)

                                                                         tttCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTa                    >  1:2641170/1‑93 (MQ=255)

                                                                         tttCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTa                    >  2:5702821/1‑93 (MQ=255)

                                                                         tttCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTa                    <  2:1858073/93‑1 (MQ=255)

                                                                           tCATCGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAAc                  <  1:3363608/93‑1 (MQ=255)

                                                                              tCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGa               <  1:7419761/93‑1 (MQ=255)

                                                                               cGTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGAt              <  2:244058/93‑1 (MQ=255)

                                                                                gTCCTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATg             >  2:5901516/1‑93 (MQ=255)

                                                                                  ccTGATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCa           <  1:4315720/93‑1 (MQ=255)

                                                                                     gATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCaann        >  1:1808233/1‑91 (MQ=255)

                                                                                     gATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCaaaa        <  1:3178859/93‑1 (MQ=255)

                                                                                     gATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCaaaa        <  1:3495145/93‑1 (MQ=255)

                                                                                     gATGCAGGCGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCaaaa        >  2:4312337/1‑93 (MQ=255)

                                                                                       tGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCAAAATc      <  2:1729907/93‑1 (MQ=255)

                                                                                           gACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCAAAATCAGGc  <  2:4026335/93‑1 (MQ=255)

TTTTCGTTCTGTTCTTCATAATGGCGCTGCAGCAATTGCTGGTTGAGGTTGTAACTGCGCACCAACATACCGATTTCATCGTCCTGATGCAGACGCGGTAATGCCAGTTGGTGACCAACAAGCTCCTTGGCTGGGATGGCGTTAAGTTCGCGAGCAATATTGCGTAACGGATGCAAAATCAGGC  >  NC_000913/3682527‑3682710

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