BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,038,976	C→T	E39K (GAA→AAA)	ygfY ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,038,976	0	C	T	88.1%	858.4 / 75.4	278	E39K (GAA→AAA)	ygfY	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (16/17); new base T (118/127); total (134/144)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

Only 100 of 285 total aligned reads displayed.

GTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCAATGAATG  >  NC_000913/3038886‑3039067

  tCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCa                                                                                         <  1:1008956/93‑1 (MQ=255)

    tGCATCGGCCAGTTTACCGTGATTCATCAGCCAGTTAACCAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATc                                                                                       <  2:2488284/93‑1 (MQ=255)

       aTCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCt                                                                                    <  1:4281963/93‑1 (MQ=255)

       aTCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCt                                                                                    >  2:2373675/1‑93 (MQ=255)

       aTCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTCGTCATCGCt                                                                                    <  1:4053383/93‑1 (MQ=255)

         cGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTa                                                                                  >  1:1018117/1‑93 (MQ=255)

          ggCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTaa                                                                                 <  2:1390665/93‑1 (MQ=255)

           gCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAg                                                                                >  2:5261198/1‑93 (MQ=255)

           gCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAg                                                                                >  2:1905027/1‑93 (MQ=255)

           gCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAg                                                                                >  1:3054750/1‑93 (MQ=255)

           gCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTCGTCATCGCTTAAg                                                                                >  2:2123059/1‑93 (MQ=255)

             cGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCt                                                                              >  1:1201561/1‑93 (MQ=255)

             cGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCt                                                                              >  2:1100555/1‑93 (MQ=255)

             cGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCt                                                                              <  2:6356986/93‑1 (MQ=255)

               gTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGt                                                                            >  2:2594153/1‑93 (MQ=255)

                tttACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGTAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTc                                                                           <  1:4809448/93‑1 (MQ=255)

                tttACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTgac                                                                           >  1:57021/1‑91 (MQ=255)

                 ttACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCg                                                                          <  2:4157819/93‑1 (MQ=255)

                  tACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGt                                                                         <  1:6445342/93‑1 (MQ=255)

                   aCCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTa                                                                        >  2:1721887/1‑93 (MQ=255)

                    ccGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTAc                                                                       >  1:2488185/1‑93 (MQ=255)

                      gTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTc                                                                     <  1:1835540/93‑1 (MQ=255)

                       tGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCa                                                                    >  2:4411608/1‑93 (MQ=255)

                       ggATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCa                                                                    <  2:1642062/92‑1 (MQ=255)

                           tcatcaGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGtt                                                                <  2:4826106/93‑1 (MQ=255)

                           tcatcaGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGtt                                                                >  2:6450797/1‑93 (MQ=255)

                            catcaGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTc                                                               <  2:916559/93‑1 (MQ=255)

                             atcaGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCg                                                              >  2:1322764/1‑93 (MQ=255)

                                aGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGaaa                                                           <  1:5302265/93‑1 (MQ=255)

                                 gCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGaaaa                                                          <  1:2851408/93‑1 (MQ=255)

                                 gCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGaaaa                                                          <  1:6384834/93‑1 (MQ=255)

                                  ccAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGaaaaa                                                         >  2:3838285/1‑93 (MQ=255)

                                  ccAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGaaaaa                                                         <  1:4134497/93‑1 (MQ=255)

                                   cAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAAc                                                        <  2:5871684/93‑1 (MQ=255)

                                   cAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAAc                                                        >  2:5609539/1‑93 (MQ=255)

                                   cAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAAc                                                        <  2:2995300/93‑1 (MQ=255)

                                   cAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAAc                                                        >  1:658870/1‑93 (MQ=255)

                                     gTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACgg                                                      >  2:6481370/1‑93 (MQ=255)

                                      ttAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGc                                                     <  1:4378155/93‑1 (MQ=255)

                                       tAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCa                                                    <  2:3328859/93‑1 (MQ=255)

                                         aaCAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATg                                                  >  1:4528164/1‑93 (MQ=255)

                                            aGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCTTGTTCGAAAAACGGCATGAtt                                               <  1:2200170/93‑1 (MQ=255)

                                            aGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGAtt                                               >  1:3221524/1‑93 (MQ=255)

                                             ggTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTg                                              >  2:3337846/1‑93 (MQ=255)

                                             ggTCCGGATCGCCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTg                                              >  2:4545862/1‑93 (MQ=255)

                                              gTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGa                                             >  2:4943947/1‑93 (MQ=255)

                                              gTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGa                                             >  1:1637972/1‑93 (MQ=255)

                                               tCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGaa                                            >  1:1948371/1‑93 (MQ=255)

                                                 cGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAat                                          <  2:4572746/93‑1 (MQ=255)

                                                 cGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAat                                          >  1:3141899/1‑93 (MQ=255)

                                                  ggATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAata                                         <  1:3490169/93‑1 (MQ=255)

                                                   gATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAatat                                        <  2:6214545/93‑1 (MQ=255)

                                                    aTCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATc                                       <  1:1400177/93‑1 (MQ=255)

                                                     tCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCg                                      <  2:2779602/93‑1 (MQ=255)

                                                     tCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCg                                      >  1:3931451/1‑93 (MQ=255)

                                                      cGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGa                                     <  1:3358250/93‑1 (MQ=255)

                                                        tCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGt                                   <  2:6614642/93‑1 (MQ=255)

                                                        tCACATTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGt                                   <  1:2928708/93‑1 (MQ=255)

                                                         cacaTTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGtt                                  >  2:1485160/1‑93 (MQ=255)

                                                         cacaTTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGtt                                  >  1:823669/1‑93 (MQ=255)

                                                         cacaTTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGtt                                  <  2:4732805/93‑1 (MQ=255)

                                                         cacaTTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGtt                                  <  1:4583723/93‑1 (MQ=255)

                                                          acaTTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTc                                 <  1:5495762/93‑1 (MQ=255)

                                                           caTTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTcg                                >  1:5024481/1‑93 (MQ=255)

                                                           caTTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTcg                                >  1:6003333/1‑93 (MQ=255)

                                                            aTTCCAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTcgc                               >  1:6323940/1‑93 (MQ=255)

                                                               ccAGCAGACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCa                            <  1:6543843/93‑1 (MQ=255)

                                                                cagcagACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCAt                           <  2:3687873/93‑1 (MQ=255)

                                                                cagcagACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCAt                           <  2:104698/93‑1 (MQ=255)

                                                                   cagACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATAcc                        <  1:2588173/93‑1 (MQ=255)

                                                                   cagACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATAcc                        <  1:514078/93‑1 (MQ=255)

                                                                   cagACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATAcc                        >  1:3741772/1‑93 (MQ=255)

                                                                    agACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCa                       <  1:3453441/93‑1 (MQ=255)

                                                                    agACGAATAAAGATGCGTTTTTTGTCATAGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGGTCGCGCATACCa                       >  2:4286937/1‑93 (MQ=255)

                                                                     gACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGTCATGATTGAAATATCGAGTTCGCGCATACCAc                      >  2:4041083/1‑93 (MQ=255)

                                                                     gACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCAc                      <  1:472589/93‑1 (MQ=255)

                                                                     gACGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCAc                      >  1:4889748/1‑93 (MQ=255)

                                                                     gACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCAc                      >  1:6231022/1‑93 (MQ=255)

                                                                       cGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCAcgg                    >  2:5111775/1‑93 (MQ=255)

                                                                       cGAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGGCGTACTCATGTTCGGAAAACGGCATGATTGAAATATCGAGTTCGCCCATCACAcgg                    >  2:6291404/1‑93 (MQ=255)

                                                                        gAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTTGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCAcggc                   <  2:3046555/93‑1 (MQ=255)

                                                                        gAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCAcggc                   >  2:1784204/1‑93 (MQ=255)

                                                                        gAATAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCAcggc                   <  1:2092555/93‑1 (MQ=255)

                                                                         aaTAAAGATGCGTTTTTTGTCATCGGTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAACTATCGAGTGCGCGCATACCAcggcg                  >  2:6094045/1‑93 (MQ=255)

                                                                           tAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACggcggc                >  2:354763/1‑93 (MQ=255)

                                                                           tAAAGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACggcggc                <  2:1984881/93‑1 (MQ=255)

                                                                            aaaGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCa               <  2:575001/93‑1 (MQ=255)

                                                                             aaGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCAt              <  2:4439213/93‑1 (MQ=255)

                                                                             aaGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCAt              >  2:2732342/1‑93 (MQ=255)

                                                                              aGATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATg             <  2:3970633/93‑1 (MQ=255)

                                                                               gATGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGc            >  1:2577248/1‑93 (MQ=255)

                                                                                 tGCGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGccc          >  1:5116089/1‑93 (MQ=255)

                                                                                 tGCGTTGTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGccc          <  2:1292945/93‑1 (MQ=255)

                                                                                   cGTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCaa        <  2:2892964/93‑1 (MQ=255)

                                                                                    gTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCaat       >  1:5651042/1‑93 (MQ=255)

                                                                                    gTTTTTTGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCaat       >  2:2970867/1‑93 (MQ=255)

                                                                                     ttttttGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCaatg      >  1:2980007/1‑93 (MQ=255)

                                                                                     ttttttGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCaatg      >  1:4233498/1‑93 (MQ=255)

                                                                                       ttttGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCaatgaa    <  1:465663/93‑1 (MQ=255)

                                                                                         tCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCaatgaatg  <  1:2722132/93‑1 (MQ=255)

GTTCTGCATCGGCCGGTTTACCGTGATTCATCAGCCAGTTAAACAGGTCCGGATCGTCACATTCCAGCAGACGAATAAAGATGCGTTTTTCGTCATCGCTTAAGCTGTCGTACTCATGTTCGAAAAACGGCATGATTGAAATATCGAGTTCGCGCATACCACGGCGGCATGCCCAATGAATG  >  NC_000913/3038886‑3039067

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