BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,804,546	T→C	E50G (GAG→GGG)	rfaG ←	glucosyltransferase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,804,546	0	T	C	99.8%	1768.1 / ‑5.3	535	E50G (GAG→GGG)	rfaG	glucosyltransferase I
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (245/287); total (246/289)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

Only 100 of 545 total aligned reads displayed.

GGATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCTCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATA  >  NC_000913/3804456‑3804635

ggATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCa                                                                                         <  1:4359222/93‑1 (MQ=255)

ggATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTAAATAAGCCCa                                                                                         >  1:1830346/1‑93 (MQ=255)

 gATGCTCTTTGAGATGATTTTGTTCCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGTCACCTGAATAAGCCCaa                                                                                        <  2:7288555/93‑1 (MQ=255)

  aTGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCaaa                                                                                       <  2:3953043/93‑1 (MQ=255)

  aTGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCaaa                                                                                       <  2:1075393/93‑1 (MQ=255)

  aTGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCaaa                                                                                       >  2:7784352/1‑93 (MQ=255)

  aTGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCaaa                                                                                       <  2:8680676/93‑1 (MQ=255)

  aTGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCaaa                                                                                       <  2:9184826/93‑1 (MQ=255)

  aTGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCaaa                                                                                       >  1:4867443/1‑93 (MQ=255)

  aTGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCaaa                                                                                       >  1:3188692/1‑93 (MQ=255)

   tGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAAt                                                                                      >  1:6413070/1‑93 (MQ=255)

    gCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATg                                                                                     <  2:6585154/93‑1 (MQ=255)

      tctTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCt                                                                                   >  1:1497198/1‑93 (MQ=255)

       ctTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCtt                                                                                  >  1:1000796/1‑93 (MQ=255)

       ctTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCtt                                                                                  <  2:8870752/93‑1 (MQ=255)

        tttGATATGATGTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 <  1:7843701/93‑1 (MQ=255)

        tttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGGCTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 <  2:4664400/93‑1 (MQ=255)

        tttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 <  1:7774575/93‑1 (MQ=255)

        tttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 <  1:3013848/93‑1 (MQ=255)

        tttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 >  1:7008374/1‑93 (MQ=255)

        tttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 <  2:2542102/93‑1 (MQ=255)

        tttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 <  2:3604488/93‑1 (MQ=255)

        tttGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCttt                                                                                 <  2:3276051/93‑1 (MQ=255)

          tGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCg                                                                               <  2:9323383/93‑1 (MQ=255)

           gagaTGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCgg                                                                              <  1:1283387/93‑1 (MQ=255)

            agaTGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGATATGGCACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCggg                                                                             <  1:1750251/93‑1 (MQ=255)

              atgatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCa                                                                           <  2:5196824/93‑1 (MQ=255)

                gatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAAt                                                                         <  1:3997484/93‑1 (MQ=255)

                gatTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAAt                                                                         <  2:5963504/93‑1 (MQ=255)

                 atTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATc                                                                        >  1:9084674/1‑93 (MQ=255)

                  ttttGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCg                                                                       >  2:6365976/1‑93 (MQ=255)

                    ttGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGcc                                                                     >  1:1397685/1‑93 (MQ=255)

                      gTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCtt                                                                   <  2:6031354/93‑1 (MQ=255)

                      gTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCtt                                                                   <  2:4537219/93‑1 (MQ=255)

                       tACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTc                                                                  <  1:163310/93‑1 (MQ=255)

                            aGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACg                                                             >  2:7887093/1‑93 (MQ=255)

                             ggCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGa                                                            <  2:7617134/93‑1 (MQ=255)

                             ggCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGa                                                            <  1:231426/93‑1 (MQ=255)

                             ggCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGa                                                            >  2:8532601/1‑93 (MQ=255)

                             ggCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGa                                                            >  1:7517914/1‑93 (MQ=255)

                              gCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCGCCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGAc                                                           <  2:5031262/93‑1 (MQ=255)

                              gCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGAc                                                           <  2:82786/93‑1 (MQ=255)

                              gCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGAc                                                           <  1:8634699/93‑1 (MQ=255)

                                ataataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCGAATGCTTTCGGGCAATCGCCTTCCCACGACtg                                                         >  2:6953859/1‑93 (MQ=255)

                                ataataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtg                                                         <  2:1699180/93‑1 (MQ=255)

                                 taataTTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACtgt                                                        >  2:7506831/1‑93 (MQ=255)

                                   atatTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTgtgt                                                      >  2:9635061/1‑93 (MQ=255)

                                   atatTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTgtgt                                                      <  1:528434/93‑1 (MQ=255)

                                   atatTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTgtgt                                                      >  2:4331541/1‑93 (MQ=255)

                                      ttCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtata                                                   <  2:538546/93‑1 (MQ=255)

                                      ttCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtata                                                   <  1:2664196/93‑1 (MQ=255)

                                       tCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTatatat                                                  >  1:9233805/1‑93 (MQ=255)

                                        cTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtatata                                                 >  1:3650297/1‑93 (MQ=255)

                                        cTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtatata                                                 <  1:5593593/93‑1 (MQ=255)

                                        cTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGtatata                                                 <  2:2400232/93‑1 (MQ=255)

                                         tGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATAc                                                <  1:8182716/93‑1 (MQ=255)

                                          gCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATAcc                                               >  1:6502903/1‑93 (MQ=255)

                                           cATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATAccc                                              >  1:6820232/1‑93 (MQ=255)

                                            aTTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCg                                             >  2:3472905/1‑93 (MQ=255)

                                             ttGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTTCCACGACTGTGTATATACCCGa                                            >  1:7736088/1‑93 (MQ=255)

                                               gCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAAc                                          <  2:3002126/93‑1 (MQ=255)

                                                cGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACa                                         >  2:1341360/1‑93 (MQ=255)

                                                 gTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAt                                        >  2:4129768/1‑93 (MQ=255)

                                                 gTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAt                                        <  1:7125733/93‑1 (MQ=255)

                                                 gTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAt                                        <  2:5715331/93‑1 (MQ=255)

                                                 gTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGCGTATATACCCGAACAt                                        <  2:9556923/93‑1 (MQ=255)

                                                  tCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtg                                       >  2:3895481/1‑93 (MQ=255)

                                                   ccATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtgg                                      >  1:2568539/1‑93 (MQ=255)

                                                   ccATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtgg                                      >  1:9530866/1‑93 (MQ=255)

                                                     aTGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtg                                    >  1:579091/1‑93 (MQ=255)

                                                     aTGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtg                                    <  2:513599/93‑1 (MQ=255)

                                                      tggttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtgg                                   >  2:9812063/1‑93 (MQ=255)

                                                      tggttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACAtggtgg                                   <  1:10347/93‑1 (MQ=255)

                                                         ttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGccc                                <  1:8415252/93‑1 (MQ=255)

                                                         ttggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGccc                                >  1:9930876/1‑93 (MQ=255)

                                                           ggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCg                              >  1:9706944/1‑93 (MQ=255)

                                                           ggtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCg                              <  1:2834796/93‑1 (MQ=255)

                                                            gtATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGt                             >  2:2754440/1‑93 (MQ=255)

                                                             tATGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTg                            >  1:2022025/1‑93 (MQ=255)

                                                              aTGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGc                           <  1:8891864/93‑1 (MQ=255)

                                                              aTGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGc                           <  1:7360860/93‑1 (MQ=255)

                                                              aTGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGc                           <  1:5112358/93‑1 (MQ=255)

                                                               tGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCg                          >  2:7954364/1‑93 (MQ=255)

                                                               tGGGACTTAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCg                          <  1:8030504/93‑1 (MQ=255)

                                                                     ttAACTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACt                    >  2:9446067/1‑93 (MQ=255)

                                                                       aaCTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGt                  <  1:9349144/93‑1 (MQ=255)

                                                                         cTGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTg                <  2:44751/93‑1 (MQ=255)

                                                                          tGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGa               <  2:2122316/93‑1 (MQ=255)

                                                                          tGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGa               <  2:7037376/93‑1 (MQ=255)

                                                                          tGGCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGa               <  2:4812601/93‑1 (MQ=255)

                                                                            gCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATg             <  2:4896100/93‑1 (MQ=255)

                                                                            gCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATg             <  1:4608561/93‑1 (MQ=255)

                                                                            gCACCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATg             <  2:1916277/93‑1 (MQ=255)

                                                                              aCCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCa           <  1:5904376/93‑1 (MQ=255)

                                                                              aCCTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCa           >  2:6211057/1‑93 (MQ=255)

                                                                                cTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAAt         <  2:8117069/93‑1 (MQ=255)

                                                                                cTGAATAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAAt         <  2:9870488/93‑1 (MQ=255)

                                                                                   aaTAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATgcg      >  1:6074252/1‑93 (MQ=255)

                                                                                    aTAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCCATgcgc     >  2:5503299/1‑93 (MQ=255)

                                                                                     tAAGCCCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCa    <  1:335177/93‑1 (MQ=255)

GGATGCTCTTTGAGATGATTTTGTACCCAGGCATAATATTCTGCATTGCGTCCATGGTTGGTATGGGACTTAACTGGCACCTGAATAAGCTCAAATGCTTTCGGGCAATCGCCTTCCCACGACTGTGTATATACCCGAACATGGTGGCCCCGTGCGGCAACTGTTGATGCAATGCGCATA  >  NC_000913/3804456‑3804635

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